Crystal Structure of Barley 1,3-1,4-beta -Glucanase at 2.0-Å Resolution and Comparison with Bacillus 1,3-1,4-beta -Glucanase*

Jürgen J. MüllerDagger §, Karl K. Thomsen, and Udo HeinemannDagger par

From the Dagger  Forschungsgruppe Kristallographie, Max-Delbrück-Centrum für Molekulare Medizin, Robert-Rössle-Strasse 10, D-13122 Berlin, Germany, the  Department of Physiology, Carlsberg Laboratory, Gamle Carlsberg Vej 10, DK-2500 Copenhagen-Valby, Denmark, and the par  Institut für Kristallographie, Freie Universität Berlin, Takustrasse 6, D-14195 Berlin, Germany

    ABSTRACT
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Both plants and bacteria produce enzymes capable of degrading the mixed-linked beta -glucan of the endosperm cell walls of cereal grains. The enzymes share the specificity for beta -1,4 glycosyl bonds of O-3-substituted glucose units in linear polysaccharides and a similar cleavage mechanism but are unrelated in sequence and tertiary structure. The three-dimensional structure of the 1,3-1,4-beta -glucanase isoenzyme EII from barley was determined from monoclinic crystals at a resolution of 2.0 Å. The protein is folded into a beta alpha 8 barrel structure as has been shown previously (Varghese, J. N., Garrett, T. P. J., Colman, P. M., Chen, L., Høj, P. B., and Fincher, G. B. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 2785-2789) by diffraction analysis at lower resolution of tetragonal crystals. It contains one N-glycosylation site which is described in detail with the sugar moieties attached to residue Asn190. The geometry and hydration of the barley 1,3-1,4-beta -glucanase is analyzed; a model beta -glucan fragment is placed into the binding site by molecular dynamics simulation, and the beta -glucan binding grooves of the plant and bacterial enzymes are compared. Their active sites are shown to have a small number of common features in generally dissimilar geometries that serve to explain both the identical substrate specificity and the observed differences in inhibitor binding.

    INTRODUCTION
Top
Abstract
Introduction
Procedures
Results & Discussion
References

The 1,3-1,4-beta -glucanase (1,3-1,4-beta -D-glucan 4-glucanohydrolase, EC 3.2.1.73) from barley (Hordeum vulgare, cv. Alexis) degrades mixed linked beta -glucans in the cell walls of the starchy endosperm in grains (1). It acts as an endohydrolase on a long linear polysaccharide with about 30% beta -1,3 and 70% beta -1,4 glycosyl linkages (2). The scissile bond is the beta -1,4-linkage of an O-3-substituted glucose in the mixed linked beta -glucan. The crystal structure of the tetragonal form of the 1,3-1,4-beta -glucanase (isoenzyme EII)1 is known at 2.2-Å resolution as well as that of the homologous 1,3-beta -glucanase from barley (3). EII is glycosylated at Asn190 by a branched oligosaccharide of mixed composition (4). This oligosaccharide may function in increasing the thermostability of the protein (5).

1,3-1,4-beta -Glucanases occur also in bacterial strains (6). Both native and engineered enzymes have been structurally characterized (7, 8) including the 1,3-1,4-beta -glucanases MAC from Bacillus macerans and LIC from Bacillus licheniformis (9, 10). The crystal structure of H(A16-M), a hybrid beta -glucanase with residues 1-16 derived from the Bacillus amyloliquefaciens enzyme (AMY) and residues 17-214 from MAC, was solved at 1.6-Å resolution (11). The crystal structures of MAC, LIC, and H(A16-M) are sufficiently similar overall and around the active sites to suggest that these 1,3-1,4-beta -glucanases, as well as AMY and SUB, the enzyme from Bacillus subtilis, share a common mode of substrate binding and hydrolysis.

There is neither sequential nor structural homology between the plant and the bacterial enzymes, and following the universally adopted nomenclature (6), they belong to the families 17 and 16 of glycosyl hydrolases, respectively. Despite this difference they cleave the same substrate at the same cutting site and are inhibited by the same covalently binding inhibitors, 3,4-epoxyalkyl-beta -D-cellobiosides (12). However, the barley 1,3-1,4-beta -glucanase binds the epoxide preferentially with a propyl linker, whereas SUB prefers a butyl linker (13). Because both endohydrolases follow the same stereochemical pathway in glycosyl bond cleavage with retention of the beta -configuration at the anomeric carbon (14, 15), the differences in inhibitor binding are surprising at first sight and warrant a structural explanation. A cleavage mechanism with overall retention of configuration requires the presence in the catalytic site of a general acid separated by about 5.5 Å from a nucleophilic residue (6). These roles have been assigned to Glu288 and Glu232 of EII, respectively (3).

