Probing the Molecular Basis of Substrate Specificity, Stereospecificity, and Catalysis in the Class II Pyruvate Aldolase, BphI

Computational insights into the binding modes, keto-enol tautomerization and stereo-electronically controlled decarboxylation of oxaloacetate in the active site of macrophomate synthase

Oxaloacetic acid (OAA) is a β-ketocarboxylic acid, which plays an important role as an intermediate in some metabolic pathways, including the tricarboxylic acid cycle, gluconeogenesis and fatty acid biosynthesis. Animal studies have indicated that supplementing oxaloacetic acid shows an increase of lifespan and other substantial health benefits including mitochondrial DNA protection, and protection of retinal, neural and pancreatic tissues. Most of the chemical transformations of OAA in the metabolic pathways have been extensively studied; however, the understanding of decarboxylation of OAA at the atomic level is relatively lacking. Here, we carried out MD simulations and combined quantum mechanical/molecular mechanical (QM/MM) calculations as an example to systematically elucidate the binding modes, keto-enol tautomerization and decarboxylation of OAA in the active site of macrophomate synthase (MPS), which is a Mg(II)-dependent bifunctional enzyme that catalyzes both the decarboxylation of OAA and [4+2] cycloaddition of 2-pyrone with the decarboxylated intermediate of OAA (pyruvate enolate). On the basis of our calculations, it was found that the Mg2+-coordinated oxaloacetate may exist in enol forms and keto forms. The four keto forms can be transformed into each other by simply rotating the C2-C3 single bond, nevertheless, the keto-enol tautomerization strictly requires the assistance of pocket water molecules. In addition, the decarboxylation is stereo-electronically controlled, i.e., it is the relative orientation of the terminal carboxyl anion that determines the rate of decarboxylation. As such, the chemistry of oxaloacetate in the active site of MPS is complex. On one hand, the most stable binding mode (K-I) may undergo enol-keto tautomerization to isomerize to the enol form, which may further react with the second substrate; on the other hand, K-I may isomerize to another binding mode K-II to proceed decarboxylation to generate pyruvate enolate and CO2. Starting from K-I, the enol-keto tautomerization corresponds to a barrier of 16.2 kcal mol-1, whereas the decarboxylation is associated with an overall barrier of 19.7 kcal mol-1. These findings may provide useful information for understanding the chemistry of OAA and the catalysis of related enzymes, and they are basically in agreement with the available experimental kinetic data.

Catalytic and structural insights into a stereospecific and thermostable Class II aldolase HpaI from Acinetobacter baumannii

Aldolases catalyze the reversible reactions of aldol condensation and cleavage and have strong potential for the synthesis of chiral compounds, widely used in pharmaceuticals. Here, we investigated a new Class II metal aldolase from the p-hydroxyphenylacetate degradation pathway in Acinetobacter baumannii, 4-hydroxy-2-keto-heptane-1,7-dioate aldolase (AbHpaI), which has various properties suitable for biocatalysis, including stereoselectivity/stereospecificity, broad aldehyde utilization, thermostability, and solvent tolerance. Notably, the use of Zn²⁺ by AbHpaI as a native cofactor is distinct from other enzymes in this class. AbHpaI can also use other metal ion (M²⁺) cofactors, except Ca²⁺, for catalysis. We found that Zn²⁺ yielded the highest enzyme complex thermostability (T_m of 87 °C) and solvent tolerance. All AbHpaI•M²⁺ complexes demonstrated preferential cleavage of (4R)-2-keto-3-deoxy-D-galactonate ((4R)-KDGal) over (4S)-2-keto-3-deoxy-D-gluconate ((4S)-KDGlu), with AbHpaI•Zn²⁺ displaying the highest R/S stereoselectivity ratio (sixfold higher than other M²⁺ cofactors). For the aldol condensation reaction, AbHpaI•M²⁺ only specifically forms (4R)-KDGal and not (4S)-KDGlu and preferentially catalyzes condensation rather than cleavage by ∼40-fold. Based on 11 X-ray structures of AbHpaI complexed with M²⁺ and ligands at 1.85 to 2.0 Å resolution, the data clearly indicate that the M²⁺ cofactors form an octahedral geometry with Glu151 and Asp177, pyruvate, and water molecules. Moreover, Arg72 in the Zn²⁺-bound form governs the stereoselectivity/stereospecificity of AbHpaI. X-ray structures also show that Ca²⁺ binds at the trimer interface via interaction with Asp51. Hence, we conclude that AbHpaI•Zn²⁺ is distinctive from its homologues in substrate stereospecificity, preference for aldol formation over cleavage, and protein robustness, and is attractive for biocatalytic applications.

