HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

J. Biol. Chem., Vol. 275, Issue 40, 31086-31092, October 6, 2000

This Article Abstract Full Text (PDF) All Versions of this Article: 275/40/31086    most recent M004281200v1 Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Request Permissions Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Yano, J. K. Articles by Poulos, T. L. Search for Related Content PubMed PubMed Citation Articles by Yano, J. K. Articles by Poulos, T. L. Social Bookmarking                      What's this?

Crystal Structure of a Thermophilic Cytochrome P450 from the Archaeon Sulfolobus solfataricus^*

Jason K. Yano, Laura S. Koo^§, David J. Schuller, Huiying Li, Paul R. Ortiz de Montellano^§, and Thomas L. Poulos^¶

From the  Department of Molecular Biology and Biochemistry and Program in Macromolecular Structure, University of California, Irvine, California 92697-3900 and the ^§ Department of Pharmaceutical Chemistry, School of Pharmacy, University of California, San Francisco, California 94143-0446

Received for publication, May 18, 2000, and in revised form, June 15, 2000

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

The structure of the first P450 identified in Archaea, CYP119 from Sulfolobus solfataricus, has been solved in two different crystal forms that differ by the ligand (imidazole or 4-phenylimidazole) coordinated to the heme iron. A comparison of the two structures reveals an unprecedented rearrangement of the active site to adapt to the different size and shape of ligands bound to the heme iron. These changes involve unraveling of the F helix C-terminal segment to extend a loop structure connecting the F and G helices, allowing the longer loop to dip down into the active site and interact with the smaller imidazole ligand. A comparison of CYP119 with P450cam and P450eryF indicates an extensive clustering of aromatic residues may provide the structural basis for the enhanced thermal stability of CYP119. An additional feature of the 4-phenylimidazole-bound structure is a zinc ion tetrahedrally bound by symmetry-related His and Glu residues.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

Understanding the structural basis for the enhanced stability of proteins from thermophilic organisms relative to their mesophilic counterparts is a challenging problem. The increasing availability of high-resolution crystal structures from both thermophiles and mesophiles has provided an important database for determining the structural basis of thermal stability. For example, Gromiha et al. (1) analyzed the thermal stability of 16 different families of proteins. The collective results suggest that the stability of thermally stable proteins may result from a better balance between packing and solubility. Other researchers have attributed various factors such as changes in loop flexibility (2), salt bridge interactions (3), and amino acid substitutions (i.e. Ile for Ala) (1) as the origin of increased thermal stability. It is important to note that despite all efforts to define specific traits, no single trait has been universally linked to the increase in thermal stability for all proteins. It is therefore important to use comparative analyses to determine the factors that lead to an increase in the thermal stability of structurally related proteins.

In addition to the fundamental question of how structure relates to stability, there is considerable practical significance in the industrial and medical application of thermally stable enzymes. (4). In this regard, reactions catalyzed by cytochromes P450 are of particular interest. P450s are heme enzymes that catalyze the following reaction, where R represents a variety of substrates, normally non-polar aromatic or aliphatic molecules.

One of the primary functions of P450s is to aid in the clearance of toxic hydrocarbons by rendering such molecules more soluble through hydroxylation. In addition to this catabolic function, P450s also participate in the biosynthesis of essential molecules such as steroids (5) and various natural products including antibiotics (6).