Here we present the structure analysis of monoclinic crystals of the barley 1,3-1,4-beta -glucanase (isoenzyme EII) at 2.0-Å resolution. We discuss structural differences between three molecules of EII in two space groups and compare the active center structures of plant and Bacillus 1,3-1,4-beta -glucanases. Structural reasons for inhibitor binding preferences are considered, and possible conformations and positions of a hexameric beta -glucan fragment are investigated by molecular dynamics calculations.

    EXPERIMENTAL PROCEDURES
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Diffraction Data and Structure Analysis by Molecular Replacement-- The purification, crystallization, and x-ray diffraction data of barley 1,3-1,4-beta -glucanase, isoenzyme EII, have been described (16). The diffraction data set was collected from one crystal at room temperature with a 180-mm MarResearch imaging plate system on an Enraf Nonius FR571 rotating anode x-ray generator (45 kV, 90 mA). The images were evaluated with MOSFLM, AGROVATA, and ROTAVATA from the CCP4 suite (20). Relevant crystallographic parameters in the monoclinic space group P21 are summarized in Table I. Space group and unit cell dimensions are consistent with the presence of two protein molecules per asymmetric unit. The 40,727 unique reflections correspond to 99.3% of the observations expected for 13.65 to 2-Å resolution, and the completeness in the outermost resolution shell from 2.09 to 2.0 Å is 96.6%. The crystal and x-ray diffraction data for EII in the tetragonal space group P43212 (Ref. 3; PDB entry 1ghr) are added to Table I for comparison.

                              
View this table:
[in this window]
[in a new window]
 
Table I
Crystallographic data

Since the structure of EII is known for the space group P43212, the phases of the structure factors could be determined by molecular replacement using XPLOR (17). Two molecules in the asymmetric unit (Mol1 and Mol2) were found by rotation and translation search. The search model 1ghr (3) gave two single peaks 10sigma and 9sigma above the average in the Patterson correlation function. After translation search (15.7sigma and 14.8sigma ), rigid body, and packing refinement against data between 13.65 and 3 Å the R value was 31.4%.

Structure Refinement-- For refinement, 5% of the reflections were set aside as test set for the calculation of Rfree (18). Rfree was used to monitor convergence throughout. After 100 steps of minimization with the standard repel nonbonded energy function followed by 50 steps of Powell Lennard-Jones minimization (Rwork = 28.3%, Rfree = 32.4%), a simulated annealing run following the standard slow-cooling protocol (19) with structure amplitudes F > 2sigma (F) between 13.65 and 2.0 Å, restrictions for non-crystallographic symmetry (NCS) converged with Rwork = 27.2% and Rfree = 31.2% (no NCS: Rwork = 26.3%, Rfree = 31.1%). In the following steps, NCS restrictions were used for all but those residues making crystallographic or non-crystallographic intermolecular contacts. Subsequent B-value refinement resulted in Rwork = 25.2% and Rfree = 28.5%. Electron density maps were calculated with XPLOR and CCP4 (20) programs and displayed in O (21). Further refinement was performed with XPLOR/O combining positional refinement with overall and atomic B-value refinement. After several cycles of positional and B-value refinement and manual revision of the model via electron density and difference density maps, water molecules were added until most significant peaks (>4sigma ) disappeared from the difference electron density map. Ending up with 418 water molecules in the asymmetric unit and reduced NCS restraints for both molecules, the final Rwork was 17.1% and Rfree 21.2%.

Molecular Dynamics Modeling of Substrate Binding-- Substrate binding to a bacterial 1,3-1,4-beta -glucanase has been modeled (11). Because the substrate specificity and the stereochemical course of beta -glucan hydrolysis are identical for Bacillus and barley 1,3-1,4-beta -glucanases, we used the hexameric beta -glucan fragment proposed earlier (11) as start model in simulating substrate binding to EII. The glucose moieties are numbered from 1 to 6 beginning at the non-reducing end. They bind to enzyme subsites p4, p3, p2, p1, p1', and p2', and cleavage occurs at the catalytic site between p1 and p1'. Manual primary positioning of the hexamer in the binding cleft using O (21) took into account restraints deduced from chemical probing with barley 1,3-1,4-beta -glucanase (22) and inhibitor binding to the structurally very similar 1,3-beta -glucanase of barley (1ghs; Ref. 23).

To start the molecular dynamics simulation, a FORTRAN program written for this purpose but not specific for the glucanase molecule first filled a P1 cell with dimensions 66.379 × 60.241 × 52.643 Å with TIP3 waters provided by the XPLOR package (24). After this procedure, molecule Mol1 of EII was oriented with its inertia equivalent ellipsoid axes parallel to the P1 cell axes. The overlapping water molecules were removed so that Mol1 finally was embedded in 5,069 non-overlapping water molecules. The water box is large enough to cover the substrate by a water shell of at least 13-Å thickness. The charged residues of the protein are made net neutral to mimic the effects of solvation and counterions. The nonbonded interaction cutoff was specified to 8.5 Å. The molecular dynamics calculation, done by XPLOR, preceded by 120 steps of standard repel non-bonded energy minimization and 80 steps of Powell Lennard-Jones minimization in vacuum, started with the fixed enzyme-substrate complex to relax the water molecules around the macromolecule. The solvent molecules were relaxed further during 120 steps of energy minimization followed by molecular dynamics over 5 ps at 100 K, 5 ps at 200 K, 20 ps at 300 K with 0.2-ps reassignment using the CHARMM force field. The energy was then minimized for the solute (100 steps Powell) followed by several steps of molecular dynamics for the solute (5 ps at 100 K, 5 ps at 200 K, 5 ps at 300 K). Finally, molecular dynamics were run for 80 ps at 300 K (with temperature coupling) for the whole system, followed by 80 steps of Powell Lennard-Jones minimization. Over the entire simulation the EII molecule remained fixed, but those side chains in the binding cleft containing atoms in a 4-Å surface shell, the hexaglucan, and the waters were allowed to move. Alternative protocols without water relaxation at 100, 200, and 300 K and without temperature coupling resulted in similar sugar conformations and positions.

    RESULTS AND DISCUSSION
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Refinement Results-- After combining working and test data sets, final positional and B-value refinement yielded an R-value of 17.0%. The R-value was 16.5% for all observations with a low resolution cutoff at 6 Å as used in the previous structure analysis (3). Two residues of N-acetyl-D-glucosamine (NAG) per molecule of the barley 1,3-1,4-beta -glucanase EII were located at the glycosylation site, Asn190. The side chains of Asp183 and Ile304 in Mol1 and Gln41 in Mol2 were refined in two alternative positions with 2/3 and 1/3 occupancy, and two acetate molecules were identified in non-related positions. In the Ramachandran diagram, more than 91% of the non-glycine and non-proline residues are within the most favored regions, and all other residues are within additional allowed regions. Two cis-peptides preceding Pro137 and Ala276, identical to those in 1ghr, are present in both Mol1 and Mol2. The average correlation of the calculated Fc map with the 2Fo - Fc electron density map is better than 0.9 for both molecules according to an analysis with O (21). The correlation coefficient drops below 0.7 for the highly flexible side chains of Arg100, Arg197, and Arg261, which have B values of about 50 Å2 and are not sterically restricted by hydrogen bonds or van der Waals contacts in both molecules.