Nucleophile Promiscuity of Engineered Class II Pyruvate Aldolase YfaU from E. Coli

Pyruvate-dependent aldolases exhibit a stringent selectivity for pyruvate, limiting application of their synthetic potential, which is a drawback shared with other existing aldolases. Structure-guided rational protein engineering rendered a 2-keto-3-deoxy-l-rhamnonate aldolase variant, fused with a maltose-binding protein (MBP-YfaU W23V/L216A), capable of efficiently converting larger pyruvate analogues, for example, those with linear and branched aliphatic chains, in aldol addition reactions. Combination of these nucleophiles with N-Cbz-alaninal (Cbz=benzyloxycarbonyl) and N-Cbz-prolinal electrophiles gave access to chiral building blocks, for example, derivatives of (2S,3S,4R)-4-amino-3-hydroxy-2-methylpentanoic acid (68 %, d.r. 90:10) and the enantiomer of dolaproine (33 %, d.r. 94:6) as well as a collection of unprecedented α-amino acid derivatives of the proline and pyrrolizidine type. Conversions varied between 6-93 % and diastereomeric ratios from 50:50 to 95:5 depending on the nucleophilic and electrophilic components.

Building Bridges: Biocatalytic C-C-Bond Formation toward Multifunctional Products

Carbon-carbon bond formation is the key reaction for organic synthesis to construct the carbon framework of organic molecules. The review gives a selection of biocatalytic C-C-bond-forming reactions which have been investigated during the last 5 years and which have already been proven to be applicable for organic synthesis. In most cases, the reactions lead to products functionalized at the site of C-C-bond formation (e.g., α-hydroxy ketones, aminoalcohols, diols, 1,4-diketones, etc.) or allow to decorate aromatic and heteroaromatic molecules. Furthermore, examples for cyclization of (non)natural precursors leading to saturated carbocycles are given as well as the stereoselective cyclopropanation of olefins affording cyclopropanes. Although many tools are already available, recent research also makes it clear that nature provides an even broader set of enzymes to perform specific C-C coupling reactions. The possibilities are without limit; however, a big library of variants for different types of reactions is required to have the specific enzyme for a desired specific (stereoselective) reaction at hand.

Binding and channeling of alternative substrates in the enzyme DmpFG: a molecular dynamics study

DmpFG is a bifunctional enzyme comprised of an aldolase subunit, DmpG, and a dehydrogenase subunit, DmpF. The aldehyde intermediate produced by the aldolase is channeled directly through a buried molecular channel in the protein structure from the aldolase to the dehydrogenase active site. In this study, we have investigated the binding of a series of progressively larger substrates to the aldolase, DmpG, using molecular dynamics. All substrates investigated are easily accommodated within the active site, binding with free energy values comparable to the physiological substrate 4-hydroxy-2-ketovalerate. Subsequently, umbrella sampling was utilized to obtain free energy surfaces for the aldehyde intermediates (which would be generated from the aldolase reaction on each of these substrates) to move through the channel to the dehydrogenase DmpF. Small substrates were channeled with limited barriers in an energetically feasible process. We show that the barriers preventing bulky intermediates such as benzaldehyde from moving through the wild-type protein can be removed by selective mutation of channel-lining residues, demonstrating the potential for tailoring this enzyme to allow its use for the synthesis of specific chemical products. Furthermore, positions of transient escape routes in this flexible channel were determined.

Mechanistic and bioinformatic investigation of a conserved active site helix in α-isopropylmalate synthase from Mycobacterium tuberculosis, a member of the DRE-TIM metallolyase superfamily

