Dimeric Crystal Structure of Rabbit l-Gulonate 3-Dehydrogenase/λ-Crystallin: Insights into the Catalytic Mechanism

l-Gulonate 3-dehydrogenase (GDH) is a bifunctional dimeric protein that functions not only as an NAD⁺-dependent enzyme in the uronate cycle but also as a taxon-specific λ-crystallin in rabbit lens. Here we report the first crystal structure of GDH in both apo form and NADH-bound holo form. The GDH protomer consists of two structural domains: the N-terminal domain with a Rossmann fold and the C-terminal domain with a novel helical fold. In the N-terminal domain of the NADH-bound structure, we identified 11 coenzyme-binding residues and found 2 distinct side-chain conformers of Ser124, which is a putative coenzyme/substrate-binding residue. A structural comparison between apo form and holo form and a mutagenesis study with E97Q mutant suggest an induced-fit mechanism upon coenzyme binding; coenzyme binding induces a conformational change in the coenzyme-binding residues Glu97 and Ser124 to switch their activation state from resting to active, which is required for the subsequent substrate recruitment. Subunit dimerization is mediated by numerous intersubunit interactions, including 22 hydrogen bonds and 104 residue pairs of van der Waals interactions, of which those between two cognate C-terminal domains are predominant. From a structure/sequence comparison within GDH homologues, a much greater degree of interprotomer interactions (both polar and hydrophobic) in the rabbit GDH would contribute to its higher thermostability, which may be relevant to the other function of this enzyme as λ-crystallin, a constitutive structural protein in rabbit lens. The present crystal structures and amino acid mutagenesis studies assigned the role of active-site residues: catalytic base for His145 and substrate binding for Ser124, Cys125, Asn196, and Arg231. Notably, Arg231 participates in substrate binding from the other subunit of the GDH dimer, indicating the functional significance of the dimeric state. Proper orientation of the substrate-binding residues for catalysis is likely to be maintained by an interprotomer hydrogen-bonding network of residues Asn196, Gln199, and Arg231, suggesting a network-based substrate recognition of GDH.     

The apo-holo hybrid of cytosolic aspartate aminotransferase, preparation and studies on subunit interactions

The apo/holo hydrid dimer of cytosolic aspartate aminotransferase from pig heart has been isolated by isoelectric focusing from a semi-reconstituted apo enzyme preparation. Comparison of the catalytic center activity and of the optical features of the coenzyme chromophore of the hybrid containing one active center with that of the homomer containing two revealed no differences in these properties and suggests that the subunits function independently from one another. However, differential thermal inactivation studies of hybrids and of homomers showed that the subunit carrying the coenzyme markedly stabilizes the structure of the neighboring subunit. Thus, coenzyme binding elicits functionally mute conformational changes which extend to and across the subunit interface.     

Smyd3 open & closed lock mechanism for substrate recruitment: The hinge motion of C-terminal domain inferred from μ-second molecular dynamics simulations

BackgroundThe human lysine methyltransferase Smyd3, a member of the SET and MYND domain containing protein family, harbors methylation activity on both histone and non-histone targets in a tightly regulated manner. The mechanism of how Smyd3 dynamically regulates substrate recognition is still not fully unveiled.MethodsHere, we employed molecular dynamics simulations on full length human Smyd3, performed to a total of 1.2 μ-second, in the presence (holo) and absence (apo) of the S-Adenosyl methionine (AdoMet) cofactor. The dynamical features of Smyd3 in apo and holo states have been examined and compared via examining geometrical and electrostatic properties.ResultsThe results show a distinct dynamics of the C-terminal domain (CTD) in the two states. In the apo state, the CTD undergoes a large hinge like motion and samples more opened configurations, thus acting like a loosened clamp and resulting in expanded substrate binding crevice. In the holo state, the CTD exhibits a restricted motion while the overall structure remains compact, mimicking a closed clamp. This leads to a localized increase in the negative potential at the substrate binding cleft. Further, solvent accessibility of critical residues at the target lysine access channel, important for methylation activity, is increased.ConclusionsWe postulate that AdoMet cofactor acts like a key and locks Smyd3 in a closed conformation. In effect, the cofactor binding restricts the elasticity of the CTD, presenting a compact substrate binding cleft with high negative potential, which may have implications on substrate recruitment via long range electrostatics.General significanceThe deletion of the CTD from Smyd3 has been shown to abolish the basal histone methylation activity. Our study highlights the importance of the CTD elasticity in shaping the substrate binding site for recognition and supports the previously proposed role of the CTD in stabilizing the active site for methylation activity.     