Stereochemical parameters of the main and side chains are strongly restrained by standard weights and are better than or inside of the bandwidth defined in PROCHECK (25) for structures with comparable resolution. The final atomic coordinate set representing the monoclinic form of EII was scrutinized with WHATCHECK (26), PROCHECK, and several programs of the CCP4 suite. The results are summarized in Tables I and II. The experimental data and the refined atomic coordinates of both molecules in the asymmetric unit were submitted to the Brookhaven Protein Data Bank (27) entry code 1aq0.

                              
View this table:
[in this window]
[in a new window]
 
Table II
Structure refinement

Overall Structure and Crystal Packing-- The overall structure and the active site of EII have been described by Varghese et al. (3) and are very similar in the monoclinic crystal form. The overall shape of the 1,3-1,4-beta -glucanase can be approximated by its inertia equivalent ellipsoid with half-axes of 29.1 × 25.2 × 16.4 Å. The global folding pattern belongs to the beta alpha 8 barrel type, although beta -strand number 8 of the barrel is truncated to just two residues (Fig. 1). The active site is located in an open cleft at the bottom of the barrel defined by the C-terminal ends of the parallel intra-barrel beta -strands. It is about 36 Å long and 8 to 9 Å deep, allowing the binding of oligosaccharide substrates.


View larger version (61K):
[in this window]
[in a new window]
 
Fig. 1.   Stereographic cartoon of the molecular structure of the barley 1,3-1,4-beta -glucanase, isoenzyme EII. The view is through the active site channel of the enzyme. Arrows indicate beta -strands, and wound ribbons indicate alpha -helices. No arrow is drawn for beta -strand number 8 consisting of only two residues. The two catalytic glutamic acid residues 232 and 288 of EII are drawn in ball-and-stick representation, as well as the carbohydrate attached to residue Asn190. Drawn with MOLSCRIPT (39).

The two molecules of EII in the asymmetric unit make numerous intermolecular contacts with distances between 2.5 and 4 Å. 13 amino acids of Mol1 are in contact with 14 amino acids of Mol2 through non-crystallographic symmetry. 39 and 43 amino acids of Mol1 and Mol2, respectively, engage in additional crystallographic contacts. These contacts result from a tighter packing of the molecules in the monoclinic cell (Vm = 2.41 Å3 Da-1) in comparison with tetragonal 1ghr (Vm = 3.26 Å3 Da-1), where 44 residues make intermolecular contacts. The contacts do not include any residue of the catalytic site and of the potential substrate-binding site (see below) of Mol1 and Mol2, as well as of 1ghr, and therefore do not exert a major influence on the protein structure within the binding cleft. This is suggested by the analysis of distances (LSQKAB, Ref. 28) between 45 equivalent residues within the 4-Å surface of the substrate binding cleft. A comparison of Mol1, Mol2, and 1ghr shows the residues in the cleft to be structurally conserved to a higher degree than those in the rest of the molecules (Table III). The differences between the atom positions of the residues in the catalytic and substrate-binding site (Glu288, Glu232, Tyr33, Glu93, Asn92, Val134, Phe135, Asn168, Tyr170, Leu173, Phe275, Glu280, Lys283, and Trp291) are below the average coordinate error for Mol1 and Mol2 determined to 0.2 Å from the Luzzati (29) analysis.

                              
View this table:
[in this window]
[in a new window]
 
Table III
Statistics measuring the degree of similarity of molecules Mol1, Mol2, and 1ghr

Structural Heterogeneity, Subdomains, and Hydration-- The most prominent difference between the main chains of the three molecules is found for amino acid residues 190-200 of Mol2 compared with the other two (Fig. 2). In this region, the structure of Mol2 is influenced by the NCS interactions of Gly199 and Ala200 in a beta -turn of Mol1 with Val196 in Mol2, and by numerous crystallographic contacts to Mol2 itself. The shift of the Calpha backbone preceding the small beta -sheet by a maximum of 4.8 Å at Thr194 indicates a certain flexibility of this region. Statistics measuring the degree of similarity of molecules (30) show clearly that Mol1 is more similar to 1ghr than is Mol2 (Table III). Therefore, all considerations of general aspects of the structures will refer to Mol1.