The characterization of functionally diverse enzyme superfamilies provides the opportunity to identify evolutionarily conserved catalytic strategies, as well as amino acid substitutions responsible for the evolution of new functions or specificities. Isopropylmalate synthase (IPMS) belongs to the DRE-TIM metallolyase superfamily. Members of this superfamily share common active site elements, including a conserved active site helix and an HXH divalent metal binding motif, associated with stabilization of a common enolate anion intermediate. These common elements are overlaid by variations in active site architecture resulting in the evolution of a diverse set of reactions that include condensation, lyase/aldolase, and carboxyl transfer activities. Here, using IPMS, an integrated biochemical and bioinformatics approach has been utilized to investigate the catalytic role of residues on an active site helix that is conserved across the superfamily. The construction of a sequence similarity network for the DRE-TIM metallolyase superfamily allows for the biochemical results obtained with IPMS variants to be compared across superfamily members and within other condensation-catalyzing enzymes related to IPMS. A comparison of our results with previous biochemical data indicates an active site arginine residue (R80 in IPMS) is strictly required for activity across the superfamily, suggesting that it plays a key role in catalysis, most likely through enolate stabilization. In contrast, differential results obtained from substitution of the C-terminal residue of the helix (Q84 in IPMS) suggest that this residue plays a role in reaction specificity within the superfamily.

Engineering aldolases as biocatalysts

Aldolases are seen as an attractive route to the production of biologically important compounds due to their ability to form carbon-carbon bonds. However, for many industrial reactions there are no naturally occurring enzymes, and so many different engineering approaches have been used to address this problem. Engineering methods have been used to alter the stability, substrate specificity and stereospecificity of aldolases to produce excellent enzymes for biocatalytic processes. Recently greater understanding of the aldolase mechanism has allowed many successes with both rational engineering approaches and computational design of aldolases. Rational engineering approaches have produced desired enzymes quickly and efficiently while combination of computational design with laboratory methods has created enzymes with activity approaching that of natural enzymes.

Mechanism of the dehydrogenase reaction of DmpFG and analysis of inter-subunit channeling efficiency and thermodynamic parameters in the overall reaction

The bifunctional, microbial enzyme DmpFG is comprised of two subunits: the aldolase, DmpG, and the dehydrogenase, DmpF. DmpFG is of interest due to its ability to channel substrates between the two spatially distinct active sites. While the aldolase is well studied, significantly less is known about the dehydrogenase. Steady-state kinetic measurements of the reverse reaction of DmpF confirmed that the dehydrogenase uses a ping-pong mechanism, with substrate inhibition by acetyl CoA indicating that NAD(+)/NADH and CoA/acetyl CoA bind to the same site in DmpF. The Km of DmpF for exogenous acetaldehyde as a substrate was 23.7 mM, demonstrating the necessity for the channel to deliver acetaldehyde directly from the aldolase to the dehydrogenase active site. A channeling assay on the bifunctional enzyme gave an efficiency of 93% indicating that less than 10% of the toxic acetaldehyde leaks out of the channel into the bulk media, prior to reaching the dehydrogenase active site. The thermodynamic activation parameters of the reactions catalyzed by the aldolase, the dehydrogenase and the DmpFG complex were determined. The Gibb's free energy of activation for the dehydrogenase reaction was lower than that obtained for the full DmpFG reaction, in agreement with the high kcat obtained for the dehydrogenase reaction in isolation. Furthermore, although both the DmpF and DmpG reactions occur with small, favorable entropies of activation, the full DmpFG reaction occurs with a negative entropy of activation. This supports the concept of allosteric structural communication between the two enzymes to coordinate their activities.

Characterization of an aldolase-dehydrogenase complex from the cholesterol degradation pathway of Mycobacterium tuberculosis