Crystal structures of the Chromobacterium violaceumω-transaminase reveal major structural rearrangements upon binding of coenzyme PLP

The bacterial ω-transaminase from Chromobacterium violaceum (Cv-ωTA, EC2.6.1.18) catalyses industrially important transamination reactions by use of the coenzyme pyridoxal 5'-phosphate (PLP). Here, we present four crystal structures of Cv-ωTA: two in the apo form, one in the holo form and one in an intermediate state, at resolutions between 1.35 and 2.4 ?. The enzyme is a homodimer with a molecular mass of ～ 100 kDa. Each monomer has an active site at the dimeric interface that involves amino acid residues from both subunits. The apo-Cv-ωTA structure reveals unique 'relaxed' conformations of three critical loops involved in structuring the active site that have not previously been seen in a transaminase. Analysis of the four crystal structures reveals major structural rearrangements involving elements of the large and small domains of both monomers that reorganize the active site in the presence of PLP. The conformational change appears to be triggered by binding of the phosphate group of PLP. Furthermore, one of the apo structures shows a disordered 'roof?' over the PLP-binding site, whereas in the other apo form and the holo form the 'roof' is ordered. Comparison with other known transaminase crystal structures suggests that ordering of the 'roof' structure may be associated with substrate binding in Cv-ωTA and some other transaminases. DATABASE: The atomic coordinates and structure factors for the Chromobacterium violaceumω-transaminase crystal structures can be found in the RCSB Protein Data Bank (http://www.rcsb.org) under the accession codes 4A6U for the holoenzyme, 4A6R for the apo1 form, 4A6T for the apo2 form and 4A72 for the mixed form STRUCTURED DIGITAL ABSTRACT: ? -transaminases and -transaminases bind by dynamic light scattering (View interaction) ? -transaminase and -transaminase bind by x-ray crystallography (View interaction) ? -transaminase and -transaminase bind by x-ray crystallography (View interaction).     

Computational study of conformational changes in human 3-hydroxy-3-methylglutaryl coenzyme reductase induced by substrate binding

Abstract

The enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is mainly involved in the regulation of cholesterol biosynthesis. HMGR catalyses the reduction of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) to mevalonate at the expense of two NADPH molecules in a two-step reversible reaction. In the present study, we constructed a model of human HMGR (hHMGR) to explore the conformational changes of HMGR in complex with HMG-CoA and NADPH. In addition, we analysed the complete sequence of the Flap domain using molecular dynamics (MD) simulations and principal component analysis (PCA). The simulations revealed that the Flap domain plays an important role in catalytic site activation and substrate binding. The apo form of hHMGR remained in an open state, while a substrate-induced closure of the Flap domain was observed for holo hHMGR. Our study also demonstrated that the phosphorylation of Ser872 induces significant conformational changes in the Flap domain that lead to a complete closure of the active site, suggesting three principal conformations for the first stage of hHMGR catalysis. Our results were consistent with previous proposed models for the catalytic mechanism of hHMGR.

Communicated by Ramaswamy H. Sarma     

Structural and functional role of nickel ions in urease by molecular dynamics simulation

Nickel ions play several roles in the biological processes of microorganisms and plants. Urease has a nickel-containing active site and catalyzes the hydrolysis of urea to yield ammonia and carbamate. In the present study, the role of nickel ions is examined using molecular dynamics simulations of the holo and apo forms. Nonbonded models used for the nickel ions provide good reproduction of the active-site structure as indicated in the crystallized structure. The results confirm that urease has a rigid active site in either its holo or its apo form. A new conformation of the flap is observed in apo urease. The connection between the metal center and His^α323 is proposed to be responsible for maintaining the flap conformation. The binding free energy of acetohydroxamic acid and urease is estimated using the molecular mechanics–generalized Born/surface area method. The binding free energy is primarily driven by electrostatic interactions in the presence of nickel ions. Normal mode analysis is employed to characterize the movements of the flap in apo urease.     

Temperature-induced conformational change at the catalytic site of Sulfolobus solfataricus alcohol dehydrogenase highlighted by Asn249Tyr substitution. A hydrogen/deuterium exchange, kinetic, and fluorescence quenching study