View larger version (29K):
[in this window]
[in a new window]
 
Fig. 2.   Least-squares superpositions of the two molecules of EII related by non-crystallographic symmetry in the P21 crystal and 1ghr, the tetragonal form of the enzyme. The view is onto the beta alpha 8 barrel. The Calpha traces are shown as bold solid line for Mol1, dashed line for Mol2, and thin solid line for 1ghr (3). The 1ghr trace is almost completely hidden behind the solid line. Note the two nearly perpendicular depressions on this face of the molecule, the beta -glucan-binding channel running horizontally and the domain boundary crossing in vertical orientation (see text).

By using the program PUU (31) two structural domains were detected with a boundary crossing the beta alpha 8 barrel perpendicular to the binding cleft. Residues Ile1 to Ser126 and Tyr273 to Phe306 constitute domain I; residues Val127 to Thr272 belong to the glycosylated domain II (right side in Fig. 3). Only two strands link both domains at Ser126/Val127 and Thr272/Tyr273. This assignment of structural domains differs from that proposed earlier (3) and is supported by a least-squares superposition of domains I of Mol1 and 1ghr yielding a rotation by about 1° between their domains II. The rotation axis traverses the center of mass of domain II and is inclined by about 50° against the cleft axis marked by a modeled hexaglucan substrate (see below) in Fig. 3. Interestingly, the two beta -strands crossing the inter-domain boundary are not regularly structured. Strand beta 8 is truncated to just two residues in Mol1 and Mol2 and is next to the cis-peptide preceding Ala276. The long strand beta 5 is interrupted by a 310 turn at Gln129 to Ile131 where it crosses from domain I into domain II also reflecting some structural irregularity. Functional implications of this domain structure will be discussed below.


View larger version (146K):
[in this window]
[in a new window]
 
Fig. 3.   Domain structure of barley 1,3-1,4-beta -glucanase. The domains were estimated by the program PUU (31). The red net defines the domain boundary calculated by ANY_PLANE (JJM) and CURVE (32). A hexameric beta -glucan model oriented such that the reducing end points to the right is shown together with the catalytic site residues Glu232 and Glu288, Tyr33, and the glycosylation site with NAG11 and NAG21. Presumably, the catalytic event is coupled with a movement of the left domain against the glycosylated right domain, thereby distorting the scissile bond of the beta -glucan positioned directly at the domain boundary. Drawn with SETOR (40).

Water molecules in 130 pairs related by NCS were identified allowing for a maximal separation of 1.5 Å. These water-binding sites present in both Mol1 and Mol2 are expected to have biological relevance. 37 out of these 130 sites have identical positions within 1.5 Å with any of the 48 water sites present in 1ghr indicating that their positions are independent of crystal contacts and preparation conditions. The water distribution around EII is quite asymmetric, since the molecular surface near the C and N termini shows significantly higher hydration than other parts of the protein molecule. This finding is true also for the waters detected in the molecule 1ghr. Interestingly, only 12 waters out of the 130 are found within the substrate-binding cleft of EII. The lack of complete cleft hydration may be due to the presence of exposed hydrophobic residues (see below).

The Glycosylation Site-- Barley 1,3-1,4-beta -glucanase, isoenzyme EII, has only one site for potential N-glycosylation, the sequence Asn190-Ala191-Ser192. The attached carbohydrate moieties were analyzed by Harthill and Thomsen (4). They identified five different branched N-glycans comprised of different sugars contributing with different relative amounts between 15 and 30%. All have a common core sequence starting with two beta -1,4-linked NAG molecules. These two residues are represented by clear electron density (Fig. 4). They have been identified in difference density maps at a late stage of refinement, because they were absent in 1ghr. The sugars in the two NCS-related protein molecules are differently well defined. The average temperature factors are 30.0 and 60.7 Å2 for NAG11 and NAG21 of Mol1, and 49.1 and 63.6 Å2 for the corresponding residues in Mol2. Whereas NAG11 and NAG21 are well defined, a third sugar residue, alpha -L-fucose, attached to O-3 of NAG12 is seen at very low electron density (0.3 e Å-3) but was not modeled. Both glycosylation sites are sterically accessible to different degrees.


View larger version (76K):
[in this window]
[in a new window]
 
Fig. 4.   Stereographic drawings of electron density around the carbohydrate moiety attached to Asn190. The 2Fo - Fc omit map with a cushion of 6 Å around the carbohydrate is contoured at 0.62 e Å-3 around NAG11 (top) and at 0.34 e Å-3 around NAG21 (bottom). Drawn with OPLOT (21).

Both disaccharides are in the energetically preferred 4C1 chair conformation, and the torsions of the beta -1,4 glycosyl bonds Phi (H-1-C-1-O-4-C-4) and Psi (C-1-O-4-C-4-H-4) are 49.8° and 15.2° in NAG11/NAG21, and 44.9° and 3.9° in NAG12/NAG22. These dihedrals correspond to Phi  and Psi  angles also found in saccharides bound to lysozyme (PDB entry 1lzr) or influenza virus neuraminidase (PDB entry 1nnc).