HsaF and HsaG are an aldolase and dehydrogenase from the cholesterol degradation pathway of Mycobacterium tuberculosis. HsaF could be heterologously expressed and purified as a soluble dimer, but the enzyme was inactive in the absence of HsaG. HsaF catalyzes the aldol cleavage of 4-hydroxy-2-oxoacids to produce pyruvate and an aldehyde. The enzyme requires divalent metals for activity, with a preference for Mn(2+). The Km values for 4-hydroxy-2-oxoacids were about 20-fold lower than observed for the aldolase homologue, BphI from the polychlorinated biphenyl degradation pathway. Acetaldehyde and propionaldehyde were channeled directly to the dehydrogenase, HsaG, without export to the bulk solvent where they were transformed to acyl-CoA in an NAD(+) and coenzyme A dependent reaction. HsaG is able to utilize aldehydes up to five carbons in length as substrates, with similar catalytic efficiencies. The HsaF-HsaG complex was crystallized and its structure was determined to a resolution of 1.93 Å. Substitution of serine 41 in HsaG with isoleucine or aspartate resulted in about 35-fold increase in Km for CoA but only 4-fold increase in Km dephospho-CoA, suggesting that this residue interacts with the 3'-ribose phosphate of CoA. A second protein annotated as a 4-hydroxy-2-oxopentanoic acid aldolase in M. tuberculosis (MhpE, Rv3469c) was expressed and purified, but was found to lack aldolase activity. Instead this enzyme was found to possess oxaloacetate decarboxylase activity, consistent with the conservation (with the 4-hydroxy-2-oxoacid aldolases) of residues involved in pyruvate enolate stabilization.

Catalytic properties of a bacterial acylating acetaldehyde dehydrogenase: evidence for several active oligomeric states and coenzyme A activation upon binding

Until the last decade, two unrelated aldehyde dehydrogenase (ALDH) superfamilies, i.e. the phosphorylating and non-phosphorylating superfamilies, were known to catalyze the oxidation of aldehydes to activated or non-activated acids. However, a third one was discovered by the crystal structure of a bifunctional enzyme 4-hydroxy-2-ketovalerate aldolase/acylating acetaldehyde dehydrogenase (DmpFG) from Pseudomonas sp. strain CF600 (Manjasetty et al., Proc. Natl. Acad. Sci. USA 100 (2003) 6992-6997). Indeed, DmpF exhibits a non-phosphorylating CoA-dependent ALDH activity, but is structurally related to the phosphorylating superfamily. In this study, we undertook the characterization of the catalytic and structural properties of MhpEF from Escherichia coli, an ortholog of DmpFG in which MhpF converts acetaldehyde, produced by the cleavage of 4-hydroxy-2-ketovalerate by MhpE, into acetyl-CoA. The kinetic data obtained under steady-state and pre-steady-state conditions show that the aldehyde dehydrogenase, MhpF, is active as a monomer, a unique feature relative to the phosphorylating and non-phosphorylating ALDH superfamilies. Our results also reveal that the catalytic properties of MhpF are not dependent on its oligomeric state, supporting the hypothesis of a structurally and catalytically independent entity. Moreover, the transthioesterification is shown to be rate-limiting and, when compared with a chemical model, its catalytic efficiency is increased 10(4)-fold. Therefore, CoA binding to MhpF increases its reactivity and optimizes its positioning relative to the thioacylenzyme intermediate, thus enabling the formation of an efficient deacylation complex.

Crystal structure of reaction intermediates in pyruvate class II aldolase: substrate cleavage, enolate stabilization, and substrate specificity