A combination of hydrogen/deuterium exchange, fluorescence quenching, and kinetic studies was used to acquire experimental evidence for the crystallographically hypothesized increase in local flexibility which occurs in thermophilic NAD(+)-dependent Sulfolobus solfataricus alcohol dehydrogenase (SsADH) upon substitution Asn249Tyr. The substitution, located at the adenine-binding site, proved to decrease the affinity for both coenzyme and substrate, rendering the mutant enzyme 6-fold more active when compared to the wild-type enzyme [Esposito et al. (2003) FEBS Lett. 539, 14-18]. The amide H/D exchange data show that the wild-type and mutant enzymes have similar global flexibility at 22 and 60 degrees C. However, the temperature dependence of the Stern-Volmer constant determined by acrylamide quenching shows that the increase in temperature affects the local flexibility differently, since the K(SV) increment is significantly higher for the wild-type than for the mutant enzyme over the range 18-45 degrees C. Interestingly, the corresponding van't Hoff plot (log K(SV) vs 1/T) proves nonlinear for the apo and holo wild-type and apo mutant enzymes, with a break at approximately 45 degrees C in all three cases due to a conformational change affecting the tryptophan microenvironment experienced by the quencher molecules. The Arrhenius and van't Hoff plots derived from the k(cat) and K(M) thermodependence measured with cyclohexanol and NAD(+) at different temperatures display an abrupt change of slope at 45-50 degrees C. This proves more pronounced in the case of the mutant enzyme compared to the wild-type enzyme due to a conformational change in the structure rather than to an overlapping of two or more rate-limiting reaction steps with different temperature dependencies of their rate constants. Three-dimensional analysis indicates that the observed conformational change induced by temperature is associated with the flexible loops directly involved in the substrate and coenzyme binding.     

Molecular modeling of formate dehydrogenase: the formation of the Michaelis complex

The formation of the reactive enzyme–substrate complex of formate dehydrogenase has been investigated by molecular dynamics techniques accounting for different conformational states of the enzyme. Simulations revealed that the transport of substrate to the active site through the substrate channel proceeds in the open conformation of enzyme due to the crucial role of the Arg284 residue acting as a vehicle. However, formate binding in the active site of the open conformation leads to the formation of a nonproductive enzyme–substrate complex. The productive Michaelis complex is formed only in the closed enzyme conformation after the substrate and coenzyme have bound, when required rigidity of the binding site and reactive formate orientation due to interactions with Arg284, Asn146, Ile122, and His332 residues is attained. Then, the high occupancy (up to 75%) of the reactive substrate–coenzyme conformation is reached, which was demonstrated by hybrid quantum mechanics/molecular mechanics simulations using various semiempirical Hamiltonians.     

Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase

The enzyme adenylate kinase (ADK) features two substrate binding domains that undergo large-scale motions during catalysis. In the apo state, the enzyme preferentially adopts a globally open state with accessible binding sites. Binding of two substrate molecules (AMP + ATP or ADP + ADP) results in a closed domain conformation, allowing efficient phosphoryl-transfer catalysis. We employed molecular dynamics simulations to systematically investigate how the individual domain motions are modulated by the binding of substrates. Two-dimensional free-energy landscapes were calculated along the opening of the two flexible lid domains for apo and holo ADK as well as for all single natural substrates bound to one of the two binding sites of ADK. The simulations reveal a strong dependence of the conformational ensembles on type and binding position of the bound substrates and a nonsymmetric behavior of the lid domains. Altogether, the ensembles suggest that, upon initial substrate binding to the corresponding lid site, the opposing lid is maintained open and accessible for subsequent substrate binding. In contrast, ATP binding to the AMP-lid induces global domain closing, preventing further substrate binding to the ATP-lid site. This might constitute a mechanism by which the enzyme avoids the formation of a stable but enzymatically unproductive state.     

Effects of single nucleotide changes on the binding and activity of RNA aptamers to human papillomavirus 16 E7 oncoprotein