The Active Site: Structure and Hydration-- The substrate binding cleft of barley 1,3-1,4-beta -glucanase was geometrically defined using SURFNET (32), and residue accessibilities were characterized by NACCESS (33). The largest cleft found by SURFNET (1 Å < sphere radius < 4 Å; cutoff distance between atoms and mask region 4 Å) has a volume of 2031 Å3 (Fig. 5, top) and comprises the catalytic site and the potential substrate binding site as described recently (3, 23). The global dimensions are 7.5 Å in width, 8 to 9 Å in depth around the catalytic center, and 36 Å in length. All accessible atoms within a distance <= 4 Å from the cleft surface were selected, and their polarity was estimated by NACCESS. The resulting patterns are given in Table IV. The "left" side of the cleft (in the orientation of Fig. 5) is mainly decorated with apolar residues, whereas the "right" side, especially in the vicinity of the catalytic site, is covered with polar residues. The residues Tyr170, Glu232, Tyr33, Glu288, and Trp291 in or near the catalytic center are accessible, also Leu173, Phe275, Asn92, Glu93, and Asn168 which possibly interact with the substrate (3) as partially proven by chemical probing (22).


View larger version (137K):
[in this window]
[in a new window]
 
Fig. 5.   Stereographic drawings of the substrate binding cleft of barley 1,3-1,4-beta -glucanase and of the hybrid Bacillus 1,3-1,4-beta -glucanase H(A16-M). Top, barley 1,3-1,4-beta -glucanase, active site atoms of Glu232, Tyr33, and Glu288 are labeled. Bottom, H(A16-M). The active site residues Glu105, Tyr94, and Glu109 are labeled. Water molecules binding to the corresponding cleft are shown as blue spheres. The maps defining the clefts were estimated by SURFNET (32). Drawn with BOBSCRIPT (42) and RASTER3D (43).

                              
View this table:
[in this window]
[in a new window]
 
Table IV
Active site residues in barley and Bacillus 1.3-1.4-beta -glucanases
Active site residues of EII and H(A16-M) possibly interacting with bound substrate. Left, and right is according to Fig. 5. Attached water molecules are indicated by w; a denotes the apolar part of the residue; and p refers to the polar part. The catalytic site glutamic acids and the tyrosine defining subsite p2 are emphasized by boldface.

Water molecules are mainly located at the polar right side of the cleft floor (Fig. 5) where several residues form hydrogen bonds to five water molecules. Only the three water molecules at the lowest end of the cleft are also present in the 1ghr model (3). Three water molecules form H bonds to the Ala174 backbone and the Tyr172 and Glu280 side chains at the left side. Thus, the mostly polar right half of the cleft is covered by a water layer. A comparison with the water structure in the active site of the Bacillus 1,3-1,4-beta -glucanase will follow below.

Molecular Dynamics Simulation of Substrate Binding-- A model of a hexaglucan substrate bound to the active site of EII was constructed to assist our understanding of the enzymology of the barley 1,3-1,4-beta -glucanase. The beta -glucan before and after the simulated docking procedure is shown in Fig. 6 (bottom). In the start configuration, C-1 of glucose residue Glc4 points to the nucleophile Glu232 of EII; Glc4 is situated halfway between Glu232 and the general acid Glu288; Glc3 O-6 forms hydrogen bonds to Glu280 and Lys283 (as the inhibitor in 1ghs); Glc3 O-2 makes an H bond to Asn92 (as in 1ghs), and the face of Glc3 is in hydrophobic contact with Tyr33 (as in 1ghs). The water molecules localized in the cleft were removed completely, because no fixed waters could be detected between inhibitor and 1,3-1,4-beta -glucanase in the crystal structure of an enzyme-inhibitor complex (34).


View larger version (99K):
[in this window]
[in a new window]
 
Fig. 6.   Structural comparison of barley (left) and Bacillus 1,3-1,4-beta -glucanases (right). Top, the EII folding pattern belongs to the beta alpha 8 barrel type; H(A16-M) adopts the jellyroll fold. The key residues for the catalytic event are marked in both enzymes. Middle, molecular surface drawings of the 1,3-1,4-beta -glucanases. The hexameric beta -glucan models are oriented such that the reducing ends point upwards and shown together with the catalytic site residues Glu232 and Glu288 (EII, left) and Glu105 and Glu109 (H(A16-M), right). Tyrosine residues Tyr33 (EII) and Tyr94 (H(A16-M)) stacking against Glc3 of the substrate at subsite p2 are also shown. Surface colors indicate electrostatic potential (blue, positive, and red, negative) as calculated with DELPHI (41). Note the different rotational setting of the substrate chain in the active site cleft and the different orientation of the protein side chains. Bottom, superposition of beta -glucan model substrates for barley and hybrid Bacillus 1,3-1,4-beta -glucanases. The molecules are superimposed with the glucose moiety Glc3 interacting with a tyrosine in subsite p2. The substrate strand fitting into the active site of H(A16-M) (Ref. 11, green) defines the start conformation used in molecular dynamics modeling of the substrate bound to EII (red).