A long sought practical goal in P450 research is to capitalize on the exquisite specificity of P450s in regio- and stereo-selective hydroxylation reactions (7). One possible example is in the preparation of important hydroxylated molecules that are difficult to prepare using traditional organic synthesis methods. In addition, P450s could prove useful in the oxidative or reductive elimination of environmental pollutants (8). A current limitation is the requirement for relatively expensive reducing equivalents in the form of NADH or NADPH. It might be possible to circumvent this restriction by the use of H₂O₂ as the oxidant rather than O₂, because H₂O₂-supported P450 reactions do not require the transfer of electrons from a redox partner (9). A second limitation is the instability of the P450 catalysts. There are basically two approaches to the resolution of this problem. One is to seek out enzymes from thermophiles that catalyze the reaction of interest, and the second is to use protein engineering to convert a mesophilic enzyme into a thermally stable enzyme. Utilization of the latter approach requires knowledge of those structural features that are critical for imparting thermal stability. Until now this has not been possible with P450s because the available structures are all from mesophiles. However, quite recently, Wright et al. (10) accidentally discovered a P450 while attempting to clone the thymidylate synthase gene from the acidothermophilic archaeon Sulfolobus solfataricus. Based on the accepted P450 nomenclature (5), this P450 has been termed CYP119.¹ S. solfactaricus has optimal growth conditions at 85 °C and pH 3.5 (11), which operationally defines it as a thermophile. As expected, CYP119 exhibits unusual thermal stability with a melting temperature of ~90 °C compared with ~55 °C for bacterial P450cam (12, 13). The function of CYP119 remains unknown, but considering that S. solfataricus can utilize elemental sulfur as an energy source (11) and that P450s are known to participate in sulfoxide formation (14), the oxidative metabolism of sulfur compounds is a possibility. Currently, the genome of S. solfataricus is being sequenced (15) and once completed may provide some hints on the function and redox partners of CYP119. Here we describe the structure determination of CYP119 and a comparison of the structure with those of its mesophilic counterparts in an attempt to understand the structural basis for thermal stability.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

Expression and Purification of CYP119-- BL21(DE3) Escherichia coli cells were transformed with a pCWori plasmid containing the gene for CYP119 (12). Single colonies were selected for overnight growth in 5-ml cultures, which were used to inoculate 1-liter cultures. Two overnight 1-liter cultures were used to inoculate a 40-liter fermenter. At an A_{600 nm} of 0.8, 500 µg/liter of isopropyl -D-thiopyranogalactoside (USB) was used to induce expression of CYP119. The conditions for cell growth were 2× YT (tryptone, 16 g, yeast extract, 10 g, NaCl, 5 g, per liter), 0.1 g/liter ampicillin (Sigma) at 37 °C with shaking at 220 rpm for the starter cultures and slow stirring (80-100 rpm) at 30° C for the fermenter. Thirty-six hours after induction the cells were harvested, and cell paste was stored at 70 °C.

Cells were lysed in a French pressure cell in 50 mM Tris-Cl, pH 7.4, 0.1% Triton X-100, and 1 mM phenylmethylsulfonyl fluoride (Sigma). Cell debris was pelleted by centrifugation, and the cell-free extract was heated to 65 °C for 1 h. Precipitated proteins were removed by centrifugation. The supernatant was brought to 80% ammonium sulfate, and the resulting precipitant was resuspended in minimal volume and dialyzed against 50 mM bis-Tris, pH 7.2 overnight (three times buffer exchanges). The dialyzed protein was loaded onto a QAE-Sepharose column (2.5 cm × 25 cm), washed with 2-3 column volumes of buffer, then eluted with a linear gradient, 0-250 mM NaCl in 50 mM bis-Tris, pH 7.2. Fractions with an absorbance ratio (A_{415 nm}/A_{280 nm}) of > 1.40 were pooled, desalted, and exchanged into 50 mM bis-Tris, pH 7.0. A chromatofocusing column (PBE94, Amersham Pharmacia Biotech) was equilibrated with 50 mM bis-Tris, pH 6.0. The protein was loaded, washed, and then eluted with Polybuffer 74 pH 5.0 (Amersham Pharmacia Biotech) according to the manufacturer's instructions. Fractions with an absorbance ratio (A_{415 nm}/A_{280 nm}) above 1.60 were pooled, desalted, and then exchanged to 50 mM bis-Tris, pH 6.0 and used for crystallization. P450 concentrations were determined with the extinction coefficient _{415 nm} = 104 mM¹ as reported by McLean et al. (13).

Crystallization-- Crystals belonging to the monoclinic space group P2₁ (unit cell dimensions a = 76.97 Å, b = 70.26 Å, c = 107.54 Å,  = 102.34°) were grown by free interface diffusion in 200-µl capillaries with a protein to precipitant ratio of 2:3. The initial protein concentration for capillary growth was 40 mg/ml. The precipitant solution consisted of 35% polyethylene glycol (PEG) 3350, 100 mM imidazole pH 8.0, and 0.2 M Li₂SO₄. Crystals belonging to the orthorhombic space group C222₁ (unit cell dimensions a = 76.96 Å, b = 114.60 Å, c = 185.20 Å) were grown by hanging-drop vapor diffusion over a well of 20% PEG 3350, 100 mM EPPS, pH 8.5, 0.2 M MgSO₄. 1 mM 4-phenylimidazole was mixed with 20 mg/ml CYP119 (final concentrations) and used to set up 5 µl + 5 µl drops on siliconized glass coverslips. All crystals were grown at room temperature. The single crystals used for data collections were carefully separated from larger crystal clusters. Cryogenic data collection conditions for the C222₁ space group consisted of a four-step transfer to artificial precipitant solution with increased ethylene glycol concentration up to 20%. Cryogenic conditions for the P2₁ space group were unsuccessful, leading to severe anisotropy in crystal diffraction. There were two molecules per asymmetric unit for both space groups.