Crystal structures of divalent metal-dependent pyruvate aldolase, HpaI, in complex with substrate and cleavage products were determined to 1.8-2.0 Å resolution. The enzyme·substrate complex with 4-hydroxy-2-ketoheptane-1,7-dioate indicates that water molecule W2 bound to the divalent metal ion initiates C3-C4 bond cleavage. The binding mode of the aldehyde donor delineated a solvent-filled capacious binding locus lined with predominantly hydrophobic residues. The absence of direct interactions with the aldehyde aliphatic carbons accounts for the broad specificity and lack of stereospecific control by the enzyme. Enzymatic complex structures formed with keto acceptors, pyruvate, and 2-ketobutyrate revealed bidentate interaction with the divalent metal ion by C1-carboxyl and C2-carbonyl oxygens and water molecule W4 that is within close contact of the C3 carbon. Arg(70) assumes a multivalent role through its guanidinium moiety interacting with all active site enzymatic species: C2 oxygen in substrate, pyruvate, and ketobutyrate; substrate C4 hydroxyl; aldehyde C1 oxygen; and W4. The multiple interactions made by Arg(70) stabilize the negatively charged C4 oxygen following proton abstraction, the aldehyde alignment in aldol condensation, and the pyruvate enolate upon aldol cleavage as well as support proton exchange at C3. This role is corroborated by loss of aldol cleavage ability and pyruvate C3 proton exchange activity and by a 730-fold increase in the dissociation constant toward the pyruvate enolate analog oxalate in the R70A mutant. Based on the crystal structures, a mechanism is proposed involving the two enzyme-bound water molecules, W2 and W4, in acid/base catalysis that facilitates reversible aldol cleavage. The same reaction mechanism promotes decarboxylation of oxaloacetate.

Rational approaches for engineering novel functionalities in carbon-carbon bond forming enzymes

Enzymes that catalyze carbon-carbon bond formation can be exploited as biocatalyst for synthetic organic chemistry. However, natural enzymes frequently do not possess the required properties or specificities to catalyze industrially useful transformations. This mini-review describes recent work using knowledge-guided site-specific mutagenesis of key active site residues to alter substrate specificity, stereospecificity and reaction specificity of these enzymes. In addition, examples of de novo designed enzymes that catalyze C-C bond reactions not found in nature will be discussed.

Substrate specificity, substrate channeling, and allostery in BphJ: an acylating aldehyde dehydrogenase associated with the pyruvate aldolase BphI

BphJ, a nonphosphorylating acylating aldehyde dehydrogenase, catalyzes the conversion of aldehydes to form acyl-coenzyme A in the presence of NAD(+) and coenzyme A (CoA). The enzyme is structurally related to the nonacylating aldehyde dehydrogenases, aspartate-β-semialdehyde dehydrogenase and phosphorylating glyceraldehyde-3-phosphate dehydrogenase. Cys-131 was identified as the catalytic thiol in BphJ, and pH profiles together with site-specific mutagenesis data demonstrated that the catalytic thiol is not activated by an aspartate residue, as previously proposed. In contrast to the wild-type enzyme that had similar specificities for two- or three-carbon aldehydes, an I195A variant was observed to have a 20-fold higher catalytic efficiency for butyraldehyde and pentaldehyde compared to the catalytic efficiency of the wild type toward its natural substrate, acetaldehyde. BphJ forms a heterotetrameric complex with the class II aldolase BphI that channels aldehydes produced in the aldol cleavage reaction to the dehydrogenase via a molecular tunnel. Replacement of Ile-171 and Ile-195 with bulkier amino acid residues resulted in no more than a 35% reduction in acetaldehyde channeling efficiency, showing that these residues are not critical in gating the exit of the channel. Likewise, the replacement of Asn-170 in BphJ with alanine and aspartate did not substantially alter aldehyde channeling efficiencies. Levels of activation of BphI by BphJ N170A, N170D, and I171A were reduced by ≥3-fold in the presence of NADH and ≥4.5-fold when BphJ was undergoing turnover, indicating that allosteric activation of the aldolase has been compromised in these variants. The results demonstrate that the dehydrogenase coordinates the catalytic activity of BphI through allostery rather than through aldehyde channeling.

Protein-protein interactions and substrate channeling in orthologous and chimeric aldolase-dehydrogenase complexes