The molecular dynamics simulation of substrate binding was equilibrated after 70 ps. For further discussion we use the middle structure of the inhibitor as determined by LSQMAN (35). Sugar residues Glc2, Glc3, Glc4, and Glc5 are situated within the binding cleft, their root mean square distances from the middle structure are 0.38, 0.31, 0.72, and 0.66 Å, respectively. The first residue, Glc1, is outside of the binding cleft as defined by SURFNET (Fig. 5). The root mean square distances of Glc1 and Glc6 are about 0.66 and 0.68 Å, respectively. Because the position of the substrate is strongly influenced by Tyr33 and the contacts are built in manually, most of the interactions found in the inhibitor complex of the 1,3-beta -glucanase (1ghs; Ref. 3) are present in this model, too. The global course of the hexamer is visible in Fig. 3. The smallest root mean square Delta value for Glc3 coincides with the longest retention of any glucose residue in one position during the simulation. This is in agreement with the results of the subsite mapping of barley 1,3-beta -glucanase isoenzyme GII (36), where the binding affinity at the second subsite toward the non-reducing end relative to the cutting point is at the maximum. In the Bacillus as well as in the barley 1,3-1,4-beta -glucanases the first and sixth residues of the beta -glucan model substrates are not tightly bound to the ends of the clefts.

All glucose units of the modeled substrate are in the preferred 4C1 chair pucker, and the torsion angles are in the low energy regions between 28 and 54° for Phi  and -52 to +14° for Psi  in beta -1,3 glycosyl bonds (37), and around Phi  = 40°, Psi  = -20° for beta -1,4-bonds (38). Only the torsion angle Phi  = 72.5° at the scissile beta -1,4-bond is outside the minimum, indicating strain which may be relieved by bond cleavage. A deviation from the low energy conformation of similar magnitude but with opposite sign was observed in the modeled beta -glucan substrate of the Bacillus 1,3-1,4-beta -glucanase (11).

In rationalizing substrate binding and conversion, the domain structure and possible structural rearrangements must be taken into account. The general acid, Glu288, and the substrate fix point Tyr33 are located at the same domain I, whereas the nucleophile Glu232 at the bottom of the cleft belongs to the glycosylated domain II. The substrate residues Glc1 to Glc4 are attached to domain I, and the scissile bond O-4 of the modeled hexamer is positioned at the inter-domain boundary, and residues Glc5 and Glc6 are attached to domain II (Fig. 3). A slight rotation of the domains relative to each other may support a rearrangement of the catalytic residues and/or the substrate position, thus providing a fine tuning of the catalytic event. A domain rearrangement of this kind is observed between Mol1 and Mol2 of EII (see above).

Comparison of Active Site Structures and Modeled Substrate Binding by Barley and Bacillus 1,3-1,4-beta -Glucanases-- Active site directed inhibition of 1,3-1,4-beta -glucanases from B. subtilis and barley (13) by covalent modification with epoxyalkyl cellobiosides reveals subtle differences in their active site geometries. Both enzymes with unrelated primary, secondary, and tertiary structures (Fig. 6) have identical substrate specificity, but the optimum aglycon length of the cellobiosides for maximum inactivation is C-4 (SUB) and C-3 (EII), respectively, and SUB is much more efficiently inhibited by the (S)-epoxybutyl-beta -cellobioside than by the R-isomer, whereas the reverse is found for EII. From this, it was concluded that the differences in inhibitor binding might be related to a different cleavage mechanism (13). It could be shown recently, however, that beta -glucan hydrolysis follows the same stereochemical course in the bacterial and plant enzymes (14, 15). In addition, the covalent complex of the barley 1,3-beta -glucanase with 2,3-epoxypropyl cellobioside shows binding of the S-isomer to the nucleophile of this close homolog of EII.

A common motif in all inhibitor binding studies (14, 34) and substrate binding models (Ref. 11; this work) is the hydrophobic "stacking" interaction of the sugar residue at the p2 subsite with a tyrosine side chain (Tyr33 in EII and Tyr94 in H(A16-M)). The distances of the plane midpoints are between 4.3 and 4.7 Å in all complexes. In H(A16-M) (2ayh), as well as in the barley glucanases (Mol1 of EII, 1ghr, 1ghs), the distances between oxygens Oepsilon -1 and Oepsilon -2 of the catalytic nucleophile and the corresponding tyrosine midpoint are identical (about 4.5 and 6.5 Å). The distances between the active Oepsilon atoms and the tyrosine midpoint are changed to 5.7 Å by the inhibitor binding to H(A16-M) and to 5.4 Å by substrate binding modeled for EII (Mol1). Taking this distance of about 5.5 Å into account, the additional length of 1.3 Å of the butyl linker in comparison to the propyl linker requires a bulging out of the linker. The molecular environment of the butyl linker of the inhibitor bound to H(A16-M) consists of residues Trp192, Val88, Trp184, and Phe30 at one side, and of Tyr123, Trp103, and Phe92 at the other. The butyl linker interacts by hydrophobic contacts with the plane of Phe92 and covers the hydrophobic surface. Such a large hydrophobic area does not exist in the case of EII. The only hydrophobic residues in the vicinity of the linker are Phe275 and Trp291, and therefore, preserving the strong hydrophobic interaction between sugar residue Glc3 and Tyr33, the 1.3 Å shorter propyl linker of the epoxyalkyl cellobioside inhibitor is preferentially bound.

