|
|
||||||||
1 Department of Biochemistry and Molecular Biology, University of Miami School of Medicine, Miami, Florida 33101, USA
2 Department of Biology and Biochemistry, University of Bath, Bath BA2 7AY, United Kingdom
Reprint requests to: K. Brew, Department of Biomedical Sciences, Florida Atlantic University, 777 Glades Road, P.O. Box 3091, Boca Raton, Florida 33431, USA; e-mail: kbrew{at}fau.edu; fax: (561) 297-2221.
(RECEIVED March 12, 2001; FINAL REVISION August 14, 2001; ACCEPTED August 16, 2001)
3 Present address: L.H. Greene, University of Oxford-Oxford Centre for Molecular Sciences, New Chemistry Laboratory, South Parks Road, Oxford OX1 3QT, United Kingdom. ![]()
4 Present address: A.C. Papageorgiou, Turku Centre for Biotechnology, University of Turku and Åbo Akademi University, BioCity, Turku 20521, Finland. ![]()
5 Present address: K. Brew, Department of Biomedical Sciences, Florida Atlantic University, Boca Raton, Florida 33431, USA. ![]()
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1101/ps.22901.
| Abstract |
|---|
|
|
|---|
Keywords: Protein stability; mutagenesis; evolution; conserved residues; lipocalins; protein folding, serum retinol-binding protein; crystal structures
Abbreviations: RBP, retinol-binding protein rRBP, recombinant human retinol-binding protein rRBPmX, single site mutants (indicating the substitution of amino acid X for residue number m) rRBPmX/nY, multi-site mutants (indicating the substitution of X for residue m and Y for residue n) MG, molten globule state PCR, polymerase chain reaction GndHCl, guanidine hydrochloride
| Introduction |
|---|
|
|
|---|
-helix, formed by a region close to the C terminus (Fig. 1
|
The folding mechanisms of ß-sheet proteins and, in particular, ß-barrel proteins are less understood than those of other structural classes. In particular, the conversion of an unstructured polypeptide to form an open barrel would require major structural reorganization after hydrophobic collapse (Clark et al. 1996). The structures of partially folded states may be relevant to this process, and in RBP a molten globule (MG) state (Ohgushi and Wada 1983) is stabilized at low pH (Bychkova et al. 1992); the MG of another lipocalin, ß-lactoglobulin, is reported to be identical with an early intermediate in the refolding process (Hamada et al. 1996). In ß-lactoglobulin it appears that nonnative structures form during folding because the MG has higher helix content than the native structure (Ragona et al. 1997). Despite their low global levels of sequence identity, different lipocalins show remarkable conservation of five patches of sequence that are associated with a structurally conserved mixed nonpolar/polar core at the base of the ß-barrel. This region appears to be a key structural element of the lipocalin fold and has been suggested to play a pivotal role in folding (Greene and Brew 1995; Brew and Greene 1997; Greene 1998).
We describe here initial studies of the role of the four tryptophans of RBP and two additional residues that are adjacent to the highly conserved Trp24 (see Fig. 1
) in the structure and stability of the apo-protein. Expression systems have been previously reported for apo-forms of human and porcine rRBP (Muller and Skerra 1993; Sivaprasadarao and Findlay 1993; Wang et al. 1993), but the recombinant proteins were not subject to detailed structural and biophysical characterization. We present the high-resolution crystal structure of recombinant RBP, folded from inclusion bodies by a new process, and a variant with substitutions for the surface tryptophans 67 and 91. Spectroscopic and stability measurements are reported for human apo-rRBP and 10 mutants. The properties of these proteins idey unique contributions made by different tryptophans to the CD and fluorescence spectrum of the native protein, indicating approaches and spectroscopic signals to investigate the development and loss of structure during folding and unfolding processes. The effects of mutations underscore the uniquely important role of interactions involving Trp24 in molecular stability.
| Results |
|---|
|
|
|---|
The structures of rRBP and rRBP67L/91H were determined at 1.7 Å and 2.0 Å resolution, respectively. Although the complete amino acid sequence of RBP comprises 182 amino acids (Rask et al. 1979), clear electron density was visible only for 174 out of 182 residues in both structures. The overall structures of the two proteins are very similar (root mean square deviation, 0.32 Å) apart from the flexible loop formed by residues 6268 and minor conformational changes of some surface residues (Fig. 2
).
|
|
Flexible loop (residues 6268)
The loop region 6268 in the rRBP structure was highly disordered, and the incorporation of the excluded residues, which were represented as alanines in the model, was possible only after applying the automated refinement procedure of Lamzin and Wilson (1993). Subsequent cycles of refinement considerably improved the quality of the electron density map, and the amino acids forming the flexible loop were introduced into the model, according to the sequence, only toward the final stage of refinement. The rRBP67L/91H structure had features similar to rRBP, except the loop region is better defined in the electron density.
This flexible loop and, in particular, Trp 67 is known to be involved in molecular packing interactions at the interface of the human RBP-TTR (transthyretin) complex (Monaco et al. 1995; Naylor and Newcomer 1999). In the two structures reported for this complex, Trp67 and Trp91 have critical roles in heterodimer stabilization, but a more detailed examination of these particular residues is not possible because of the low resolution (3.1 Å). Thus, when the rRBP67L/91H structure is compared with the structure of RBP in the complex, only gross structural differences can be ideied. The most profound differences between the recombinant apo-RBP and RBP-TTR complex are those around residue 62 and at the C terminus (root mean square, 1.05 Å; Fig. 2
); both regions are implicated in interactions with TTR in the complex. In the rRBP structure, Trp67 is disordered, as are the rest of the residues that form the flexible loop, whereas in the rRBP67L/91H structure, the effect of the sequence substitution at position 67 was not investigated in detail because of the poor electron density in this region.
ß-Barrel structure
The hydrophobic ß-barrel structure is very similar to those in previously reported RBP structures (Figs. 1, 2![]()
). The structural similarity of RBP and ß-lactoglobulin, which is isolated from milk as an apo-protein (Sawyer et al. 1985), also shows that the absence of a ligand does not cause a collapse of the ß-barrel. The structures of rRBP and rRBP67L/91H and the apo-form of RBP obtained by depleting the natural holo-protein of retinol (Zanotti et al. 1993c) confirm this suggestion. The barrel of RBP accommodates one retinol molecule (Cowan et al. 1990), located in the interior with the ß-ionone ring pointing to the deepest region of the barrel and the isoprene tail stretching almost to the barrel entrance. The apo-form of RBP derived from protein obtained from human plasma was determined at 2.5 Å resolution by Zanotti et al. (1993c), who observed that in the absence of retinol, the internal cavity contained significant portions of FoFc electron density. Molecules responsible for this density could not be ideied because of insufficient resolution of the X-ray data, and it was suggested that the unexplained density might be attributed to solvent molecules. The structure of rRBP at the higher resolution (1.7 Å) presented here clarifies the tentative conclusions of Zanotti and coworkers by revealing the presence of both water and what appear to be glycerol molecules in the ß-barrel in the rRBP and rRBP67L/91H structures. Glycerol molecules were incorporated in the model only during the final stage of refinement and are well positioned in the density, being stabilized by hydrogen bonds and van der Waals interaction with residues lining the interior of the ß-barrel.
The hydrophobic nature of the interior of the ß-barrel is evident from the presence of a large number of aromatic rings located within the lumen of the barrel. Phe36, which is located at the entrance to the interior, appears to block the opening of the barrel in the apo-RBP structure (Zanotti et al. 1993c) but adopts a different conformation when retinol is bound that opens the entrance of the barrel (Cowan et al. 1990; Naylor and Newcomer 1999). Thus, the side-chain of Phe36 seems to act like a tollgate for retinol binding at the entrance of the barrel. Although a similar feature is observed in the structures of rRBP and rRBP67L/91H (Fig. 3
), Phe36 does not prevent the binding of small organic molecules such as glycerol in the barrel. Interestingly, the position of two glycerol molecules in the barrel in the rRBP and the variant structures appear to mimic the interactions made by the retinol molecule in holo-RBP structure.
Environments of tryptophans within the three-dimensional structure and conservation in sequence
The four tryptophans vary in their degree of conservation in a comparison with 32 members of the lipocalin superfamily (Fig. 4
), with Trp24 and the proximal Gly22 (numbering in accordance with recombinant human RBP) being the most highly conserved residues. Lipocalins differ in function and structure-function relationships, and the conservation of Trp24 does not reflect a conserved role in function because it does not directly interact with bound ligands in currently known structures. Two other tryptophans, Trp67 and Trp91, are not conserved but Trp105 is conserved in all 12 known RBP sequences (data not shown).
|
|
|
N of Lys85 in the present rRBP structure is
3.5 Å from the center of the indole ring of Trp105 (measured using Hyperchem), reflecting their involvement in an amine-aromatic electrostatic interaction.
Design and characterization of mutants
Ten single and multisite mutants were expressed and characterized to investigate the contributions of the different tryptophans to structure, stability, and spectroscopic properties of RBP. Therefore, the substitutions were designed to introduce minimal changes in structure and stability using, where possible, amino acids at corresponding sites in closely related lipocalins or to introduce a structurally similar amino acid. Substitutions of Leu for Trp67 and of His for Trp91 reflect a sequence change in Xenopus RBP, and Phe for Trp105 was based on structural conservation. Trp24 was replaced by Tyr, based on the previously discussed substitution in late lactation proteins from marsupials (Collet et al. 1989) and by Phe and Leu to introduce a weakly fluorescent or nonfluorescent residue. These changes will also eliminate the H-bond between the indole nitrogen and the carbonyl oxygen of Phe20. Gln was substituted for Arg139 as an alternative method of disrupting the Trp24-Arg139 interaction (Table 2
). A mutant with substitutions of Phe for Trp24 and of Ala for the highly conserved Gly22 was also produced by a polymerase chain reaction (PCR)-derived mutation during the construction of the Trp24-to-Phe mutant. Designed multisite mutants were rRBP67L/91H, rRBP67L/91H/105F, and rRBP24Y/67L/91H. All rRBP mutants gave single bands of similar size to wild-type RBP on sodium dodecyl sulfat (SDS) polyacrylamide gel electrophoresis. Retinol-binding studies with selected tryptophan mutants indicated that the mutations did not eliminate the ability of the protein to bind all-trans retinol (data not shown). After folding and purification, the mutants were isolated in yields of
20 mg/L, except for those with substitutions for Trp24 or Arg139, in which the yields were 4- to 20-fold lower. To obtain small quantities of pure forms of rRBP24L and rRBP139Q, it was necessary to perform an additional separation by fast protein liquid anion exchange chromatography.
Contributions of tryptophans to the CD spectrum
Near-ultraviolet (UV) CD spectra of proteins contain features that reflect the fixed asymmetric environments of aromatic amino acids and cystines in the tertiary structure (Strickland 1974). Far-UV CD spectra largely arise from secondary structure, but aromatic residues can also contribute to this region (Woody 1995). Figure 5A and 5D![]()
shows near- and far-UV CD spectra of rRBP under native conditions (pH 7.4 and 20°C) and in acid-denatured (pH 2.0) and solvent-denatured (6 M guanidine hydrochloride [GndHCl]) states. The near- and far-UV CD spectra of rRBP in the native state are very similar to those previously reported by Bychkova et al. (1992) for human serum apo-RBP. rRBP unfolds at low pH and elevated temperatures to produce partially denatured conformers that have no near-UV CD spectrum and characteristic far-UV CD spectra that differ from the spectrum of the fully denatured protein. Both acid-denatured and thermally denatured rRBP (spectrum not shown) have greater negative ellipticity between 210 and 240 nm than the native protein, and the spectrum at pH 2.0 is similar to that previously described as a MG by Bychkova et al. (1992).
|
-helix and ß-sheet (Table 3
-helix formation in the primary structure in comparison with the helix content present in the native structure determined by X-ray crystallography. However, the helix content of RBP predicted using several sequence-based secondary structure prediction algorithms is in good agreement with the secondary structure content of the native structure (data not shown), indicating that the formation of MG with enhanced helix content may be a feature of the folding process in both lipocalins.
|
280 nm (Fig. 5B
275 nm (Fig. 5C
228 nm (Fig. 5F
228 nm (Fig. 5E,F). This region of the far-UV CD spectra is only affected by substitutions for Trp24 or Trp105, and the ellipticity at 228 nm is therefore a signal that can be used to monitor the tertiary structure interactions of Trp24 and Trp105. Substitutions for Trp24 do not affect the trough at
208 nm (Fig. 5F
Roles of different tryptophans in fluorescence spectra of native and unfolded rRBP
RBP contains 4 tryptophan, 8 tyrosine, and 10 phenylalanine residues. The fluorescence emission spectra of rRBP and mutants missing tryptophan 24 (rRBP24F), tryptophan 105 (rRBP105F), or tryptophans 67, 91, and 105 (rRBP67L/91H/105F) in native and unfolded states when excited at 295 nm are shown in Figure 6
. Wild-type rRBP shows a red shift in the emission maximum from 345 to 355 nm and a 20% decrease in fluorescence intensity on unfolding with 6 M GndHCl. Emission maxima for the denatured states of most mutants were similar (352 to 357 nm), but those of mutants with substitutions of Trp24 by Phe or Tyr (data not shown) were shifted to lower wavelength (351 and 349 nm). The blue-shifted contributions of Trp 24 and Trp 105 mutants in the native state fluorescence spectra compared with the denatured state (Fig. 6
) are consistent with their low levels of solvent-accessibility (Table 1
) and location in a hydrophobic environment (Table 1
). The tryptophan fluorescence spectra indicate that Trp 24 makes the largest contribution to the decrease of fluorescence on unfolding, which may be partly attributable to a change in polarity of its environment on folding (Burnstein et al. 1973). However, this tryptophan is also part of a cluster of aromatic residues, its indole ring being 6 to 7 Å from the side-chains of tyrosines 25, 111, and 114 in the structure of rRBP, and it is also possible that the enhanced fluorescence on folding associated with Trp24 arises from fluorescence energy transfer from the adjacent tyrosine side-chains. These results indicate that Trp24 would serve as a good fluorescent probe for monitoring the formation of the core during studies of folding and unfolding kinetics.
|
|
|
GH2O) and at Gnd50 (Table 4b
|
| Discussion |
|---|
|
|
|---|
Although all rRBP mutants were expressed in similar yields in Escherichia coli as inclusion bodies, after in vitro folding the native forms of all Trp24 mutants were isolated in lower yields (
25%) compared with those of the wild-type protein and other mutants; in the case of the Trp24Leu and Arg139Gln variants, the yields were about fourfold lower. This suggests that these mutations increase the stability of misfolded forms relative to the native structure so that a larger proportion of protein molecules fold to nonnative structures. This view is reinforced by the lack of secretion of the Trp19 to Ala and Phe mutants of ß-lactoglobulin by the eukaryotic expression system because such systems tend to degrade wrongly folded proteins. These findings are conceptually significant because they link mutations in conserved residues to protein misfolding, which is the basis of many diseases (Dobson 2001).
The folding behavior of RBP has a specific biological interest because only the holo form of the protein is secreted after biosynthesis in mammalian and other cells; RBP molecules that do not acquire a retinol ligand within the cell are retained in the endoplasmic reticulum (Melhus et al. 1992) and appear not to be fully folded (Kaji and Lodish 1993). Previously determined structures for human and bovine apo-RBP show close similarity to the holo-protein. These apo-proteins were prepared from natural holo-protein after extraction with ethyl ether to remove the bound ligand (Zanotti et al. 1993a). Here we find that recombinant human apo-RBP produced by in vitro folding of material extracted from inclusion bodies has a structure and spectroscopic properties that are closely similar to those of apo and holo forms of natural human holo-RBP (Cowan et al. 1990), holo and apo bovine RBP (Zanotti et al. 1993a,c), and RBP in complexes with different retinoids (Zanotti et al. 1993b). Because our preparations of apo-rRBP have never bound a retinol ligand, we can conclude that although retinol may enhance folding yields, the ligand is not necessary for an irreversible maturation step in folding. Thus, the degradation of apo-RBP in vivo must be linked to some specific structural feature or property of the apo- versus holo-protein.
The results described here are necessary for the design and interpretation of unfolding and folding kinetics of rRBP because they establish conditions in which the folded and unfolded conformers are most populated. They also allow comparisons of the effects of mutations on the stability of the transition states for folding and unfolding with those of the native and unfolded states. The spectroscopic properties of the mutants indicate signals that provide information about the structure formation in different parts of the RBP molecule during folding processes. A major focus of our work is to determine if there is a relationship between folding and sequence conservation in functionally divergent paralogous proteins. This is a concept that is recently receiving attention from experimentalists (Martinez and Serrano 1999; Hamill et al. 2000; Nishimura et al. 2000; Plaxco et al. 2000). rRBP, as an experimental model for the large lipocalin superfamily, represents an ideal system to further our understanding between evolution and folding. Toward this end, kinetic and additional X-ray crystallographic studies are in progress with recombinant RBP and mutants.
| Materials and methods |
|---|
|
|
|---|
and transformed into BL21(DE3) cells for expression. DNA sequencing of a selected clone showed that it is identical to that published by Colantuoni and coworkers (1983) apart from two silent substitutions, TCC for TCG in the codon for serine 46 and GTG for GTA in the codon for valine 116.
Construction of RBP mutants
Mutations were introduced using the "megaprimer" method (Sarkar and Sommer 1990). The megaprimer was generated by amplification with T7 promoter primer or T7 terminator primer together with an appropriate mutagenic primer (Grobler et al. 1994). The megaprimer was purified by agarose gel electrophoresis and used in the second amplification with the same template and the cognate T7 primer. The final amplification product was purified and cloned into pET3a as described for RBP. The complete sequences of the inserts in the recombinant vectors were determined by automated DNA sequencing using both the T7 promoter and T7 terminator primers. Oligonucleotide synthesis and automated DNA sequencing were performed in the DNA Core Facility, University of Miami School of Medicine.
Protein expression and initial purification of rRBP
Cultures of E. coli strain BL21(DE3) transformed with the recombinant vector (pRBP) or variant were grown at 37°C with 100 mg/mL ampicillin, induced with 0.4 mM isopropylthiogalactopyranoside (IPTG) when the optical density reached 0.6, and harvested after 3 h. Cells were suspended in lysis buffer (0.05 M Tris, 0.1 M NaCl, and 1 mM EDTA at pH 8.3) and lysed by treatment with lysozyme and deoxycholic acid and treated with DNAse (Sambrook et al. 1989). Inclusion bodies containing recombinant RBP (rRBP) were isolated by centrifugation and washed initially by resuspension in lysis buffer and then with double distilled water. The protein was extracted with 8 M urea containing 20 mM Tris base and 20 mM dithiothreitol (
50 mL per 6 L of culture), centrifuged at 15,000 rpm for 10 min, and loaded on a column (2 x 10 cm) of Macroprep 50Q equilibrated with 20 mM Tris-HCl (pH 8.5) containing 8 M urea. The column was washed with equilibration buffer until the absorbance at 280 nm was nearly to baseline, and the protein was eluted with a linear gradient from 0 to 0.5 M NaCl in 20 mM Tris-HCl (pH 8.5) containing 8 M urea at a flow rate of 2 mL/min. Fractions of 6.0 mL were collected and monitored by absorbance at 280 nm. SDS gel electrophoresis of aliquots of different fractions indicated that the major UV-absorbing peak, which eluted at
25 min, contained rRBP (molecular weight 21,000) in a relatively pure form.
In vitro folding and purification of folded rRBP
Fractions containing rRBP or a mutant were pooled, and the concentration of protein was estimated from the absorbance at 280 nm; the extinction coefficients for RBP and tryptophan mutants were calculated from their contents of tryptophan, tyrosine, and cysteine (Mach et al. 1992). E280nm0.1% values are as follows: wild-type rRBP and rRBP139Q, 1.61; rRBP24Y, 1.44; rRBP24F, rRBP24L, rRBP91H, rRBP105F, and rRBP22A/24F, 1.34; rRBP67L/91H, 1.06; rRBP67L/91H/105F, 0.79; and rRBP24Y/67L/105F, 0.80. Native protein was generated by adding the denatured material in 8 M urea drop-wise to folding buffer (20 mM Tris HCl at pH 8.5 containing 10 mM ß-mercaptoethanol, 1 mM 2-hydroxyethyldisulfide, and 1% glycerol) at a rate of
30 drops/min. The final concentrations of protein and urea were 50 to 100 µg/mL and 0.3 M, respectively. Folding was allowed to proceed for 16 h at 4 °C, and the solution was then concentrated using an Amicon Spiral Filter Ultrafiltration System. The protein was loaded on a column (2 x 10 cm) of DE52 DEAE cellulose equilibrated with 20 mM Tris-HCl (pH 8.5) at 22°C. After washing with 20 mM Tris HCl (pH 8.5), the sample was eluted with a linear gradient from zero to 0.5 M NaCl in the same buffer over 85 min at a flow rate of 2 mL/min. A single peak containing folded protein eluted at
0.3 M NaCl and was dialyzed against water at 4°C and lyophilized. Further purification was performed by gel filtration with a column of Superdex G75 (7 x 21cm) equilibrated with 20 mM Tris-HCl (pH 7.4) and 0.2 M NaCl at a flow rate of 1.5 mL/min. CD spectroscopy showed that the second, larger protein peak contained native protein.
Folded rRBP was obtained by this method in a yield of 20 mg/L of bacterial culture. The purified protein was homogeneous on SDS polyacrylamide gel electrophoresis. Fourteen cycles of Edman degradation using an Applied Biosystems model 470A Protein/Peptide Sequencer gave a single N-terminal sequence identical to that of human RBP, except for the presence of an N-terminal methionyl extension that reflects its intracellular expression in E. coli. The ligand-binding properties of RBP were measured using retinoic acid, the interaction of which with RBP quenches tryptophan fluorescence. The calculated dissociation constant (Kd) for rRBP was 1.9 ± 0.2 x 10-7 M (data not shown), which is within experimental error of the value of 2.1 x 10-7 M reported for natural human apo-RBP (Cogan et al. 1976).
X-ray crystallographic analysis
Crystallization rRBP
Large needle shaped crystals of native recombinant human RBP were grown by the vapor diffusion method at 16°C. Equal volumes (2.0 µL) of protein (
12 mg/mL in 10 mM sodium cacodylate and 3.9 M sodium chloride at pH 6.8) and reservoir solution (20 mM sodium cacodylate and 4.5 M sodium chloride at pH 6.8) were mixed on siliconized coverslips and left to equilibrate against the reservoir solution.
rRBP67L/91H variant
rRBP double mutant was crystallized by the vapor diffusion method at 16°C. Single crystals were obtained by applying the seeding technique from crystals of the native recombinant protein, and their size was optimized by addition of 0.05% ß-octyl-glucoside. Equal volumes (2.0 µL) of protein (
12 mg/mL) and reservoir solution (4.5 M sodium chloride, 50 mM Tris-HCl at pH 9.1, and 0.05% ß-octyl-glucoside) were used.
Diffraction data collection-rRBP
Data were collected to 1.7 Å resolution from one crystal at the Synchrotron Radiation Source, Daresbury on station PX 9.5 (oscillation range, 2.0°;
= 0.9 Å) under cryogenic conditions (100 K) using a 30-cm MAR Research image plate. Before data collection, crystals were transferred to buffer solution (same as the reservoir solution) containing 25% (v/v) glycerol (cryoprotectant) for
30 sec. The data set was indexed in rhombohedral spacegroup R3, with one molecule in the asymmetric unit and solvent content of 68.3% (Matthews coefficient of
3.9 Å3/Da). Data integration and reduction were performed with the programs DENZO and SCALEPACK (Otwinowski and Minor 1997). The data collection statistics are presented in Table 5
.
|
= 1.09 Å). The data set was collected under cryogenic conditions (100 K) using 20% (v/v) glycerol as cryoprotectant. Both DENZO and SCALEPACK programs were used for data processing (Otwinowski and Minor 1997). The double variant was crystallized in rhombohedral spacegroup R3, like the native protein, with one molecule in the asymmetric unit and solvent content of 69.5% (Matthews coefficient of
4.0 Å3/Da; Table 5
Structure determination rRBP
The structure of human serum RBP (hRBP) determined previously at 2.0 Å resolution (Cowan et al. 1990) was used for molecular replacement with the program AMoRe (Navaza 1994). Superposition of hRBP with rat androgen-dependent epididymal retinoic acidbinding protein (E-RABP; Newcomer 1995) indicated that residues 1430, 70, 80, 8486, 106115, 135139, 146151, 158160, 165, 172, and 174 are conserved in the RBP structure; hence, only these were included in the search model. Residues that were part of flexible loops or other disordered regions were excluded from the starting model.
rRBP67L/91H variant
The structure of the double variant at 2.0 Å resolution was determined by molecular replacement with X-PLOR (Brünger 1992b) using the rRBP structure determined at 1.7 Å resolution as a search model.
Structure refinement
The resultant models of rRBP and rRBP67L/91H were first subjected to restrained refinement with maximum likelihood method using the program REFMAC (CCP4 1994) followed by the automated refinement procedure (ARP; Lamzin and Wilson 1993) and to rigid body refinement using the program X-PLOR, respectively. The combination of REFMAC and ARP refinement protocol, applied to the starting model of rRBP, significantly improved the quality of the electron density map and enabled the building of the residues that had been initially excluded from the molecular replacement solution (only conserved residues). Cycles of refinement were performed using the slowcool protocol, at first with the program X-PLOR (Brünger 1992b) and during the final stages of refinement with CNS (Brünger et al. 1998), for both structures. The progress of refinement was monitored through both free and conventional R-factors (Brünger 1992a). Alternate cycles of manual rebuilding with the graphics program O (Jones et al. 1991) and refinement cycles using the standard protocol improved the quality of the model. Simulated annealing omit maps calculated using either X-PLOR or CNS for the two structures indicated density in the core of the ß-barrel and on the surface sufficient to accommodate one or more glycerol molecules. Six glycerol molecules were incorporated into the rRBP model (five in the rRBP67L/91H model) at the very final stage of refinement. In each case, water molecules were inserted in the model only if there were peaks in the FoFc electron density maps with heights >3
and they were at hydrogen bondforming distances from appropriate atoms. 2FoFc electron density maps were also used to check the consistency in peaks. The final refinement statistics and model quality are summarized in Table 5
. Analysis of the Ramachandran plot, calculated with the program PROCHECK (Laskowski et al. 1993), showed that all residues lay in the allowed regions except Tyr 111 in both the rRBP and rRBP67L/91H structures, for which the
and
values are
= 69.1°,
= `39.8° and
= 64.6°,
= -36.7°, respectively. The details of refinement statistics and model accuracy are listed in Table 5
. Structural superpositions were performed with the program MAPS (Lu 2000) and the accessibility of the protein surface was calculated with the program DSSP (Kabsch and Sanders 1983). The atomic coordinates have been deposited at the RCSB Protein Data Bank (accession codes 1JYD for apo-rRBP and 1JYJ for rRBP67L/91H).
Circular dichroism
Near- and far-UV CD spectra of rRBP and mutants were determined with a Jasco J-710/720 spectropolarimeter. Twelve spectra were scanned for each sample at a speed of 50 nm/min and were averaged. Near-UV CD spectra (240 to 320 nm) were determined using a cell with a path length of 1.0 cm; far-UV spectra (200 to 240 nm), using a path length of 0.1 cm. Protein samples were dissolved in 5 mM sodium phosphate buffer (pH 7.4) at concentrations of 0.2 to 0.5 mg/mL. A constant temperature of 20°C was maintained with a circulating NESLab RTE-111 water bath. MG spectra were determined with protein samples dissolved in 50 mM phosphate buffer (pH 2.0). Secondary structure analysis of far-UV CD spectra was performed using the k2d neural network program (Andrande et al. 1993).
Thermal unfolding
Thermal unfolding data were analyzed by fitting to the following equation:
![]() |
Hm is the enthalpy of unfolding, Tm is the midpoint of the unfolding transition, yn is the ellipticity of the native protein and yd is the ellipticity of the unfolded protein (both extrapolated to 0 K), and mn and md are the baseline slopes for the native and denatured states, respectively. The approximate effects of mutations on stability were calculated from the changes in Tm (
Tm) produced by the following mutation:
![]() |

G(WT-mutant) is the stability difference between wild-type protein and mutant, Tm(WT) and Tm(mutant) are the mid-points of the unfolding transition of the wild-type protein and mutant, and
Sm(WT) is entropy of unfolding of the wild-type protein.
Sm is calculated from the relationship
Sm =
Hm/Tm.
Chemical-induced equilibrium unfolding
Fluorescence spectra of native and unfolded forms of protein and equilibrium GndHCl unfolding were determined using a Perkin Elmer Luminescence LS-50 Fluorimeter at 23°C. Protein samples (50 to 130 µg/mL) in 5 mM phosphate buffer (pH 7.0) containing 0 to 6 M GndHCl were preincubated for 24 h at 23°C, and fluorescence spectra were measured from 300 to 400 nm using an excitation wavelength of 295 nm. For unfolding measurements at different denaturant concentrations, samples were excited at 290 nm, and the average fluorescence emission between 329 and 331 nm was determined. An excitation wavelength of 290 nm was used in the latter experiments to maximize the fluorescence yield and difference between native and unfolded states. Each data point is the average of nine scans. The concentrations of GndHCl in sample solutions were determined by refractive index at 25°C (Nozaki 1972) using a Milton Roy ABBE-3L refractometer. In all unfolding studies, solutions were prepared using Hamilton syringes to ensure accuracy.
The fluorescence spectra were analyzed by fitting to appropriate equations using the curve fitter program Sigma Plot, version 3 (Jandel Scieic). Baseline data from the linear native and unfolded regions of the curve (ideied from initial plots of fluorescence against denaturant) were analyzed by linear regression to obtain baseline corrections reflecting the response of the native and unfolded states to denaturant. Subtraction of these across the unfolding curve allowed the calculation of the fraction of folded and unfolded protein across the transition (Pace 1986). Free energies of stability were calculatedas
G = -RTln(Kf), where Kf = fraction folded/fraction unfolded. Linear plots of
G[denaturant] versus [denaturant] allowed the calculation of
G(water) and [denaturant]50%the free energy of stability in the absence of denaturant and the denaturant concentration at the midpoint of the unfolding transition, respectivelytogether with their standard errors.
| Acknowledgments |
|---|
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
Beg, O.U. and Shaw, D.C. 1994. The complete primary structure of late lactation protein from quokka (Setonix brachyurus). J. Protein Chem. 13: 513516.[CrossRef][Medline]
Beste, G., Schmidt, F.S., Stibora, T., and Skerra, A. 1999. Small antibody-like proteins with prescribed ligand specificities derived from the lipocalin fold. Proc. Natl. Acad. Sci. 96: 18981903
Bianchet, M.A., Bains, G., Pelosi, P., Pevsner, J., Snyder, S.H., Monaco, H.L., and Amzel, L.M. 1996. The three-dimensional structure of bovine odorant binding protein and its mechanism of odor recognition. Nat. Struct. Biol. 3: 934939.[CrossRef][Medline]
Bishop, R.E. and Weiner, J.H. 1996. Outlier lipocalins more than peripheral. Trends Biochem. Sci. 21: 127.[CrossRef][Medline]
Bocskei, Z., Groom, C.R., Flower, D.R., Wright, C.E., Phillips, E.V., Cavaggioni, A., Findlay, J.B.C., and North, A.C.T. 1992. Pheromone binding to two rodent urinary proteins revealed by X-ray crystallography. Nature 360: 186188.[CrossRef][Medline]
Brew, K. and Greene, L. 1997. Evolution, folding and flexibility. Protein Eng. 10: S44.
Brownlow, S., Morais Cabral, J.H., Cooper, R., Flower, D.R., Yewdall, S.J., Polikarpov, I., North, A.C.T., and Sawyer, L. 1997. Bovine ß-lactoglobulin at 1.8 Å resolution: Still an enigmatic lipocalin. Structure 5: 481495.[Medline]
Brünger, A.T. 1992a. Free R value: A novel statistical quantity for assessing the accuracy of crystal structures. Nature 355: 472475[CrossRef]
Brünger, A.T. 1992b. X-PLOR Version 3.1 Manual: A system for X-ray Crystallography and NMR, Yale University Press, New Haven.
Brünger, A.T., Adams, P.D., Clore, G.M., DeLano, W.L., Gros, P., GrosseKunstleve, R.W., Jiang, J.S., Kuszewski, J., Nilges, M., Pannu, N.S., Read, R.J., Rice, L.M., Simonson, T., and Warren, G.L. 1998. Crystallography and NMR system: A new software suite for macromolecular structure determination. Acta Crystallogr. D 54: 905921.[CrossRef][Medline]
Bugos, R.C. and Yamamoto, H.Y. 1996. Molecular cloning of violaxanthin de-epoxidase from romaine lettuce and expression in Eschericia coli. Proc. Natl. Acad. Sci. 93: 63206325.
Burnstein, E.A., Vedenkina, N.S., and Ivkova, M.N. 1973. Fluorescence and the location of tryptophan residues in protein molecules. Photochem. and Photobiol. 18: 263279.[Medline]
Bychkova, V.E., Berni, R., Rossi, G.L., Kutyshenko, V.P., and Ptitsyn, O.B. 1992. Retinol-binding protein is in the molten globule state at low pH. Biochemistry 31: 75667571.[CrossRef][Medline]
CCP4. 1994. The CCP4 suite: programs for protein crystallography. Acta Crystallogr. D 50: 760763.[CrossRef][Medline]
Clark, P.L., Liu, Z., Zhang, J., and Gierasch, L.M. 1996. Intrinsic tryptophans of CRABPI as probes of structure and folding. Protein Sci. 5: 11081117.[Abstract]
Cogan, U., Kopelman, M., Mokady, S., and Shinitsky, M. 1976. Binding affinities of retinol and related proteins to retinol binding proteins. Eur. J. Biochem. 65: 7178.[Medline]
Colantuoni, V., Romano, V., Bensi, G., Santoro, C., Costanzo, F., Raugei, G., and Cortese, R. 1983. Cloning and sequencing of a full-length cDNA for human retinol-binding protein. Nucl. Acids Res. 11: 77697777.