Data Collection-- Room temperature data were collected from the P2₁ crystal form using an in-house R-Axis IV imaging plate equipped with a Rigaku rotating copper anode x-ray generator. High-resolution data collection (50-1.93 Å) at 100 K on the C222₁ crystal form was performed at the Stanford Synchrotron Radiation Laboratory, beamline 7-1 with a Mar345 imaging plate. Low-resolution data (50-3.1 Å) on the C222₁ form were also collected in-house at 100 K. Optimization of data collection was guided by the STRATEGY function of MOSFLM (16). All data were reduced using DENZO and SCALEPACK (17), and rejections were performed with ENDHKL (Louis Sanchez, California Institute of Technology) in conjunction with SCALEPACK.

Molecular Replacement-- Molecular replacement in the P2₁ space group was carried out with a version of BRUTE modified to use log-likelihood scoring² using data to 4.5 Å with a polyalanine model of P450eryF (PDB accession number 1OXA) (18) as the search model. The top rotation solution was correct and was further optimized using data to 3 Å, which was followed by a translational search. The second rotation solution was generated through application of the self-rotation peak and followed with another translational search. The final solution had an R-factor of 0.53 after rigid body refinement. Density modification in MAGICSQUASH (19) was used to improve phases. Electron density map fitting beginning from the polyalanine P450eryF model was carried out with TOM (20). Several sections of the electron density map, especially the core of the structure near the heme, clearly showed the identity of CYP119 side chains. Regions of the molecule that lacked clear backbone electron density were deleted prior to further refinement with CNS (21). With successive rounds of refinement and model building the entire CYP119 sequence was fit to the electron density map except for the last two residues. In the first 12 rounds of refinement, the entire chain of one molecule was traced and non-crystallographic symmetry operators were used to generate the second molecule in the asymmetric unit. In the last two rounds of refinement, each chain was fit separately. The final R-factor was 0.219 (R-free = 0.253). Backbone geometry was analyzed in PROCHECK (22) and only two residues (alanine 152 in both molecules) were in the disallowed region. There were 29 water and 2 sulfate molecules found per asymmetric unit.

To determine phases for the C222₁ space group the imidazole-bound (P2₁) structure was used as a search model in AMoRe (23). The first molecule produced a solution that stood out above the background and was fixed to find the second molecule giving a R-factor of 0.436. After one round of both positional and temperature factor refinement in CNS (21) the R-factor dropped to 0.314 (R-free = 0.348) at a resolution of 1.93 Å. The F and G helical regions of the protein had to be extensively remodeled. The final R-factor was 0.225 (R-free = 0.274). Backbone geometry was checked in PROCHECK (22), and none of the residues were in the disallowed region. All of the residues except the last residue, 368, were present in the electron density maps. The final model contains 622 water molecules, 2 sulfates, and one zinc per asymmetric unit. Relevant statistics are given in Table 1. 

The coordinates have been deposited to the Protein Data Bank and the PDB codes are 1F4U (imidazole-bound) and 1F4T (4- phenylimidazole-bound).

	RESULTS AND DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