By analogy to the cleft determination and characterization described above for EII, the substrate binding site of H(A16-M) was investigated. Here, the form of the cleft is quite different from that of EII (Fig. 5, bottom), and the volume is about 1,640 Å3. The width is comparable with the cleft width of 7-8 Å of EII near the catalytic site, but the groove is considerably deeper (about 12.5 Å) and shorter (about 29 Å). In both cases the clefts are somewhat branched out at the molecular surface (Fig. 5). All residues within the active site (Asp107, Glu105, Glu109, Trp103, Glu119, and Trp192) or interacting with the cellobioside inhibitor (Tyr24, Asn26, Glu63, and His99) as described (34) and additional potential substrate binding partners (Arg65, Tyr94, Ser90, Glu131, and Asn121) have accessible surface fractions. These residues are summarized in Table IV for comparison with EII. No obvious similarity between both patterns is found. Whereas the branched cleft in EII has a nearly isometric cross-section in the substrate binding region, the cleft of H(A16-M) is narrow and deep. The hydration pattern of the H(A16-M) groove is quite different from the water structure in the EII cleft (Fig. 5).

The six glucose residues cover the whole length of the cleft in H(A16-M), but eight can be accommodated in EII (36). In any case, only the innermost four appear tightly bound within the reaction center. In Fig. 6 (middle), the surfaces, modeled substrates, catalytic residues, and the hydrophobic binding centers at subsite p2 (Tyr33 of EII and Tyr94 of H(A16-M)) are shown. The basic differences between both families of 1,3-1,4-beta -glucanases are the relative positioning of the general acid, the nucleophile, and these tyrosine residues. Given the same orientation of the substrate with the reducing end of Glc6 facing upward in Fig. 6 (middle), the general acid, Glu288, of EII is at the left side of the hexameric beta -glucan, but the corresponding Glu109 is at the cleft bottom for H(A16-M). To follow the same stereochemical pathway a global rotation is necessary between both substrates by about 90°.

In conclusion, the active sites of the plant and Bacillus 1,3-1,4-beta -glucanases are surprisingly different in view of their close functional similarity. The enzymes show nicely how the same catalytic activity can evolve on completely different protein folds and in dissimilar local geometries.

    ACKNOWLEDGEMENTS

Kim Henrick is thanked for helpful discussions regarding the carbohydrate moieties and Yves Muller for critically reading the manuscript.

    FOOTNOTES

* This work was supported by the Deutsche Forschungsgemeinschaft through He 1318/9-3 and the Fonds der Chemischen Industrie.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 structure factors (code 1AQ0) have been deposited in the Protein Data Bank, Brookhaven National Laboratory, Upton, NY.

§ To whom correspondence should be addressed. Tel.: 49-30-94063421; Fax: 49-30-94062548; E-mail: jjm{at}mdc-berlin.de; http:// www.mdc-berlin.de/~crystal.

1 The abbreviations used are: EII, 1,3-1,4-beta -glucanase isoenzyme EII from barley; Mol1, Mol2, molecules 1 and 2 from the asymmetric unit of monoclinic EII crystals, respectively; AMY, MAC, LIC, SUB, 1,3-1,4-beta -glucanases from B. amyloliquefaciens, B. macerans, B. licheniformis, and B. subtilis, respectively; H(A16-M), hybrid beta -glucanase with residues 1-16 derived from AMY and residues 17-214 from MAC; NAG, N-acetyl-D-glucosamine; PDB, Brookhaven Protein Data Bank; 1ghr, 1ghs, 2ayh, 1byh, PDB entries for the tetragonal crystal form of EII, 1,3-beta -glucanase isoenzyme GII from barley, H(A16-M), and H(A16-M) with 3,4-epoxybutyl beta -D-cellobioside inhibitor, respectively; NCS, non-crystallographic symmetry.

    REFERENCES
