|Home | About | Journals | Submit | Contact Us | Français|
Bilirubin oxidase (BOD), a multicopper oxidase found in Myrothecium verrucaria, catalyzes the oxidation of bilirubin to biliverdin. Oxygen is the electron acceptor and is reduced to water. BOD is used for diagnostic analysis of bilirubin in serum and has attracted considerable attention as an enzymatic catalyst for the cathode of biofuel cells that work under neutral conditions. Here, the crystal structure of BOD is reported for the first time. Blue bipyramid-shaped crystals of BOD obtained in 2-methyl-2,4-pentanediol (MPD) and ammonium sulfate solution were merohedrally twinned in space group P63. Structure determination was achieved by the single anomalous diffraction (SAD) method using the anomalous diffraction of Cu atoms and synchrotron radiation and twin refinement was performed in the resolution range 33–2.3 Å. The overall organization of BOD is almost the same as that of other multicopper oxidases: the protein is folded into three domains and a total of four copper-binding sites are found in domains 1 and 3. Although the four copper-binding sites were almost identical to those of other multicopper oxidases, the hydrophilic Asn residue (at the same position as a hydrophobic residue such as Leu in other multicopper oxidases) very close to the type I copper might contribute to the characteristically high redox potential of BOD.
Bilirubin oxidase (BOD; EC 126.96.36.199), a member of the multicopper oxidase family, is found in Myrothecium verrucaria. It catalyzes the oxidation of bilirubin to biliverdin and then to a purple pigment in vitro (Shimizu, Kwon et al., 1999 ; Shimizu, Sasaki et al., 1999 ). In medical examinations of the liver, BOD is used for diagnostic measurement of the bilirubin concentration in serum (Kosaka et al., 1987 ; Shimizu, Kwon et al., 1999 ; Shimizu, Sasaki et al., 1999 ). Laccase, the most studied enzyme in this family, has great industrial potential in processes such as the decolorization of dyes, biopulping, organic synthesis and hair dyeing (Couto & Herrera, 2006 ).
Multicopper oxidases contain one type I, one type II and two (one pair of) type III coppers. The type II and III coppers form a trinuclear centre that reduces oxygen to water molecules and the type I copper functions as the electron mediator from the substrate to the trinuclear centre (Shleev et al., 2005 ; Enguita et al., 2004 ). Recently, BOD has been considered to be one of the best enzymes for reactions in biofuel cell cathodes because of its high reactivity at neutral pH (Cracknell et al., 2008 ; Tsujimura et al., 2001 ).
The structures of several multicopper oxidases, such as the CueO laccase from Escherichia coli (27% sequence identity; Roberts et al., 2003 ), the CotA laccase from Bacillus subtilis (35% sequence identity; Enguita et al., 2004 ), laccase from Trametes versicolor (Piontek et al., 2002 ) and ascorbate oxidase (Messerschmidt et al., 1992 ) have been reported. BOD is a monomeric protein (like other multicopper oxidases) with a molecular mass of 60 kDa. The structure of BOD has been predicted to be similar to those of other multicopper oxidases, but no crystal structure of BOD has been obtained.
In the present paper, we report the crystal structure of BOD. Blue bipyramid-shaped crystals of BOD were grown in MPD and ammonium sulfate solution using the hanging-drop vapour-diffusion method and were found to be merohedrally twinned in space group P63. The crystal structure of BOD was solved at 2.3 Å resolution by the SAD method using the anomalous diffraction of Cu atoms. The structure of BOD is similar to those of other multicopper oxidases such as the laccases CueO and CotA. The protein is folded into three domains and four copper-binding sites are located in domains 1 and 3. The copper-binding sites are almost identical to those in other multicopper oxidases.
BOD from M. verrucaria (MT-1) was purified as described previously (Shimizu, Kwon et al., 1999 ; Shimizu, Sasaki et al., 1999 ). The sequence-database reference code for the protein is UniProtKB/Swiss-Prot Q12737. The sample was dialyzed against 20 mM Tris–HCl buffer pH 8.0 and further purified on a Mono Q column (Mono Q 5/50 GL, GE Healthcare, Buckinghamshire, England) equilibrated with 20 mM Tris–HCl pH 8.0 and eluted with a linear gradient of buffer containing 0–1 M KCl at a flow rate of 1.5 ml min−1. The purified sample was concentrated to about 10 mg ml−1 using Centricon-10 or Ultrafree-4 centrifugal filter units (Millipore, Bedford, Massachusetts, USA). The protein concentration was determined using a molar absorption coefficient at 600 nm of 600 = 5000 M −1 cm−1. Other chemicals were of guaranteed grade from Wako Pure Chemical Industries (Osaka, Japan), Nacalai Tesque (Kyoto, Japan) and Sigma–Aldrich Chemicals (St Louis, Missouri, USA).
The concentrated purified BOD sample (9.2 mg ml−1 in 10 mM Tris–HCl pH 8.0) was crystallized using the hanging-drop vapour-diffusion method. The solution in the crystallization drop was prepared on a siliconized cover glass by mixing identical volumes (3 µl + 3 µl) of protein solution and precipitant solution [10% 2-methyl-2,4-pentanediol (MPD), 1.44 M ammonium sulfate, 10% glycerol and 0.5 M KCl]. The droplets were equilibrated against 0.2–1.0 ml precipitant solution at 293 K. The obtained crystals were light blue in colour and bipyramidal in shape (spindle-shaped); the radius of the crystals was around 0.4 mm.
The crystal was mounted in a nylon CryoLoop (Hampton Research, Aliso Viejo, California, USA) and placed into liquid nitrogen. Cryoprotectant (10% MPD and 10% glycerol) was included in the precipitant solution used for crystallization. The X-ray fluorescence spectrum of the sample was measured for Cu atoms and the peak wavelength was determined to be 1.37751 Å. Diffraction data were collected in a cold nitrogen-gas stream (100 K) on an ADSC Quantum 4 CCD area detector (Rigaku/MSC, The Woodlands, Texas, USA) using synchrotron radiation at a wavelength of 1.37751 Å at the BL38B1 station of SPring-8 (Hyogo, Japan). The resulting data set was processed, merged and scaled using the HKL-2000 program package (HKL Research, Charlottesville, Virginia, USA) to a resolution of 2.5 Å (Table 1 ). Diffraction data from another crystal were collected at a wavelength of 1.0 Å and processed to a resolution of 2.3 Å in the same way (Table 2 ).
The initial phase determination as a single crystal (not as a twinned crystal) was performed using the SAD method implemented in the program SOLVE/RESOLVE (Terwilliger, 2003 , 2004 ; Table 1 ). Using the eight Cu atoms found in the structure, a model of B. subtilis CotA laccase (PDB code 1gsk or 1uvw; Enguita et al., 2003 , 2004 ) was roughly superimposed on the calculated electron-density map. The amino-acid sequence of CotA was changed to that of BOD and initial modelling was performed using the TURBO-FRODO program (AFMB, Marseille, France). Rigid-body refinement and restrained least-squares refinement were performed using the REFMAC program (Collaborative Computational Project, Number 4, 1994 ). Several rounds of refinement as a twinned crystal (with twin law h, −h − k, −l determined using phenix.xtriage), followed by manual model building using Coot (Emsley & Cowtan, 2004 ), were carried out using SHELXL (Sheldrick, 2008 ) and phenix.refine (Adams et al., 2002 ). Final twin refinement was performed using phenix.refine. The atomic coordinates of the BOD structure have been deposited in the Protein Data Bank (code 3abg).
Structure-similarity searches using BOD as a query were performed using the program DaliLite (Holm et al., 2008 ; Holm & Park, 2000 ). Two multicopper oxidases, the laccases CueO and CotA, which were determined by DaliLite to have the same highest Z scores (46.0), were compared with BOD using the SSM algorithm implemented in Coot and the LSQKAB least-squares fitting in the CCP4 package.
The molecular dynamics of a BOD–mediator complex were calculated using GROMACS (Van Der Spoel et al., 2005 ) and the GMX force field. The coordinates of the mediator compound 2,2′-azino-bis(3-ethylbenzthiazoline-6-sulfonic acid) (ABTS) were taken from the structure of the CotA–ABTS complex (PDB code 1uvw). A model of the BOD–ABTS complex was made by superimposing (using Coot) the main chains of BOD and CotA. Water molecules were added to the model and molecular-dynamics calculations were performed for 10.0 ps (5000 steps of 0.002 ps).
Although a clear diffraction pattern was observed for the BOD crystals, the graphs of the cumulative intensity distribution and the fourth moment of E calculated using the TRUNCATE program indicated the possibility of perfect twinning (data not shown). The program phenix.xtriage also suggested a high possibility of twinning with twin law (h, −h − k, −l) and I 2/I2 = 1.615. Therefore, the data sets were processed as P63 (not P6322 as misleadingly presumed by visual examination of precession images from the data sets). Phase determination and initial refinement calculations could be carried out using P63 twinned data applied to a single-crystal (untwinned) model and statistics for the phase determination are shown in Table 1 . Final twin refinement was calculated with the twin law (h, −h − k, −l) and the twin fraction was determined to be 0.487 (almost equal to 0.5). Statistics for the refinement are summarized in Table 2 .
The complete BOD molecule comprises 572 amino-acid residues. However, residues 509–572 of molecule A and 519–572 of molecule B were not included in the final model because no clear interpretable electron density was observed for these residues. The asymmetric unit contained two molecules and the root-mean-square deviation (r.m.s.d.) between molecules A and B was 0.55 Å for all Cα atoms. Weak (the coordinate_sigma parameter was 0.5 in phenix.refine) noncrystallographic symmetry restraints for residues 1–508 were applied during refinement. Four carbohydrate chains in molecules A and B were added to residues Asn472 and Asn482. Two GlcNAcs and one mannose were added to Asn472 and two GlcNAcs were added to Asn482. Relevant refinement statistics are given in Table 2 . The overall completeness, R factor and R free values were 96.7%, 0.205 and 0.263, respectively, for the data used in the refinement. For the highest resolution bin (2.34–2.30 Å), the completeness was 98.0% and the R factor and R free values were 0.299 and 0.287, respectively. Because twinned refinement can lead to artificially low R values, parts of the reported structure may not be free of significant errors. In the Ramachandran plot of the main-chain torsion angles, 61.9% of the residues are in the most favoured regions, with 97.5% of the residues lying within the generously allowed regions as defined in the program PROCHECK (Collaborative Computational Project, Number 4, 1994 ). The low percentage in the most favoured region is probably a consequence of the difficulty of refinement calculations using a twinned crystal.
Fig. 1 displays the overall structure of BOD (chain A in the asymmetric unit) as a ribbon representation. As in other multicopper oxidases, BOD is folded into three domains and the four metal-binding sites are located in the clefts formed by domains 1 and 3 (Fig. 1 a). The two multicopper oxidases CueO laccase and CotA laccase had the highest Z scores (46.0) from the DaliLite structure-similarity search server. The r.m.s.d.s calculated using the SSM superposition in Coot were 1.7 Å between BOD and CueO for the Cα atoms of 403 corresponding residues and 1.5 Å between BOD and CotA for 428 residues (Fig. 1 b).
The type I Cu is coordinated by four amino-acid residues (His398, Cys457, His462 and Met 467) and the type II Cu is coordinated by two amino-acid residues (His94 and His401) (listed in Table 3 ). The type III copper coordination by six amino-acid residues (His136, His403 and His456; His134, His96 and His458) and a water molecule is shown in Fig. 2 . Although a liganded water molecule was included between the two type III coppers in the model of molecule B, the B factor (40.2 Å2) of the O atom of the liganded water was considerably higher than that (26.6 Å2) in molecule A. The distances from the coppers to the liganded amino-acid residues in the molecule are summarized in Table 3 . As shown in Fig. 2 , clear electron density was observed for all four coppers in the 2F o − F c electron-density map and the anomalous map calculated using the diffraction data used for the SAD method.
A molecular-dynamics simulation was performed for the BOD–mediator (ABTS) complex (the initial model of BOD–ABTS was produced as described in §2.6). The average structure after the simulation shows the aromatic rings of the ABTS molecule and the Trp396 residue of BOD, indicating that the binding of the mediator molecule to BOD is stabilized by π–π interactions.
In the study described here, the copper-bound structure of M. verrucaria BOD was determined at 2.3 Å resolution using merohedrally twinned crystals as the first demonstration of the structure of BOD. The type I Cu site of BOD has been reported to have a high redox potential compared with other multicopper oxidases with Met as an axial ligand (Tsujimura et al., 2005 ; Christenson et al., 2006 ). This study provides important insights regarding the characteristics of the redox potential of BOD.
The overall structure of BOD was similar to those of two multicopper oxidases (CueO laccase from E. coli and CotA laccase from B. subtilis) deposited in the PDB, as shown in Fig. 1 (b). The coordination of the type I Cu by four amino-acid residues (His398, Cys457, His462 and Met 467) and the coordination of the type II Cu by two amino-acid residues (His94 and His401) were almost the same as those in CueO (Fig. 3 ; Roberts et al., 2003 ). The coordination of the type III Cu by six amino-acid residues (His136, His403 and His456; His134, His96 and His458) and a water molecule was also similar to that in CueO. The distances between the Cu atoms and the liganded amino-acid residues in BOD, CueO and CotA are compared in Table 3 . All bond lengths between coppers and liganded amino-acid residues in each molecule (A and B) are similar to those in CueO and CotA. While the coordination patterns of the principal residues are similar, BOD has a polar residue Asn459 at the opposite position to the axial ligand (Met467) of the type I Cu; the distance between Asn459 Nδ2 and the copper ion is 4.61 Å. CueO and CotA have nonpolar residues (Leu and Ile, respectively) at the position corresponding to Asn459 of BOD. The unique environment around the type I Cu might be responsible in part for the high redox potential. The redox potential of BOD is higher than those of CueO and CotA by 200–230 mV.
A molecular-dynamics simulation of the BOD–mediator (ABTS) complex suggested that the binding of the mediator molecule to Trp396 of BOD is stabilized by a π–π interaction. The corresponding Trp residue is not found in CotA or CueO. This coupling, which is not highly dependent on pH, might contribute to the higher activity at neutral pH of BOD.
The synchrotron-radiation experiments were performed at BL38B1 of SPring-8 with the approval of the Japan Synchrotron Radiation Research Institute (JASRI).