Overall Structure-- CYP119 exhibits the typical P450-fold (Fig. 1). However, at 368 residues, CYP119 is considerably shorter than P450cam (24) or P450eryF (18), at 414 and 403 residues, respectively. The majority of this difference in length is located at the N-terminal region. P450cam and P450eryF begin with relatively ill-defined N-terminal tails whereas CYP119 begins immediately with helix A. As a result, residue 1 in CYP119 corresponds to residue 42 or 16 in P450cam or P450eryF, respectively. Therefore, the shorter N terminus alone in CYP119 accounts for 41 of the 56 fewer amino acids in CYP119 compared with P450cam. The remainder of the differences in length occurs primarily in surface turns (Fig. 1). For example, the 1-1/1-2 hairpin (residues 13-24) in CYP119 is four residues shorter than that in P450cam (residues 52-66). Another difference includes the 225-234 5-turn in P450cam compared with the 191-195 loop in CYP119, a difference of five residues. This tight loop between the H and I helices in CYP119 resembles what was found in P450nor (25), which does not require a redox protein partner for electron transfer. In contrast, all other known P450 structures show an extended but ill-defined -hairpin similar to that in P450cam. Interestingly, this hairpin region makes direct contact with the FMN-binding domain in the crystal structure of the complex formed between the heme and FMN domains in P450 BM-3 (26).

View larger version (79K): [in this window] [in a new window]   Fig. 1.   Ribbon diagram of CYP119 with 4-phenylimidazole ligated to the heme group. Helices are represented by white cylinders whereas -structures are represented by thick gray arrows. All figures were prepared with the program MOLSCRIPT (37) or SETOR (38) and rendered in Raster3D (39).

The most striking differences involve the B' helix region, which is known to be important in substrate binding (27). Residues 63-66 in CYP119 corresponding to the location of B' helix in P450eryF are no longer a complete helix. Instead this region defines only one turn of the 3₁₀ helix. However, CYP119 does have a helix that we have termed the B' helix consisting of residues 49-53 (Fig. 2). In CYP119, the B' helix is preceded by a long loop that enables the B' helix to adopt a position away from the active site toward the molecular surface (Fig. 2). Such a change should leave the active site relatively exposed. However, a solvent accessibility calculation, without substrate or inhibitors included, shows that ~24 Ų of the heme is exposed in CYP119 compared with ~18 Ų in P450cam. The reason that CYP119 and P450cam are so similar in active site access by a solvent probe is because of the repositioning of the loop connecting the F and G helices (Fig. 2). In P450cam and other P450s, this F/G loop extends out on the surface. Conversely, in CYP119, this loop points completely in the opposite direction and dips into the active site, thereby occupying the space normally taken by the B' helix region. As we shall see in the next section, this loop can adopt alternate conformations depending on the identity of the ligand coordinated to the heme iron.

View larger version (21K): [in this window] [in a new window]   Fig. 2.   A stereo diagram showing the comparison of the F/G loop region in the CYP119 imidazole complex (shaded) with P450cam (clear). For this diagram, the molecules were superimposed using only the heme groups. Note that the F/G loop in P450cam points up and away from the heme whereas in CYP119 the F/G loop points down toward the distal (substrate-binding) side of the heme. Note, too, that the F/G loop in CYP119 is positioned where the B' helix is located in P450cam.

The most highly conserved regions are near the heme, which includes the I and L helices. The polypeptide conformation and local environment surrounding the Cys thiolate heme axial ligand, Cys-317 in CYP119, is strictly conserved. The I helix spans the entire molecule and is situated directly above the heme. The I helix in CYP119 is unusual, having two additional Thr residues following the conserved Thr at position 213 (Fig. 3), which is present in most P450s. The positioning of Thr-213 close to the heme surface relative to Thr-214 and Thr-215 was correctly predicted based on chemical modification data (12). This region has been postulated to be involved in a proton shuttle network considered important for delivering solvent protons required for the activation of O₂ during the catalytic cycle (28). The H-bonding network involving helix I Thr residues, however, is somewhat different in CYP119. The side chain hydroxyl group of the conserved Thr-213 does not donate an H-bond to a peptide carbonyl as it does in P450cam. Instead, Thr-214 donates an H-bond to the peptide oxygen atom of Gly-210 (Fig. 3). This is very similar to the H-bonding pattern found in P450BM-3 that also has a Thr corresponding to Thr-214 in CYP119 (29). Interestingly, the T214A (but not T213A) mutant results in a large increase in the rate at which CYP119 catalyzes H₂O₂-supported styrene epoxidation (12). Because Thr-214 but not Thr-213 H-bonds with the backbone, it might be expected that mutation of Thr-214 would cause a greater perturbation in the I helix thereby leading to a greater change in reactivity toward styrene.

View larger version (23K): [in this window] [in a new window]   Fig. 3.   A segment of helix I showing the location of the conserved Thr-213 and the additional two Thr residues in CYP119. Note that the side chain oxygen atom of Thr-213 donates an H-bond to the carbonyl oxygen of Gly-210.