Top
Abstract
Introduction
Procedures
Results & Discussion
References

  1. Fincher, G. B., Lock, P. A., Morgan, M. M., Lingelbach, K., Wettenhall, R. E. H., Mercer, J. F. B., Brandt, A., Thomsen, K. K. (1986) Proc. Natl. Acad. Sci. U. S. A. 83, 2081-2085[Abstract]
  2. Bamforth, C. W. (1982) Brew. Dig. 57, 22-35
  3. Varghese, J. N., Garrett, T. P. J., Colman, P. M., Chen, L., Høj, P. B., Fincher, G. B. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 2785-2789[Abstract]
  4. Harthill, J. E., and Thomsen, K. K. (1995) Plant Physiol. Biochem. 33, 9-18
  5. Doan, D. N. P., and Fincher, G. B. (1992) FEBS Lett. 309, 265-271[CrossRef][Medline] [Order article via Infotrieve]
  6. Davies, G., and Henrissat, B. (1995) Structure 3, 853-859[Medline] [Order article via Infotrieve]
  7. Heinemann, U., Ay, J., Gaiser, O., Müller, J. J., Ponnuswamy, M. N. (1996) Biol. Chem. Hoppe-Seyler 377, 447-454
  8. Heinemann, U., and Hahn, M. (1996) Trends Biochem. Sci. 20, 349-350[CrossRef]
  9. Hahn, M., Olsen, O., Politz, O., Borriss, R., and Heinemann, U. (1995) J. Biol. Chem. 270, 3081-3088[Abstract/Free Full Text]
  10. Hahn, M., Pons, J., Planas, A., Querol, E., and Heinemann, U. (1995) FEBS Lett. 374, 221-224[CrossRef][Medline] [Order article via Infotrieve]
  11. Hahn, M., Keitel, T., and Heinemann, U. (1995) Eur. J. Biochem. 232, 849-858[Abstract]
  12. Høj, P. B., Rodriguez, E. B., Stick, R. V., Stone, B. A. (1989) J. Biol. Chem. 264, 4939-4947[Abstract/Free Full Text]
  13. Høj, P. B., Rodriguez, E. B., Iser, J. R., Stick, R. V., Stone, B. A. (1991) J. Biol. Chem. 266, 11628-11631[Abstract/Free Full Text]
  14. Chen, L., Sadek, M., Stone, B. A., Brownlee, R. T. C., Fincher, G. B., Høj, P. B. (1995) Biochim. Biophys. Acta 1253, 112-116[Medline] [Order article via Infotrieve]
  15. Malet, C. J.-B., Bernabé, M., Brosa, C., and Planas, A. (1993) Biochem. J. 296, 753-758[Medline] [Order article via Infotrieve]
  16. Keitel, T., Thomsen, K. K., and Heinemann, U. (1993) J. Mol. Biol. 232, 1003-1004[CrossRef][Medline] [Order article via Infotrieve]
  17. Brünger, A. T. (1990) Acta Crystallogr. Sect. A 46, 46-57[CrossRef]
  18. Brünger, A. T. (1992) Nature 355, 472-474[CrossRef]
  19. Brünger, A. T., Krukowski, A., and Erickson, J. W. (1990) Acta Crystallogr. Sect. A 46, 585-592[CrossRef][Medline] [Order article via Infotrieve]
  20. Collaborative Computational Project, No. 4. (1994) Acta Crystallogr. Sect. D 50, 760-763 [CrossRef][Medline] [Order article via Infotrieve]
  21. Jones, T. A., Zou, J.-Y., Cowan, S. M., Kjeldgaard, M. (1991) Acta Crystallogr. Sect. A 47, 110-119[CrossRef][Medline] [Order article via Infotrieve]
  22. Chen, L., Fincher, G. B., and Høj, P. B. (1993) J. Biol. Chem. 268, 13318-13326[Abstract/Free Full Text]
  23. Chen, L., Garrett, T. P. J., Fincher, G. B., Høj, P. B. (1995) J. Biol. Chem. 270, 8093-8101[Abstract/Free Full Text]
  24. Brünger, A. T. (1992) X-PLOR Manual, Version 3.1, Yale University Press, New Haven
  25. Laskowski, R. A., MacArthur, M. W., Moss, D. S., Thornton, J. M. (1993) J. Appl. Crystallogr. 26, 283-291 [CrossRef]
  26. Vriend, G. (1990) J. Mol. Graph. 8, 52-56[CrossRef][Medline] [Order article via Infotrieve]
  27. Bernstein, F. C., Koetzle, T. F., Williams, G. J. B., Meyer, E. E., Jr., Brice, M. D., Rodgers, J. R., Kennard, O., Shimanouchi, T., Tasumi, M. (1977) J. Mol. Biol. 112, 535-542[Medline] [Order article via Infotrieve]
  28. Kabsch, W. (1976) Acta Crystallogr. Sect. A 32, 922-923[CrossRef]
  29. Luzzati, V. (1952) Acta Crystallogr. 5, 802-810[CrossRef]
  30. Kleywegt, G. J. (1996) Acta Crystallogr. Sect. D 52, 842-857 [CrossRef]
  31. Holm, L., and Sander, C. (1994) Proteins Struct. Funct. Genet. 19, 256-268 [Medline] [Order article via Infotrieve]
  32. Laskowski, R. A. (1995) J. Mol. Graph. 13, 323-330[CrossRef][Medline] [Order article via Infotrieve]
  33. Hubbard, S. J., and Thornton, J. M. (1993) NACCES, version 2.1, Computer Program, Department of Chemistry and Molecular Biology, University College London
  34. Keitel, T., Simon, O., Borris, R., and Heinemann, U. (1993) Proc. Natl. Acad. Sci. U. S. A. 90, 5287-5291[Abstract]
  35. Kleywegt, G. J., and Jones T. A. (1997) Methods Enzymol., in press
  36. Høj, P. B., and Fincher, G. B. (1995) Plant J. 7, 367-379[CrossRef][Medline] [Order article via Infotrieve]
  37. Noguchi, K., Okuyama, K., Kitamura, S., and Takeo, K. (1992) Carbohydrate Res. 237, 33-43[Medline] [Order article via Infotrieve]
  38. Sathyanarayana, B. K., and Rao, V. S. R. (1971) Biopolymers 10, 1605-1615
  39. Kraulis, J. P. (1991) J. Appl. Crystallogr. 5, 802-810 [CrossRef]
  40. Evans, S. V. (1993) J. Mol. Graph. 11, 134-138[CrossRef][Medline] [Order article via Infotrieve]
  41. Nicholls, A., Sharp, K. A., and Honig, B. (1991) Proteins Struct. Funct. Genet. 11, 281-296 [Medline] [Order article via Infotrieve]
  42. Esnouf, R. M. (1997) J. Mol. Graph. Mod. 15, 132-134 [CrossRef][Medline] [Order article via Infotrieve]
  43. Merritt, E. A., and Murphy, M. E. P. (1994) Acta Crystallogr. D 50, 869-873 [CrossRef][Medline] [Order article via Infotrieve]


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