ATP:Cobalamin adenosyltransferases catalyze the transfer a 5′-deoxyadenosyl moiety from ATP to cob(I)alamin in the synthesis of the Co–C bond of coenzyme B₁₂. There are three types of adenosyltransferases, CobA, PduO and EutT. Among these adenosyltransferases, the PduO-type adenosyltransferases is the most widely distributed enzyme. Structural comparisons between apo BcPduO and BcPduO in complex with MgATP revealed that the N-terminal strands of both structures were ordered, which is in contrast with the most previously available PduO-type adenosyltransferase structures. Furthermore, unlike other reported structures, apo BcPduO was bound to additional dioxane molecules causing a side chain conformational change at the Tyr30 residue, which is an important residue that mediates hydrogen bonding with ATP molecules upon binding of cobalamin to the active site. This study provides more structural information into the role of active site residues on substrate binding.     

The accessibility of the active site and conformation states of the beta 2 subunit of tryptophan synthase studied by fluorescence quenching

The rate of quenching of the fluorescence of pyridoxal 5'-phosphate in the active site of the beta 2 subunit of tryptophan synthase from Escherichia coli was measured to estimate the accessibility of the coenzyme to the small molecules iodide and acrylamide. The alpha subunit and the substrate L-serine substantially reduced the quenching rate. For iodide, the order of decreasing quenching was: Schiff's base of N alpha-acetyl-lysine with pyridoxal 5'-phosphate greater than holo beta 2 subunit greater than holo alpha 2 beta 2 complex approximately equal to holo beta 2 subunit + L-serine greater than holo alpha 2 beta 2 complex + L-serine. The coenzyme in the beta 2 subunit is apparently freely accessible to both iodide and acrylamide (kappa approximately equal to 2 X 10(9) M-1 s-1), but the alpha subunit and L-serine decrease the rate by factors of 2-5. Quenching of the fluorescence of the single tryptophan residue of the beta 2 subunit revealed that the apo and holo forms exist in different states, whereas the alpha subunit stabilizes a third conformation. As the alpha subunit binds to the beta 2 subunit, the tryptophan residue, which is within 2.2 nm of the active site of the beta 2 subunit, probably rotates with respect to the plane of the ring of the coenzyme, such that fluorescence energy transfer from tryptophan to pyridoxal phosphate is greatly reduced. The alpha subunit strongly protects the active-site ligand indole propanol phosphate from quenching with acrylamide, consistent with the active site being deep in a cleft in the protein. Iodide induces dissociation of the holo alpha 2 beta 2 complex [E. W. Miles & M. Moriguchi (1977) J. Biol. Chem. 252, 6594-6599]. The effect of iodide on the fluorescence properties of holo alpha 2 beta 2 complex allows us to estimate an upper limit for the dissociation constant for the alpha 2 beta 2 complex of 10(-8) M, in the absence of iodide.     

Prediction of interaction sites from apo 3D structures when the holo conformation is different

The predictability of catalytic and binding sites from apo structures is addressed for proteins that undergo significant conformational change upon binding. Theoretical microscopic titration curves (THEMATICS), an electrostatics-based method for the prediction of functional sites, is performed on a test set of 24 proteins with both apo and holo structures available. For 23 of these 24 proteins (96%), THEMATICS predicts the correct catalytic or binding site for both the apo and holo forms. For only one of the 24 proteins, THEMATICS makes the correct prediction for the holo structure but fails for the apo structure. The metrics used by THEMATICS to identify functional residues generally are larger in absolute value for the functional residues in the holo forms compared to the corresponding residues in the apo forms. However, even in the apo forms, these identifying metrics are still statistically significantly larger for functional residues than for residues not involved in catalysis or binding. This indicates that some of the unusual electrostatic properties of functional residues are preserved in the apo conformation. Evidence is presented that certain residues immediately surrounding the active catalytic and binding residues impart functionally important chemical and electrostatic properties to the active residues. At least parts of these microenvironments exist in the unbound conformations, such that THEMATICS is able to distinguish the functional residues even in the apo structures.     

Molecular Basis of Calcium-Sensitizing and Desensitizing Mutations of the Human Cardiac Troponin C Regulatory Domain: A Multi-Scale Simulation Study