Zinc Binding Site-- An interesting feature of the 4-phenylimidazole structure is a cation, most likely Zn²⁺, situated at the interface between the two molecules in the asymmetric unit (Figs. 4 and 5). The identification of the ion as Zn²⁺ is based on anomalous difference Fouriers calculated in PHASES (30) at different x-ray wavelengths. Zinc has an absorption edge at 1.28 Å whereas that for iron is at 1.74 Å. Using data from the copper rotating anode at a wavelength of 1.54 Å, an anomalous difference Fourier shows strong 10-11 peaks for the two iron atoms but nothing for the presumed zinc. However, using synchrotron data at a wavelength of 1.08 Å, the zinc site has a peak of 17, whereas the iron peaks are 13-14. The type of ligands and the tetrahedral coordination environment of the ion also support its identification as zinc. As shown in Fig. 5, the Zn²⁺ is coordinated by Glu-139 and His-178 in one subunit, and their symmetry mates in the other subunit. The zinc binding motif is very likely to be native to the enzyme rather than an artifact of crystal lattice formation because Zn²⁺ was not included in the crystallization buffers. The fact that Zn²⁺ was not seen in the imidazole complex may be due to the presence of 0.1 M imidazole in the buffer, which may have scavenged the low amounts of Zn²⁺ present as a contaminant. It remains unclear if the dimer and Zn²⁺ site are functionally important. Only 727 Å² per monomer of accessible surface area is buried at the interface, which is below the range for most true dimers. However, the tight coordination of Zn²⁺ might provide sufficient stabilization to favor the dimeric form under physiological conditions.

View larger version (37K): [in this window] [in a new window]   Fig. 4.   A ribbon diagram showing the location of the Zn2+ ion at the dimer interface in 4-phenylimidazole complex.

View larger version (45K): [in this window] [in a new window]   Fig. 5.   The dimeric interface in CYP119 and a close-up view of the tetrahedral ligation of the putative zinc ion in the 4-phenylmidazole-bound structure. The zinc ion is ligated by Glu-139 and His-178 and their symmetry mates forming a tetrahedral geometry common in known crystal structures with bound zinc (40).

Comparison of the Two Crystal Forms-- A comparison of the imidazole and 4-phenylimidazole structures reveals interesting and unexpected differences. The active site adjusts in size and shape depending on the ligand bound, with the active site expanding to accommodate the larger 4-phenylimidazole ligand. The F/G region undergoes the largest movements with backbone atoms shifting as much as 6 Å. The most significant change occurs in the F helix. In the 4-phenylimidazole complex, the F helix includes residues 141-155, but in the imidazole complex the helix stops at residue 151, a loss of one full turn of helix (Fig. 6). This unraveling of the helix lengthens the F/G loop, which then is able to dip down into the active site to make contacts with the smaller imidazole ligand. For example, the Leu-155 side chain moves approximately 6 Å, which enables the side chain to interact with the imidazole ligand. The loop "in" conformation is stabilized by direct contacts between the F/G loop with helix I, which are not present in the 4-phenylimidazole complex (Fig. 7).

View larger version (29K): [in this window] [in a new window]   Fig. 6.   Ribbon diagram of the F-G region of CYP119. The F-G region encompasses residues 134-183. The white strand represents the F-G region of the imidazole-bound CYP119 structure. The gray strand represents the F-G region of the 4-phenylimidazole-bound structure. This representation shows the unwrapping of the F helix to accommodate the smaller imidazole ligand.

View larger version (41K): [in this window] [in a new window]   Fig. 7.   Stereo views of the electron density surrounding a portion of the F-G region. The model appears in light gray and black with the 2Fo-Fc electron-density map in thin dark gray lines. Residues Arg-154, Leu-155, and Glu-212 are shown for both models. a, the imidazole-bound structure (loop in). Arg-154 forms a hydrogen bond with Glu-212 causing Leu-155 to protrude into the active site. b, the 4-phenylimidazole-bound structure (loop out). Arg-154 no longer forms a hydrogen bond with Glu-212 and Leu-155 moves out of the way to accommodate the larger ligand. The omit electron density maps were generated with the atoms shown removed from a round of simulated-annealing refinement.