Bacterial aldolase-dehydrogenase complexes catalyze the last steps in the meta cleavage pathway of aromatic hydrocarbon degradation. The aldolase (TTHB246) and dehydrogenase (TTHB247) from Thermus thermophilus were separately expressed and purified from recombinant Escherichia coli. The aldolase forms a dimer, while the dehydrogenase is a monomer; these enzymes can form a stable tetrameric complex in vitro, consisting of two aldolase and two dehydrogenase subunits. Upon complex formation, the K(m) value of 4-hydroxy-2-oxopentanoate, the substrate of TTHB246, is decreased 4-fold while the K(m) of acetaldehyde, the substrate of TTHB247, is increased 3-fold. The k(cat) values of each enzyme were reduced by ~2-fold when they were in a complex. The half-life of TTHB247 at 50 °C increased by ~4-fold when it was in a complex with TTHB246. The acetaldehyde product from TTHB246 could be efficiently channelled directly to TTHB247, but the channeling efficiency for the larger propionaldehyde was ~40% lower. A single A324G substitution in TTHB246 increased the channeling efficiency of propionaldehyde to a value comparable to that of acetaldehyde. Stable and catalytically competent chimeric complexes could be formed between the T. thermophilus enzymes and the orthologous aldolase (BphI) and dehydrogenase (BphJ) from the biphenyl degradation pathway of Burkholderia xenovorans LB400. However, channeling efficiencies for acetaldehyde in these chimeric complexes were ~10%. Structural and sequence analysis suggests that interacting residues in the interface of the aldolase-dehydrogenase complex are highly conserved among homologues, but coevolution of partner enzymes is required to fine-tune this interaction to allow for efficient substrate channeling.

Biological channeling of a reactive intermediate in the bifunctional enzyme DmpFG

It has been hypothesized that the bifunctional enzyme DmpFG channels its intermediate, acetaldehyde, from one active site to the next using a buried intermolecular channel identified in the crystal structure. This channel appears to switch between an open and a closed conformation depending on whether the coenzyme NAD(+) is present or absent. Here, we applied molecular dynamics and metadynamics to investigate channeling within DmpFG in both the presence and absence of NAD(+). We found that substrate channeling within this enzyme is energetically feasible in the presence of NAD(+) but was less likely in its absence. Tyr-291, a proposed control point at the channel's entry, does not appear to function as a molecular gate. Instead, it is thought to orientate the substrate 4-hydroxy-2-ketovalerate in DmpG before reaction occurs, and may function as a proton shuttle for the DmpG reaction. Three hydrophobic residues at the channel's exit appear to have an important role in controlling the entry of acetaldehyde into the DmpF active site.

Improving upon nature: active site remodeling produces highly efficient aldolase activity toward hydrophobic electrophilic substrates

The substrate specificity of enzymes is frequently narrow and constrained by multiple interactions, limiting the use of natural enzymes in biocatalytic applications. Aldolases have important synthetic applications, but the usefulness of these enzymes is hampered by their narrow reactivity profile with unnatural substrates. To explore the determinants of substrate selectivity and alter the specificity of Escherichia coli 2-keto-3-deoxy-6-phosphogluconate (KDPG) aldolase, we employed structure-based mutagenesis coupled with library screening of mutant enzymes localized to the bacterial periplasm. We identified two active site mutations (T161S and S184L) that work additively to enhance the substrate specificity of this aldolase to include catalysis of retro-aldol cleavage of (4S)-2-keto-4-hydroxy-4-(2'-pyridyl)butyrate (S-KHPB). These mutations improve the value of k(cat)/K(M)(S-KHPB) by >450-fold, resulting in a catalytic efficiency that is comparable to that of the wild-type enzyme with the natural substrate while retaining high stereoselectivity. Moreover, the value of k(cat)(S-KHPB) for this mutant enzyme, a parameter critical for biocatalytic applications, is 3-fold higher than the maximal value achieved by the natural aldolase with any substrate. This mutant also possesses high catalytic efficiency for the retro-aldol cleavage of the natural substrate, KDPG, and a >50-fold improved activity for cleavage of 2-keto-4-hydroxy-octonoate, a nonfunctionalized hydrophobic analogue. These data suggest a substrate binding mode that illuminates the origin of facial selectivity in aldol addition reactions catalyzed by KDPG and 2-keto-3-deoxy-6-phosphogalactonate aldolases. Furthermore, targeting mutations to the active site provides a marked improvement in substrate selectivity, demonstrating that structure-guided active site mutagenesis combined with selection techniques can efficiently identify proteins with characteristics that compare favorably to those of naturally occurring enzymes.