Troponin C (TnC) is implicated in the initiation of myocyte contraction via binding of cytosolic and subsequent recognition of the Troponin I switch peptide. Mutations of the cardiac TnC N-terminal regulatory domain have been shown to alter both calcium binding and myofilament force generation. We have performed molecular dynamics simulations of engineered TnC variants that increase or decrease sensitivity, in order to understand the structural basis of their impact on TnC function. We will use the distinction for mutants that are associated with increased affinity and for those mutants with reduced affinity. Our studies demonstrate that for GOF mutants V44Q and L48Q, the structure of the physiologically-active site II binding site in the -free (apo) state closely resembled the -bound (holo) state. In contrast, site II is very labile for LOF mutants E40A and V79Q in the apo form and bears little resemblance with the holo conformation. We hypothesize that these phenomena contribute to the increased association rate, , for the GOF mutants relative to LOF. Furthermore, we observe significant positive and negative positional correlations between helices in the GOF holo mutants that are not found in the LOF mutants. We anticipate these correlations may contribute either directly to affinity or indirectly through TnI association. Our observations based on the structure and dynamics of mutant TnC provide rationale for binding trends observed in GOF and LOF mutants and will guide the development of inotropic drugs that target TnC.     

Crystal structure of the apo form of a β‐transaminase from Mesorhizobium sp. strain LUK

Pyridoxal 5′‐phosphate (PLP)‐dependent β‐transaminases (βTAs) reversibly catalyze transamination reactions by recognizing amino groups linked to the β‐carbon atoms of their substrates. Although several βTA structures have been determined as holo forms containing PLP, little is known about the effect of PLP on the conversion of the apo structure to the holo structure. We determined the crystal structure of the apo form of a βTA from Mesorhizobium sp. strain LUK at 2.2 Å resolution to elucidate how PLP affects the βTA structure. The structure revealed three major disordered regions near the active site. Structural comparison with the holo form also showed that the disordered regions in the apo form are ordered and partially adopt secondary structures in the holo form. These findings suggest that PLP incorporation into the active site contributes to the structural stability of the active site architecture, thereby forming the complete active site. Our results provide novel structural insights into the role of PLP in terms of active site formation.     

Entrapment of carbon dioxide in the active site of carbonic anhydrase II

The visualization at near atomic resolution of transient substrates in the active site of enzymes is fundamental to fully understanding their mechanism of action. Here we show the application of using CO(2)-pressurized, cryo-cooled crystals to capture the first step of CO(2) hydration catalyzed by the zinc-metalloenzyme human carbonic anhydrase II, the binding of substrate CO(2), for both the holo and the apo (without zinc) enzyme to 1.1A resolution. Until now, the feasibility of such a study was thought to be technically too challenging because of the low solubility of CO(2) and the fast turnover to bicarbonate by the enzyme (Liang, J. Y., and Lipscomb, W. N. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 3675-3679). These structures provide insight into the long hypothesized binding of CO(2) in a hydrophobic pocket at the active site and demonstrate that the zinc does not play a critical role in the binding or orientation of CO(2). This method may also have a much broader implication for the study of other enzymes for which CO(2) is a substrate or product and for the capturing of transient substrates and revealing hydrophobic pockets in proteins.     

Molecular modeling of formate dehydrogenase: the formation of the Michaelis complex

The formation of the reactive enzyme-substrate complex of formate dehydrogenase has been investigated by molecular dynamics techniques accounting for different conformational states of the enzyme. Simulations revealed that the transport of substrate to the active site through the substrate channel proceeds in the open conformation of enzyme due to the crucial role of the Arg284 residue acting as a vehicle. However, formate binding in the active site of the open conformation leads to the formation of a nonproductive enzyme-substrate complex. The productive Michaelis complex is formed only in the closed enzyme conformation after the substrate and coenzyme have bound, when required rigidity of the binding site and reactive formate orientation due to interactions with Arg284, Asn146, Ile122, and His332 residues is attained. Then, the high occupancy (up to 75%) of the reactive substrate-coenzyme conformation is reached, which was demonstrated by hybrid quantum mechanics/molecular mechanics simulations using various semiempirical Hamiltonians.     

Structural insight into poplar glutaredoxin C1 with a bridging iron-sulfur cluster at the active site