In the 4-phenylimidazole complex, the F/G loop must move to make room for the larger ligand. One of the most dramatic changes is in the location of Arg-154. In the 4-phenylimidazole complex, Arg-154 is part of the F helix where it H-bonds with the peptide oxygen atom of Glu-198 (Fig. 7). In the imidazole complex, this section of the F helix unfolds enabling the Arg-154 side chain to reorient by 180^o where it can then interact with Glu-212. This new ion pair in the imidazole complex and other F/G loop-I helix contacts helps to stabilize the loop "in" conformation. It appears that favorable interactions lost upon unfolding of the F helix are partially offset by favorable interactions between the F/G loop and helix I.

Another region that experiences a large change is the 4 hairpin turn centered on Val-353, which moves about 3 Å. In this case, the hairpin squeezes in closer to the ligand in the 4-phenylimidazole complex which is in the F/G loop "out" conformation. This location of the turn is incompatible with the F/G loop "in" conformation in the imidazole complex. The F/G loop and the 4 hairpin as a whole undergoes an "induced fit" conformational change depending upon the size and shape of the ligand bound at the heme active site.

The only other P450 in which similar changes have been observed is P450BM-3. In that case, there is a large difference in the position of the F and G helices as well as the F/G loop (29). In going from the substrate-free to -bound complexes, the entry to the active site closes down around the substrate. However, secondary structural elements in P450BM-3 move as a unit with no net gain or loss in secondary structure as observed with CYP119. The various studies with P450cam provide better analogies with CYP119 because with P450cam there are several structures available of the enzyme complexed with substrate analogues as well as various imidazole complexes. In these cases, the changes are very modest compared with what we observe in CYP119 and involve primarily rearrangements in water structure, side chains, and small movements of backbone atoms. The one very dramatic exception is a chemically modified derivative of P450cam where the Cys residues were derivatized with N-(2-ferrocenylethyl)maleimide. In this case, part of the active site essentially unfolds to enable the ferrocene moiety attached to Cys-85 to enter the active site resulting in changes of as much as 10 Å in the position of some residues (31). Whereas demonstrating an unexpected level of flexibility in the active site, the introduction of a covalently attached ferrocene is a very large change compared with the difference in going from imidazole to 4-phenylimidazole in CYP119. This indicates that the CYP119 active site pocket is unusually sensitive to subtle changes in the types of molecules bound at the active site compared with P450cam. Whether or not these changes are functionally relevant must await the discovery of the true function of CYP119.

Structural Basis for Thermal Stability-- The enhanced thermal stability of proteins from extreme thermophiles has been attributed to a number of factors. For example, an increase in Arg content can aid in thermal stability because of the additional H-bonding possibilities with Arg compared with Lys as well as the partial aromatic character of the guanidinium group (32). However, CYP119 contains a normal Arg content, 7.6% of the total amino acids, compared with P450cam and P450eryF, 8.4 and 6.2%, respectively. CYP119 also does not contain an unusual proportion of aliphatic and aromatic residues. In addition, the fraction of total accessible surface area due to nonpolar (aromatic and aliphatic) residues in CYP119 is relatively high at 17%, compared with 11 and 12% for P450cam and P450eryF, respectively.

Where CYP119 does differ is in the number and type of salt linked networks and aromatic interactions. Salt bridges were computed with the program HBPLUS (33) and the results are summarized in Table 2. CYP119 has a lower total number of 2-residue salt bridges but a larger number of salt-bridged networks. Whereas the two crystal forms of CYP119 differ slightly in salt bridges (Table 2), both CYP119 structures have four, P450cam has three, and P450eryf has only one salt-bridged network. The presence of only one additional salt-linked network in CYP119 compared with P450cam is unlikely to be the primary reason for enhanced stability. There is, however, an additional interesting difference in salt bridges. In P450cam the longest distance between interacting side chains in a salt-bridged network is 78 residues whereas the longest in CYP119 is 287 residues; Arg-9 in CYP119 interacts with Glu-296, which effectively ties the N- and C-terminal regions together in a salt bridge. We doubt, however, that a single long-range salt bridge network can account for a large portion of the enhanced stability of CYP119 compared with P450cam.

