INTRODUCTION
Fusarium head blight (FHB) is one of the most destructive and economically important crop diseases. It not only causes yield loss, but the causative agent
Fusarium graminearum (sexual stage:
Gibberella zeae) also produces a variety of mycotoxins that potentially pose a serious threat to human health and animal production (
Voigt et al., 2005). With the increase of global temperatures, FHB has spread worldwide, devastating wheat and barley crops and threatening to the global food supply by significantly reducing kernel quality (
Bai and Shaner, 2004). Assessments in nine US states showed losses of ~870 million US dollars due to FHB infection of wheat and barley from 1998 to 2000 (
Nganje et al., 2001).
Recently, several extracellular lipases secreted by
G. zeae (
Voigt et al., 2005) were reported to account for its pathogenicity.
G. zeae lipase (GZEL), which is one of these secreted lipases, is encoded by the
FGL1 gene and was demonstrated to be a virulence factor (
Voigt et al., 2005). It is also known that GZEL can be repressed by glucose and inhibited by ebelactone B, a well-known lipase inhibitor (
Voigt et al., 2005).
Lipases (triacylglycerol acylhydrolases, EC3.1.1.3) catalyze the hydrolysis of long-chain triglycerides into fatty acids and glycerol and are classified as serine hydrolases due to their inhibition by dirthylp-nitrophenyl phosphate (
Hasan et al., 2009). They also represent one of the most important biocatalyst enzyme families and have broad biotechnological applications, such as detergent industry, food production, and pharmaceuticals (
Jaeger and Reetz, 1998). Previous results demonstrate that lipase activity is greatly increased at the lipid-water interface (i.e., interfacial activation) (
Sarda and Desnuelle, 1958). Structural analyses of lipases have revealed that they share a typical α/β hydrolase scaffold comprising a core of five (or more) parallel β-strands flanked by several α-helices; their catalytic centers are formed by a highly conserved, trypsin-like, Ser-His-Asp-(Glu) motif (
Ollis et al., 1992). Another distinct structural feature of the majority of lipases is the shielding of the catalytic site from the external environment by loops or helices (variously referred to as “loops”, “flaps”, or (more often) “lids”), which lay over the catalytic center (
Brzozowski et al., 2000). According to the position of the “lids”, the reported lipase crystal structures are divided into three groups: (i) lipases with lids in a closed position (
Brady et al., 1990;
Winkler et al., 1990;
Schrag et al., 1991); (ii) lipases with open lids stabilized by a covalently bound inhibitor or detergent occluding the active site (
Brzozowski et al., 1991;
van Tilbeurgh et al., 1993;
Grochulski et al., 1994); and (iii) uninhibited lipases with open lids that are not stabilized directly by any substrate or product analog (
Schrag et al., 1997).
Here, we report the crystal structure the secreted G. zeae GZEL lipase. The unique C-terminal helix of GZEL, compared to the other reported lipases, prompted us to propose a novel “double lock” enzyme mechanism for GZEL. To gain insight into substrate/inhibitor recognition, we also propose a model of GZEL in complex with substrate and the lipase inhibitor ebelactone B that reveals well-defined sites within the catalytic cleft and suggests an interesting “anti sn-l” binding mode. These results shed light on the unique enzymatic mechanism of GZEL and provide a structural basis for anti-FHB inhibitor discovery.
RESULTS AND DISCUSSION
Overall structure of GZEL
The cDNA encoding full-length GZEL was cloned and amplified using PCR. The crystal structure consists of amino acids Met1–Ser319, with residues Gly240–Ser259 missing due to a lack of electron density that was likely caused by their flexibility in the structure. The entire GZEL molecule displays a compact globular 3D structure, with a relatively separated α-helix (α9), encompassing 50 × 40 × 40 Å (Fig. 1A and 1B). Though this accessory α-helix is quite unique in the lipase family, the GZEL core shares the typical α/β-hydrolase fold (consisting of a ten-stranded β-sheet core surrounded eight α-helices) reported for other reported lipase structures (
Ollis et al., 1992).
The anti-parallel eight strands (β1–β3 and β6–β10) in the core part are assembled as a semi-barrel β-sheet, which is surrounded by five helices (α1–α2 and α4–α6) to form the main portion and catalytic center of GZEL. The interface between the core and the C-terminal accessory α-helix α9 buries a solvent-accessible surface (SAS) (calculated by using the CCP4 package (
Collaborative Computational Project, 1994)) of 860 Å
2 compared to the 3139 Å
2 total area of α9, suggesting a strong interaction. The potential interaction surface between the GZEL core and α9 also suggest that there is a positively charged deep groove at the bottom of core, which is fully occupied by α9 in the GZEL structure (Fig. 1C). Compared with previously reported crystal structures of members of the lipase family (
Brady et al., 1990;
Winkler et al., 1990;
Kohno et al., 1996;
McAuley et al., 2004), this unique structural feature distinguishes GZEL from other lipases and indicates a putatively different enzymatic mechanism for GZEL.
An unique “lockpin-hinge-lockhole” architecture for GZEL activity
The GZEL active site is composed of Ser144, Asp198 and His257, which form the conserved catalytic triad reported in other lipase structures. Many factors have been suggested to trigger the activation of lipases at a lipid-water interface (
Brockman et al., 1973) and induce alterations of the “lid” conformation from closed to open. The “lid” of the GZEL active site, which is located at α3 helix and arranged similarly to other reported lipases, is formed by residues Ile83 to Phe94 (Fig. 3).
Using multiple sequence alignments and DALI searches, several putative key residues in the GZEL active site, both in the primary sequence and in the crystal structure, were identified as highly conserved in lipases from several different species, suggesting a similar enzymatic mechanism for GZEL (Fig. 2). To investigate this mechanism, several homologous lipase structures (
Derewenda et al., 1994;
Kohno et al., 1996;
Brzozowski et al., 2000;
McAuley et al., 2004) were carefully superimposed onto the GZEL structure based on their amino acid sequences and structural similarities. We found that the Ser144, Asp198 and His257 residues of the GZEL overlapped quite nicely onto the key residues of homologous lipase active sites with a r.m.s deviation for Cα of 0.5 Å. This indicates that, as reported for other catalytic machines, the relative positions of the key active site residues are much better conserved than the general or local protein folds (Fig. 2). Because GZEL was demonstrated to possess lipase activity (
Voigt et al., 2005), it is likely that α3, i.e., the lid, must adopt two conformations to fit its different activity states.
Although attempts to crystallize GZEL in an “open-lid” conformation failed, the conservation of GZEL among lipase family members allowed us to predict an open-lid GZEL structure by comparing the experimental lid-closed conformation to the open-lid
T. lanuginosa lipase structure (PDB code: 1EIN) (
Brzozowski et al., 2000) (Fig. 3). The proposed “open-lid” conformation was validated by PROCHECK (
Laskowski et al., 1993), and the modeled residues (Ile83–Phe94) are all in the most favored area (data not shown). The “open and closed lid structures show that the lid helix, which is formed by residues Ile83 to Phe94, is translated ~7 Å and rotated ~170° about an axis nearly parallel to that of the helix, which causes it shift away from the closed conformation to expose the active site (Fig. 3A and 3B). Moreover, the hydrophilic side of the lid fragment became partially buried in a polar cavity, which is filled by well-order solvent molecules in the closed form. Therefore, the exposure of the catalytic center is accompanied by a marked increase in the nonpolarity of the surrounding surface. Interfacial activation is thus explained by the stabilization of this nonpolar surface by the lipid environment, which would (in effect) create a catalytically competent enzyme able to attack the triglyceride molecules in the lipid phase (
Brzozowski et al., 1991).
By comparison to GZEL homologs, the GZEL crystal structure displays a unique “lockpin-hinge-lockhole” configuration, which is distinctly different from the other lipase or esterase “hinge” structures due to the additional α9 helix. Generally, there are two major types of “hinge” residues in the lipase family: (i) an Arg in the hinge position, which is present in the
T. lanuginosa lipase (
Brzozowski et al., 2000;
Yapoudjian et al., 2002); and (ii) a Gly in the hinge position, which is found in feruloyl esterases from
Aspergillus niger. In sharp contrast, however, the crystal structure of GZEL shows a novel conformational switch configuration. In the lid-closed conformation, two amino nitrogen atoms from Arg86 are hydrogen bonded with the O
δ1 atom of Asp193 and the carbon atom of Met191, with distances of 2.59 and 2.64 Å, respectively. Thus, Arg86 acts as a “lockpin” in this conformation, which is locked by the “lockhole” formed by residues Met191 to Asp193 (Fig. 4A). Furthermore, when in the fully lid-open conformation, the two hydrogen bonds formed by Arg116 and Asp192 with Met191 are broken, and the Arg86 lockpin is released. Therefore, Arg86 turns toward the solution, and the lid assumes a fully open conformation (Fig. 4B).
To analyze the GZEL substrate binding mode, the partial coordinates of oleic acid (OA), which exclude the last seven methyls (referred to as OAΔ11 herein) from the previously reported
T. lanuginose lipase and oleic acid (TLL-OA) complex structure (
Yapoudjian et al., 2002), were applied by AUTODOCK (
Huey et al., 2007) to calculate the possible binding solution (Fig. 5A and 5B). The final solution, with 186 hits over 250 runs (compared to the second solution with 64 hits), has a free energy of binding of −4.41 kcal/mol.
The crystal structure of the TLL lipase (PDB code: 1GT6) was superimposed onto the GZEL structure with a r.m.s. deviation of Cα of 2.36 Ǻ. The comparison of substrate binding behaviors in TLL and GZEL reveals a relatively similar binding mode. The opened lid in GZEL (Ile83–Phe94) and its counterpart (Asp198–Ile207) form the left and right walls of the substrate binding cavity, respectively. The side chains of Tyr89–Asn176 and Val107–Gly110 form the bottom of the cavity, with the catalytic center (Ser144–His257) at its top. This substrate binding cavity is 8 Å in width, 12 Å in height, and displays electron neutrality (Fig. 5C and 5D). The bound OAΔ11 is deeply trapped in the binding groove with its carboxylic group lying at the top of the binding cavity, which is stabilized by hydrogen bonds with the Nδ1 atom of Asn81 (3.2 Å), the Nε2 atom of H257 (3.9 Å), and the amino atom of Leu145 (3.0 Å). The fact that the alkyl chains of OAΔ11 occupy the sn-1 position in the active site indicates that the lipid binds according to the conventional mode in GZEL.
In the catalytic center, His257 occupies a position where the N
ε atom forms a hydrogen bond to the O
γ atom of Ser134 (distance of 2.5 Å), and the N
δ atom hydrogen bonds to Asp198 (distances of 2.6 Å to O
δ1 and 3.5 Å to O
δ2). This arrangement of residues is characteristic of enzymes that cleave ester bonds by nucleophilic attack and is known as a catalytic triad (
Hedstrom, 2002). His257 acts as a general base by deprotonating the serine, which allows nucleophilic attack of the carbonyl C atom of the substrate by the serine, yielding a tetrahedral transition state. The resulting positive charge on the histidine is stabilized by an electrostatic interaction with the negatively charged aspartic acid residue. The tetrahedral transition state collapses to give an acyl-enzyme intermediate. The next step is an attack of the acyl-enzyme by a nucleophilic water molecule, again activated by the histidine residue, leading to release of the acid product via a second tetrahedral intermediate transition state. The negative charge on the carbonyl O atom of the tetrahedral intermediates is stabilized by the oxyanion hole formed by the backbone N atoms of (in this case) residues Ser82 and Leu145.
These results suggest that the conformation of the GZEL lid may first change by interfacial activation, and then the lipase bond is subsequently hydrolyzed by a mechanism very similar to that used by serine proteases. Additionally, the lid may also serve as a device to inhibit the proteolytic activity of the triad, thereby protecting the enzyme itself.
Accessory lock in the lipase family
Combing the structural and functional results, Trp88 and Leu145 are indicated to be crucial for substrate recognition and binding, acting as an accessory and small “lockhole-lockpin.” In the lid-closed conformation, Trp88 is located near Leu145 (3.6 Å away), which indicates a relatively strong Van der Waals interaction that may seal the top of the binding cavity like an accessory lock (Fig. 6). In contrast, this accessory lock is released (moved a distance of nearly 8 Å) and can accommodate the carboxylic group of bound substrate in the lid-open conformation. By comparing the sequences of GZEL homologs and their structures, a similar “accessory lock” also appears in other lipases. In the multiple sequence alignment, we show that Leu145 is conserved in all G. zeae lipase homologs, while Trp88 is also largely conserved in these homologs, though identically replaced by leucine in the Aspergillus oryzae, A. niger and Feruloyl esterases. The only exception is Rhiropus niveus, which has an alanine residue instead of a Trp or Leu. All of these lipase structures were superimposed (Fig. 6B), and we found that the position of Trp88 and Leu145 were also structurally conserved. The Leu145 position is only in the locked form in the R. niveus lipase (RNL) and the closed conformation of GZEL, which is consistent with the sequence alignment results. Related substrate molecular docking was applied to OAΔ11 with different Leu145 conformations in the GZEL open conformation. From this comparison, it can be seen that bound OAΔ11 moves out of the substrate binding cavity following the conformational change of the Leu145 side chain. All of these facts suggest that this “accessory lock” may exist universally in the lipase family and be crucial for substrate binding.
Possible substrate/inhibitor binding
Because the activity of GZEL can be inhibited by ebelactone B (Fig. 7A) both
in vivo and
in vitro (
Voigt et al., 2005), we attempted to simulate the enzyme-inhibitor complex on the open-lid conformation to explain a possible inhibition mechanism (Fig. 7). We performed the molecular docking in AUTODOCK4 (
Huey et al., 2007); the best solution (53 hits in 250 runs) had a −8.42 kcal/mol free energy and was used to investigate the inhibition.
The inhibitor binding model shows that ebelactone B binds to GZEL with an interesting “
anti sn-1” binding mode, which was reported previously (
Yapoudjian et al., 2002) (Fig. 7B and 7C). Through examining the structures of ebelactone B and oleic acid, we found that the tail of ebelactone B is structurally similar to the carboxylic group of oleic acid. The last methyl occupies the carboxylic group binding position, with distances of 3.2 and 3.5 Å away from the O
γ atom of Ser144 and the N
δ atom of His257, respectively. The top tetra-lactonic group of ebelactone B is stabilized by the end of the GZEL substrate binding cavity via hydrophobic interactions. In addition, the phenyl ring of Phe94, which is parallel to the tetra-lactonic group, also stabilizes ebelactone B by a conjugated effect (Fig. 7C). These facts suggest that ebelactone B stably binds to the catalytic cavity of GZEL, which is consistent with biochemistry investigations (
Voigt et al., 2005).
CONCLUSIONS
In this study, we reported the crystal structure of GZEL with a novel “lockhole-hinge-lockpin” enzyme activity switch. Multiple sequence alignments and structure comparisons with lipases/esterases from different species show that GZEL has a classic α/β-hydrolase fold and highly conserved key residues in its active site, which indicate that GZEL functions in a similar manner as other lipases/esterases. Structural analysis revealed that the lockhole-hinge-lockpin” activity switch is formed by an accessory C-terminal α-helix, as well as Arg86 and Ile83 in the lid helix, and is unique in all reported lipase/esterase structures. Furthermore, the models of GZEL in complex with substrate and inhibitor indicate that residues Arg86 to Leu145 (or their counterparts in homologs) should function as an accessory lock in most of the reported lipase structures. Hence, we purpose that enzymes such as GZEL utilize a “double lock” mechanism. The molecular docking of a known inhibitor (
Voigt et al., 2005) also supports this model.
MATERIALS AND METHODS
Protein expression, purification, and crystallization
The
FGL1 gene was amplified by PCR, and the GZEL protein was overexpressed and subsequently purified via a previously reported method (
Voigt et al., 2005) with some modifications. Briefly, the entire
FGL1 coding sequence of was PCR-amplified from
F. graminearum cDNA. Two individual PCR fragments were independently obtained by using two sets of primer pairs: Na (5'-GATCCACCATGCGTCTCC TGTCACTCCTC-3') and Cb (5'-CTGATGAGCGGCTGGCGTGAG-3'); and Nb (5'-CACCATGCGTCTCCTGTCACTCCTC-3') and Ca (5'-AATTCTGATGAGCGGCTGGCGTGAG-3'). The fragment was then ligated into the same sites in the pLIZG7 expression vector. An EFHHHHHHP tag was fused to the C terminus of the recombinant protein. Plasmid DNA containing the
FGL1 gene was then linearized with
PmeI (NEB) and transformed into
Pichia pastoris KM71 (Invitrogen) following the manufacturer’s instructions. Transformants were grown on a histidine-deficient MD plate for 2–3 d. The positive candidates were inoculated into 3 mL BMSY, and the cultures were shaken in a 28°C incubator for 3 d until the OD
600 reached 20. Methanol was then added daily to the culture to a final concentration of 0.5% for the next 4 d to induce GZEL expression. On day 4 of induction, the cultures were harvested by centrifugation, and the activity of the supernatant was tested on olive oil/Bright Green plates (pH 7.0). The positive clones showed dark green zones. Positive supernatants also showed a protein band of the appropriate size on SDS-PAGE followed by Coomassie staining.
The culture supernatant was separated from the cells by centrifugation, and the pH of the supernatant was adjusted to 7.0. The supernatant was then filtered and applied onto a Ni-Sepharose FF column equilibrated with 25 mM Tris-HCl (pH 7.0) and 0.3 M NaCl. The target protein was eluted with a 0–1 M imidazole gradient, and fractions from the column were analyzed for activity. The fractions displaying peak enzyme activity were pooled and concentrated. The samples were then loaded onto a Superdex 75 gel-filtration column equilibrated with 25 mM Tris-HCl (pH 8.0) and 0.15 M NaCl, and the eluted active lipase was concentrated and dialyzed against 25 mM Tris-HCl (pH 8.0). The lipase was analyzed by SDS-PAGE, and pure fractions were prepared for crystallization trials.
Crystallization experiments were conducted using the vapor diffusion method (
Sun et al., 2008). Briefly, crystals were grown at 18°C from a solution containing 10 mg/mL protein, 0.2 M ammonium sulfate, 0.1 M Bis-Tris (pH 5.5), and 25% (
w/v) PEG-3350.
Data collection
The native GZEL data set was collected at 100 K using an in-house Rigaku MM-007 generator with a Mar345dtb image plate detector. Because the well solution an efficient for cyro-protectant, the fresh crystals were directly flash-cooled in a stream of nitrogen gas cooled to 100 K. The resolution of data extended to 2.8 Å, with an overall
Rmerge of 10.1%. The crystal belongs to space group
P212121, with the unit-cell parameters
a = 78.4 Å,
b = 91.0 Å,
c = 195.8 Å, and α = β = γ = 90°. We assumed four GZEL molecules per asymmetric unit with a
VM of 2.6 Å
3 Da
−1, corresponding to a solvent content of 48% (
Matthews, 1968). Processing of diffraction images and scaling of the integrated intensities were performed using the HKL2000 software package (
Otwinowski and Minor, 1997). The statistics of all data collections are presented in Table 1.
Structure determination and refinement
The structure of GZEL was solved by molecular replacement, employing the crystal structure of the lipase from
Thermomyces lanuginosa (PDBank code: 1EIN) as the initial searching model, by using the program PHASER (
McCoy et al., 2005). The clear solutions in both the rotation and translation functions indicated the presence of four molecules in one asymmetric unit, which is consistent with the Matthews coefficient and solvent content (
Matthews, 1968). Residues which differ between GZEL and the searching model were manually rebuilt in the program O (
Jones et al., 1991) under the guidance of
Fo-Fc and 2
Fo-Fc electron density maps.
After the initial refinement, a distinct density was observed in the
Fo-Fc map, which should belong to the additional Thr320–Ser349 residues of the GZEL polypeptide. The program DM (
Cowtan, 1994) was then used to perform density modification in order to refine the phases of this difference density and avoid model bias that may have arisen from phasing by molecular replacement. The remaining residues were manually built into their density in the programs O and COOT (
Jones et al., 1991;
Emsley and Cowtan, 2004) to generate the final complete coordinates. After the refinement of the model using simulated annealing, energy minimization, restrained individual B factors and addition of 221 water molecules in CNS and PHENIX (
Brunger et al., 1998;
Adams et al., 2002), the respective working
R-factor and
Rfree dropped from 0.42 and 0.45 to 0.23 and 0.26, respectively, for all data from 50.0 to 2.8 Å. Refinement was monitored by calculating
Rfree based on a subset containing 5% of the total reflections. Model geometry was verified using the program PROCHECK (
Laskowski et al., 1993). Data collection and refinement statistics are detailed in Table 1. All structure figures were prepared using PYMOL (
DeLano, 2002).
Computational methods
GZEL sequence and structural homologs were identified using BLASTp searches against the sequence database at the National Center for Biotechnology Information, and by DALI searches against the structural database (
Holm et al., 2008). A structural multiple sequence alignment of these homologous protein sequences was generated using ClustalX (
Thompson et al., 1997) and plotted with ESPript (
Gouet et al., 2003). Substrate and inhibitor docking studies were performed using AutoDock4 (
Huey et al., 2007) with MGLTools1.4.5.
PROTEIN STRUCTURE ACCESSION NUMBER
Coordinates and structure factors of the GZEL crystal structure were deposited in the Protein Data Bank under accession number 3NGM.
Higher Education Press and Springer-Verlag Berlin Heidelberg 2010