Glutaredoxins are glutathione-dependent enzymes that function to reduce disulfide bonds in vivo. Interestingly, a recent discovery indicates that some glutaredoxins can also exist in another form, an iron-sulfur protein [Lillig, C. H., et al. (2005) Proc. Natl. Acad. Sci. U.S.A. 102, 8168-8173]. This provides a direct connection between glutaredoxins and iron-sulfur proteins, suggesting a possible new regulatory role of iron-sulfur clusters along with the new functional switch of glutaredoxins. Biochemical studies have indicated that poplar glutaredoxin C1 (Grx-C1) is also such a biform protein. The apo form (monomer) of Grx-C1 is a regular glutaredoxin, and the holo form (dimer) is an iron-sulfur protein with a bridging [2Fe-2S] cluster. Here, we report the structural characterizations of poplar Grx-C1 in both the apo and holo forms by NMR spectroscopy. The solution structure of the reduced apo Grx-C1, which is the first plant Grx structure, shows a typical Grx fold. When poplar Grx-C1 forms a dimer with an iron-sulfur cluster, each subunit of the holo form still retains the overall fold of the apo form. The bridging iron-sulfur cluster in holo Grx-C1 is coordinated near the active site. In addition to the iron-sulfur cluster linker, helix alpha3 of each subunit is probably involved in the direct contact between the two subunits. Moreover, two glutathione molecules are identified in the vicinity of the iron-sulfur cluster and very likely participate in cluster coordination. Taken together, we propose that the bridging [2Fe-2S] cluster is coordinated by the first cysteine at the glutaredoxin active site from each subunit of holo Grx-C1, along with two cysteines from two glutathione molecules. Our studies reveal that holo Grx-C1 has a novel structural and iron-sulfur cluster coordination pattern for an iron-sulfur protein.     

NAD-dependent formate dehydrogenase of methylotrophic bacteria Pseudomonas sp. 101. II. Enzyme conformation at 3.0 A resolution

Three heavy atom isomorphous derivatives were used for the X-ray analysis of the holo form of NAD-dependent bacterial formate dehydrogenase (ternary complex enzyme-NAD-azide) at 3.0 A resolution. The enzyme subunit contains a catalytic and a coenzyme binding domain, with the active centre and the coenzyme binding site in the cleft between the domains. The polypeptide chain's fold and the position of 393 C alpha-atoms were determined. The secondary structure of the formate dehydrogenase was resolved. The structure of the NAD-binding domain is shown to be similar to that of other NAD-dependent enzymes.     

Crystal structure of a ternary complex of the alcohol dehydrogenase from Sulfolobus solfataricus

The crystal structure of a ternary complex of the alcohol dehydrogenase from the archaeon Sulfolobus solfataricus (SsADH) has been determined at 2.3 A. The asymmetric unit contains a dimer with a NADH and a 2-ethoxyethanol molecule bound to each subunit. The comparison with the apo structure of the enzyme reveals that this medium chain ADH undergoes a substantial conformational change in the apo-holo transition, accompanied by loop movements at the domain interface. The extent of domain closure is similar to that observed for the classical horse liver ADH, although some differences are found which can be related to the different oligomeric states of the enzymes. Compared to its apo form, the SsADH ternary complex shows a change in the ligation state of the active site zinc ion which is no longer bound to Glu69, providing additional evidence of the dynamic role played by the conserved glutamate residue in ADHs. In addition, the structure presented here allows the identification of the substrate site and hence of the residues that are important in the binding of both the substrate and the coenzyme.     

Human glutathione-dependent formaldehyde dehydrogenase. Structures of apo,binary, and inhibitory ternary complexes

The human glutathione-dependent formaldehyde dehydrogenase is unique among the structurally studied members of the alcohol dehydrogenase family in that it follows a random bi bi kinetic mechanism. The structures of an apo form of the enzyme, a binary complex with substrate 12-hydroxydodecanoic acid, and a ternary complex with NAD+ and the inhibitor dodecanoic acid were determined at 2.0, 2.3, and 2.3 A resolution by X-ray crystallography using the anomalous diffraction signal of zinc. The structures of the enzyme and its binary complex with the primary alcohol substrate, 12-hydroxydodecanoic acid, and the previously reported binary complex with the coenzyme show that the binding of the first substrate (alcohol or coenzyme) causes only minor changes to the overall structure of the enzyme. This is consistent with the random mechanism of the enzyme where either of the substrates binds to the free enzyme. The catalytic-domain position in these structures is intermediate to the "closed" and "open" conformations observed in class I alcohol dehydrogenases. More importantly, two different tetrahedral coordination environments of the active site zinc are observed in these structures. In the apoenzyme, the active site zinc is coordinated to Cys44, His66 and Cys173, and a water molecule. In the inhibitor complex, the coordination environment involves Glu67 instead of the solvent water molecule. The coordination environment involving Glu67 as the fourth ligand likely represents an intermediate step during ligand exchange at the active site zinc. These observations provide new insight into metal-assisted catalysis and substrate binding in glutathione-dependent formaldehyde dehydrogenase.