The most striking difference that might be associated with enhanced thermal stability is a unique clustering of aromatic residues in CYP119. A recent homology model of CYP119 (34) correctly predicted the aromatic clustering found in the crystal structure. Fig. 8 displays all the aromatic residues in CYP119, P450cam, and P450eryF. CYP119 has two sets of aromatic clusters not found in the other P450s. The first cluster involves five residues (Tyr-2, Trp-4, Phe-5, Phe-24, and Trp-281) that span a distance of ~11.3 Å between -carbons. Tyr-15 could also be considered part of this cluster because Met-8 contacts both Tyr-15 and Phe-24. The second cluster consists of seven residues: Phe-225, Phe-228, Trp-231, Tyr-250, Phe-298, Phe-334, and Phe-338 and spans a distance of ~24 Å between -carbons. Cluster 1 ends with Tyr-2 whereas cluster 2 begins with Tyr-250. The distance between the two coplanar aromatic rings is ~7 Å. However, the guanidinium group of Arg-287, which ion pairs with Asp-296, is stacked between the two Tyr side chains thus providing a continuous  stacking interaction between the two aromatic clusters. This means a continuous aromatic/nonpolar "ladder" involving all aromatic residues, and one arginine spans the entire edge of CYP119 covering a distance of ~39 Å. P450eryF, with a single cluster of six aromatic residues that spans a distance of ~11 Å comes closest to matching this arrangement. Based on the analysis thus far, this is the one obvious difference that is most likely to be a critical part of the enhanced thermal stability properties of CYP119. It is of interest to note that a recent mutagenesis study has shown that increasing aromatic interactions by introducing only one additional aromatic side chain increased thermal stability by 9 °C in a xylanase (35). In addition, a single point mutation (F31Y) in ribonuclease P2 from S. solfactaricus disrupted an aromatic cluster which " ... led to an unprecedented loss" of thermal- and baro-stability by at least 27 K and 10 kilobar, respectively (36). Future mutagenesis work on CYP119 should reveal if the aromatic clustering revealed by the CYP119 crystal structure is the key to its increased thermal stability.

View larger version (53K): [in this window] [in a new window]   Fig. 8.   Aromatic networks in P450s. CYP119 (a) and the enlargement of the aromatic network in CYP119 (b) are shown. Arg-287 marks the transition point between the two networks. P450cam (c) and P450eryF (d) are shown for comparison.

Conclusions-- As expected, the CYP119 structure closely resembles other known P450 structures. CYP119 is shorter owing primarily to fewer residues at the N terminus and two surface -hairpins. Unlike what has been observed in other P450s, the F/G loop connecting the F and G helices can adapt dramatically different conformations to optimize interactions with active site ligands of different sizes and shapes. Such changes may foreshadow the type of flexibility required for those P450s that are able to metabolize a wide range of substrates, especially the key drug metabolizing P450s. One important structural feature that may be related to thermal stability is a large cluster of aromatic residues in CYP119 that is not found in the other P450s.

	ACKNOWLEDGEMENTS

We thank Jean-Philippe Cartailler for help with all of the figures and Irina Sevrioukova for stimulating discussions. We also thank Randy J. Read for early access to his molecular replacement program.

	FOOTNOTES

^* This work was supported by National Institutes of Health Grants GM32688 (to T. L. P.) and GM25515 (to P. O. M.).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.

The atomic coordinates and the structure factors (code 1F4U, 1F4T) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).

^¶ To whom correspondence should be addressed. Tel.: 949-824-7020; Fax: 949-824-3280; E-mail: poulos@uci.edu.

Published, JBC Papers in Press, June 19, 2000, DOI 10.1074/jbc.M004281200

² Read, R. J., International Union of Crystallography Congress, Glasgow, Scotland, 1999.

	ABBREVIATIONS

The abbreviations used are: CYP119, cytochrome P450 from Sulfolobus solfataricus; P450, cytochrome P450; P450eryF, cytochrome P450 from Saccaropolyspora erythraea; P450cam, cytochrome P450 from Pseudomonas putida; bis-Tris, (bis[2-hydroxyethy]imino-tris[hydroxymethyl]methane); EPPS, (N-[2-hydroxyethyl]piperazine-N'-[3-propanesulfonic acid]; P450 BM-3, P450 from Bacillus Megaterium.

	REFERENCES

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

CiteULike

Complore

Connotea

Del.icio.us

Digg

Reddit

Technorati

This article has been cited by other articles:

				L. Li, Z. Chang, Z. Pan, Z.-Q. Fu, and X. Wang Modes of heme binding and substrate access for cytochrome P450 CYP74A revealed by crystal structures of allene oxide synthase PNAS, September 16, 2008; 105(37): 13883 - 13888. [Abstract] [Full Text] [PDF]

				M. Makino, H. Sugimoto, Y. Shiro, S. Asamizu, H. Onaka, and S. Nagano Crystal structures and catalytic mechanism of cytochrome P450 StaP that produces the indolocarbazole skeleton PNAS, July 10, 2007; 104(28): 11591 - 11596. [Abstract] [Full Text] [PDF]

				A. Verras, A. Alian, and P. R.O. d. Montellano Cytochrome P450 active site plasticity: attenuation of imidazole binding in cytochrome P450cam by an L244A mutation Protein Eng. Des. Sel., November 1, 2006; 19(11): 491 - 496. [Abstract] [Full Text] [PDF]

				D. H. Sherman, S. Li, L. V. Yermalitskaya, Y. Kim, J. A. Smith, M. R. Waterman, and L. M. Podust The Structural Basis for Substrate Anchoring, Active Site Selectivity, and Product Formation by P450 PikC from Streptomyces venezuelae J. Biol. Chem., September 8, 2006; 281(36): 26289 - 26297. [Abstract] [Full Text] [PDF]

				J. Baudry, S. Rupasinghe, and M. A. Schuler Class-dependent sequence alignment strategy improves the structural and functional modeling of P450s Protein Eng. Des. Sel., August 1, 2006; 19(8): 345 - 353. [Abstract] [Full Text] [PDF]

				B. K. Muralidhara, S. Negi, C. C. Chin, W. Braun, and J. R. Halpert Conformational Flexibility of Mammalian Cytochrome P450 2B4 in Binding Imidazole Inhibitors with Different Ring Chemistry and Side Chains: SOLUTION THERMODYNAMICS AND MOLECULAR MODELING J. Biol. Chem., March 24, 2006; 281(12): 8051 - 8061. [Abstract] [Full Text] [PDF]

				T. L. Poulos STRUCTURAL AND FUNCTIONAL DIVERSITY IN HEME MONOOXYGENASES Drug Metab. Dispos., January 1, 2005; 33(1): 10 - 18. [Abstract] [Full Text] [PDF]

				E. E. Scott, M. A. White, Y. A. He, E. F. Johnson, C. D. Stout, and J. R. Halpert Structure of Mammalian Cytochrome P450 2B4 Complexed with 4-(4-Chlorophenyl)imidazole at 1.9-A Resolution: INSIGHT INTO THE RANGE OF P450 CONFORMATIONS AND THE COORDINATION OF REDOX PARTNER BINDING J. Biol. Chem., June 25, 2004; 279(26): 27294 - 27301. [Abstract] [Full Text] [PDF]

				G. A. Schoch, J. K. Yano, M. R. Wester, K. J. Griffin, C. D. Stout, and E. F. Johnson Structure of Human Microsomal Cytochrome P450 2C8: EVIDENCE FOR A PERIPHERAL FATTY ACID BINDING SITE J. Biol. Chem., March 5, 2004; 279(10): 9497 - 9503. [Abstract] [Full Text] [PDF]

				T. L. Poulos Cytochrome P450 flexibility PNAS, November 11, 2003; 100(23): 13121 - 13122. [Full Text] [PDF]

				S. Nagano, H. Li, H. Shimizu, C. Nishida, H. Ogura, P. R. O. de Montellano, and T. L. Poulos Crystal Structures of Epothilone D-bound, Epothilone B-bound, and Substrate-free Forms of Cytochrome P450epoK J. Biol. Chem., November 7, 2003; 278(45): 44886 - 44893. [Abstract] [Full Text] [PDF]

				L. M. Podust, Y. Kim, M. Arase, B. A. Neely, B. J. Beck, H. Bach, D. H. Sherman, D. C. Lamb, S. L. Kelly, and M. R. Waterman The 1.92-A Structure of Streptomyces coelicolor A3(2) CYP154C1. A NEW MONOOXYGENASE THAT FUNCTIONALIZES MACROLIDE RING SYSTEMS J. Biol. Chem., March 28, 2003; 278(14): 12214 - 12221. [Abstract] [Full Text] [PDF]