Analysis of the class I aldolase binding site architecture based on the crystal structure of 2-deoxyribose-5-phosphate aldolase at 0.99A resolution

Enantiodivergent Synthesis of Both (R)- and (S)-Heteroaryl Aldols by Rationally Engineered Aldolases

Aldolases, especially 2-deoxyribose-5-phosphate aldolase (DERA) enzymes, have been widely employed to access key chiral precursors for various active pharmaceutical ingredients (APIs). This has been enabled by expanding their substrate scope toward non-natural acceptors and donors via protein engineering. In this study, we endeavored to broaden the acceptor substrate scope of DERA from Geobacillus sp. (DERAGeo) toward the heteroaryl aldehydes through a rational protein engineering approach. We successfully performed iterative saturation mutagenesis of DERAGeo, resulting in two enantiocomplementary variants, viz., (R)-selective DERAGeo-S185G and (S)-selective DERAGeo-T12I/S185A, with enhanced catalytic efficiencies and enantioselectivities. Remarkably, the natural enantioselectivity of DERAGeo was reversed by a single mutation (S185G). The synthetic applicability of the mutants was demonstrated by conducting aldol reactions on a semipreparative scale, from which both (R)- and (S)-enantiomers of heteroaryl aldols were isolated with high yields (up to 99%) and excellent enantiopurities (up to 99:1 e.r.).

Characterization of a widespread sugar phosphate-processing bacterial microcompartment

Many prokaryotes form Bacterial Microcompartments (BMCs) that encapsulate segments of specialized metabolic pathways to enhance catalysis. The various functions of metabolosomes, catabolic BMCs, are dictated by the signature enzyme that processes initial substrates of the confined pathway. The components and native functions of several metabolosomes have been experimentally characterized; however one of the most prevalent across all bacteria has yet to be studied. Sugar Phosphate Utilizing (SPU) BMC loci encode enzymes predicted to be involved in sugar phosphate metabolism. The SPU genetic loci are found in organisms occupying habitats ranging from soils to hot springs, highlighting the ubiquity of the SPU BMC. We bioinformatically characterized seven SPU subtypes, all which contain an enzyme unique to SPU BMCs, a deoxyribose 5-phosphate aldolase (DERA). Here, we define the fundamental characteristics of SPU BMCs and have expressed, purified, and characterized a set of SPU core enzymes. These include a protein-protein complex formed between a SPU BMC DERA and a predicted ribose 5-phosphate isomerase. Further, we show that the SPU BMC DERA is catalytically active and propose that it acts as the universal signature enzyme for the SPU BMC, with implications for fundamental understanding and biotechnological applications of SPU BMCs.

Shotgun metagenomic insights into secondary metabolite biosynthetic gene clusters reveal taxonomic and functional profiles of microbiomes in natural farmland soil

Antibiotic resistance is a worldwide problem that imposes a devastating effect on developing countries and requires immediate interventions. Initially, most of the antibiotic drugs were identified by culturing soil microbes. However, this method is prone to discovering the same antibiotics repeatedly. The present study employed a shotgun metagenomics approach to investigate the taxonomic diversity, functional potential, and biosynthetic capacity of microbiomes from two natural agricultural farmlands located in Bekeka and Welmera Choke Kebelle in Ethiopia for the first time. Analysis of the small subunit rRNA revealed bacterial domain accounting for 83.33% and 87.24% in the two selected natural farmlands. Additionally, the analysis showed the dominance of Proteobacteria representing 27.27% and 28.79% followed by Actinobacteria making up 12.73% and 13.64% of the phyla composition. Furthermore, the analysis revealed the presence of unassigned bacteria in the studied samples. The metagenome functional analysis showed 176,961 and 104, 636 number of protein-coding sequences (pCDS) from the two samples found a match with 172,655 and 102, 275 numbers of InterPro entries, respectively. The Genome ontology annotation suggests the presence of 5517 and 3293 pCDS assigned to the "biosynthesis process". Numerous Kyoto Encyclopedia of Genes and Genomes modules (KEGG modules) involved in the biosynthesis of terpenoids and polyketides were identified. Furthermore, both known and novel Biosynthetic gene clusters, responsible for the production of secondary metabolites, such as polyketide synthases, non-ribosomal peptide synthetase, ribosomally synthesized and post-translationally modified peptides (Ripp), and Terpene, were discovered. Generally, from the results it can be concluded that the microbiomes in the selected sampling sites have a hidden functional potential for the biosynthesis of secondary metabolites. Overall, this study can serve as a strong preliminary step in the long journey of bringing new antibiotics to the market.

Complete Computational Reaction Mechanism for Foldamer-Catalyzed Aldol Condensation

Foldamers, small synthetic peptides made of α and β-amino acids, have been found to be efficient catalysts for carbon-carbon bond-forming aldol reactions; of particular interest is their ability to catalyze macrocycle ring closure reactions. These catalysts feature a pair of amine groups that are aligned by the helical conformation and act in concert. Kinetic measurements show that the rate of the reaction depends on the identity of the amine side chains present. However, such kinetic analyses and other characterization techniques (e.g. mass spectrometry) can provide only limited information regarding the overall mechanism and rate-determining step of foldamer catalysis. We use semi-empirical density functional tight binding quantum mechanics molecular mechanics metadynamics simulations to determine the free energy and barrier for all elementary steps involved in the ring closure aldol reactions. We have performed calculations for 44 elementary reaction steps to identify key trends regarding amine identity, and provide insight into the intermediates and rate-limiting step of the catalytic cycle. From our results and other known aldol catalysts, we propose foldamer mutants which simulations predict to be better catalysts.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Unlocking Asymmetric Michael Additions in an Archetypical Class I Aldolase by Directed Evolution

Class I aldolases catalyze asymmetric aldol addition reactions and have found extensive application in the biocatalytic synthesis of chiral β-hydroxy-carbonyl compounds. However, the usefulness of these powerful enzymes for application in other C-C bond-forming reactions remains thus far unexplored. The redesign of class I aldolases to expand their catalytic repertoire to include non-native carboligation reactions therefore continues to be a major challenge. Here, we report the successful redesign of 2-deoxy-d-ribose-5-phosphate aldolase (DERA) from Escherichia coli, an archetypical class I aldolase, to proficiently catalyze enantioselective Michael additions of nitromethane to α,β-unsaturated aldehydes to yield various pharmaceutically relevant chiral synthons. After 11 rounds of directed evolution, the redesigned DERA enzyme (DERA-MA) carried 12 amino-acid substitutions and had an impressive 190-fold enhancement in catalytic activity compared to the wildtype enzyme. The high catalytic efficiency of DERA-MA for this abiological reaction makes it a proficient "Michaelase" with potential for biocatalytic application. Crystallographic analysis provides a structural context for the evolved activity. Whereas an aldolase acts naturally by activating the enzyme-bound substrate as a nucleophile (enamine-based mechanism), DERA-MA instead acts by activating the enzyme-bound substrate as an electrophile (iminium-based mechanism). This work demonstrates the power of directed evolution to expand the reaction scope of natural aldolases to include asymmetric Michael addition reactions and presents opportunities to explore iminium catalysis with DERA-derived catalysts inspired by developments in the organocatalysis field.

Current state of and need for enzyme engineering of 2-deoxy-D-ribose 5-phosphate aldolases and its impact

Abstract

Deoxyribose-5-phosphate aldolases (DERAs, EC 4.1.2.4) are acetaldehyde-dependent, Class I aldolases catalyzing in nature a reversible aldol reaction between an acetaldehyde donor (C2 compound) and glyceraldehyde-3-phosphate acceptor (C3 compound, C3P) to generate deoxyribose-5-phosphate (C5 compound, DR5P). DERA enzymes have been found to accept also other types of aldehydes as their donor, and in particular as acceptor molecules. Consequently, DERA enzymes can be applied in C–C bond formation reactions to produce novel compounds, thus offering a versatile biocatalytic alternative for synthesis. DERA enzymes, found in all kingdoms of life, share a common TIM barrel fold despite the low overall sequence identity. The catalytic mechanism is well-studied and involves formation of a covalent enzyme-substrate intermediate. A number of protein engineering studies to optimize substrate specificity, enzyme efficiency, and stability of DERA aldolases have been published. These have employed various engineering strategies including structure-based design, directed evolution, and recently also machine learning–guided protein engineering. For application purposes, enzyme immobilization and usage of whole cell catalysis are preferred methods as they improve the overall performance of the biocatalytic processes, including often also the stability of the enzyme. Besides single-step enzymatic reactions, DERA aldolases have also been applied in multi-enzyme cascade reactions both in vitro and in vivo. The DERA-based applications range from synthesis of commodity chemicals and flavours to more complicated and high-value pharmaceutical compounds.

Key points

• DERA aldolases are versatile biocatalysts able to make new C–C bonds.

• Synthetic utility of DERAs has been improved by protein engineering approaches.

• Computational methods are expected to speed up the future DERA engineering efforts.

Graphical abstract

Reverse C-glycosidase reaction provides C-nucleotide building blocks of xenobiotic nucleic acids

C-Analogues of the canonical N-nucleosides have considerable importance in medicinal chemistry and are promising building blocks of xenobiotic nucleic acids (XNA) in synthetic biology. Although well established for synthesis of N-nucleosides, biocatalytic methods are lacking in C-nucleoside synthetic chemistry. Here, we identify pseudouridine monophosphate C-glycosidase for selective 5-β-C-glycosylation of uracil and derivatives thereof from pentose 5-phosphate (D-ribose, 2-deoxy-D-ribose, D-arabinose, D-xylose) substrates. Substrate requirements of the enzymatic reaction are consistent with a Mannich-like addition between the pyrimidine nucleobase and the iminium intermediate of enzyme (Lys166) and open-chain pentose 5-phosphate. β-Elimination of the lysine and stereoselective ring closure give the product. We demonstrate phosphorylation-glycosylation cascade reactions for efficient, one-pot synthesis of C-nucleoside phosphates (yield: 33 - 94%) from unprotected sugar and nucleobase. We show incorporation of the enzymatically synthesized C-nucleotide triphosphates into nucleic acids by RNA polymerase. Collectively, these findings implement biocatalytic methodology for C-nucleotide synthesis which can facilitate XNA engineering for synthetic biology applications.

Substrate specificity of 2-deoxy-D-ribose 5-phosphate aldolase (DERA) assessed by different protein engineering and machine learning methods

In this work, deoxyribose-5-phosphate aldolase (Ec DERA, EC 4.1.2.4) from Escherichia coli was chosen as the protein engineering target for improving the substrate preference towards smaller, non-phosphorylated aldehyde donor substrates, in particular towards acetaldehyde. The initial broad set of mutations was directed to 24 amino acid positions in the active site or in the close vicinity, based on the 3D complex structure of the E. coli DERA wild-type aldolase. The specific activity of the DERA variants containing one to three amino acid mutations was characterised using three different substrates. A novel machine learning (ML) model utilising Gaussian processes and feature learning was applied for the 3rd mutagenesis round to predict new beneficial mutant combinations. This led to the most clear-cut (two- to threefold) improvement in acetaldehyde (C2) addition capability with the concomitant abolishment of the activity towards the natural donor molecule glyceraldehyde-3-phosphate (C3P) as well as the non-phosphorylated equivalent (C3). The Ec DERA variants were also tested on aldol reaction utilising formaldehyde (C1) as the donor. Ec DERA wild-type was shown to be able to carry out this reaction, and furthermore, some of the improved variants on acetaldehyde addition reaction turned out to have also improved activity on formaldehyde. KEY POINTS: • DERA aldolases are promiscuous enzymes. • Synthetic utility of DERA aldolase was improved by protein engineering approaches. • Machine learning methods aid the protein engineering of DERA.

Catalytic Foldamers: When the Structure Guides the Function

Enzymes are predominantly proteins able to effectively and selectively catalyze highly complex biochemical reactions in mild reaction conditions. Nevertheless, they are limited to the arsenal of reactions that have emerged during natural evolution in compliance with their intrinsic nature, three-dimensional structures and dynamics. They optimally work in physiological conditions for a limited range of reactions, and thus exhibit a low tolerance for solvent and temperature conditions. The de novo design of synthetic highly stable enzymes able to catalyze a broad range of chemical reactions in variable conditions is a great challenge, which requires the development of programmable and finely tunable artificial tools. Interestingly, over the last two decades, chemists developed protein secondary structure mimics to achieve some desirable features of proteins, which are able to interfere with the biological processes. Such non-natural oligomers, so called foldamers, can adopt highly stable and predictable architectures and have extensively demonstrated their attractiveness for widespread applications in fields from biomedical to material science. Foldamer science was more recently considered to provide original solutions to the de novo design of artificial enzymes. This review covers recent developments related to peptidomimetic foldamers with catalytic properties and the principles that have guided their design.

Rational engineering of 2-deoxyribose-5-phosphate aldolases for the biosynthesis of (R)-1,3-butanediol

Carbon-carbon bond formation is one of the most important reactions in biocatalysis and organic chemistry. In nature, aldolases catalyze the reversible stereoselective aldol addition between two carbonyl compounds, making them attractive catalysts for the synthesis of various chemicals. In this work, we identified several 2-deoxyribose-5-phosphate aldolases (DERAs) having acetaldehyde condensation activity, which can be used for the biosynthesis of (R)-1,3-butanediol (1,3BDO) in combination with aldo-keto reductases (AKRs). Enzymatic screening of 20 purified DERAs revealed the presence of significant acetaldehyde condensation activity in 12 of the enzymes, with the highest activities in BH1352 from Bacillus halodurans, TM1559 from Thermotoga maritima, and DeoC from Escherichia coli The crystal structures of BH1352 and TM1559 at 1.40-2.50 Å resolution are the first full-length DERA structures revealing the presence of the C-terminal Tyr (Tyr²²⁴ in BH1352). The results from structure-based site-directed mutagenesis of BH1352 indicated a key role for the catalytic Lys¹⁵⁵ and other active-site residues in the 2-deoxyribose-5-phosphate cleavage and acetaldehyde condensation reactions. These experiments also revealed a 2.5-fold increase in acetaldehyde transformation to 1,3BDO (in combination with AKR) in the BH1352 F160Y and F160Y/M173I variants. The replacement of the WT BH1352 by the F160Y or F160Y/M173I variants in E. coli cells expressing the DERA + AKR pathway increased the production of 1,3BDO from glucose five and six times, respectively. Thus, our work provides detailed insights into the molecular mechanisms of substrate selectivity and activity of DERAs and identifies two DERA variants with enhanced activity for in vitro and in vivo 1,3BDO biosynthesis.

Silver-coated magnetic nanocomposites induce growth inhibition and protein changes in foodborne bacteria

Cytotoxicity concerns of nanoparticles on animal or human bodies have led to the design of iron oxide core nanocomposites, coated with elemental silver to allow their magnetic removal from bio-mixtures. Although the antimicrobial effect of silver is well-described, the effects of nanoparticles derived from silver on microorganisms remain unfolded. Here, we characterized a customized magnetic silver nanocomposite (Ag-MNP) and evaluated its effects on bacterial growth and protein changes. The Ag-MNP displayed both longitudinal and round shapes under High-Resolution Transmission Electron Microscopy imaging, while the Energy Dispersive X-ray Spectroscopy and X-ray diffraction analysis confirmed the presence of Ag, Fe₃O₄ (Magnetite) and FeO₂ (Goethite). Optical density, bioluminescence imaging, and Colony Forming Unit assessments revealed that the presence of Ag-MNP induced strong dose-dependent bacteria (Escherichia coli O157:H7, Salmonella enterica serovar Typhimurium and S. Anatum) growth inhibition. The TEM imaging showed penetration and infiltration of bacteria by Ag-MNP, leading to membrane degeneration and vacuole formation. The presence of Ag-MNP led to fifteen up-regulated and nine down-regulated proteins (P < 0.05) that are involved in cell membrane synthesis, inhibition of protein synthesis, interference with DNA synthesis, and energy metabolism inhibition. This study provides insights to develop alternative antimicrobials to treat foodborne pathogens with antibiotic resistance avoidance.

2-Deoxy-D-ribose-5-phosphate aldolase (DERA): applications and modifications

2-Deoxy-D-ribose-5-phosphate aldolase (DERA) is a class I aldolase that offers access to several building blocks for organic synthesis. It catalyzes the stereoselective C-C bond formation between acetaldehyde and numerous other aldehydes. However, the practical application of DERA as a biocatalyst is limited by its poor tolerance towards industrially relevant concentrations of aldehydes, in particular acetaldehyde. Therefore, the development of proper experimental conditions, including protein engineering and/or immobilization on appropriate supports, is required. The present review is aimed to provide a brief overview of DERA, its history, and progress made in understanding the functioning of the enzyme. Furthermore, the current understanding regarding aldehyde resistance of DERA and the various optimizations carried out to modify this property are discussed.

Conformational Sampling of the Intrinsically Disordered C-Terminal Tail of DERA Is Important for Enzyme Catalysis

2-Deoxyribose-5-phosphate aldolase (DERA) catalyzes the reversible conversion of acetaldehyde and glyceraldehyde-3-phosphate into deoxyribose-5-phosphate. DERA is used as a biocatalyst for the synthesis of drugs such as statins and is a promising pharmaceutical target due to its involvement in nucleotide catabolism. Despite previous biochemical studies suggesting the catalytic importance of the C-terminal tyrosine residue found in several bacterial DERAs, the structural and functional basis of its participation in catalysis remains elusive because the electron density for the last eight to nine residues (i.e., the C-terminal tail) is absent in all available crystal structures. Using a combination of NMR spectroscopy and molecular dynamics simulations, we conclusively show that the rarely studied C-terminal tail of E. coli DERA (ecDERA) is intrinsically disordered and exists in equilibrium between open and catalytically relevant closed states, where the C-terminal tyrosine (Y259) enters the active site. Nuclear Overhauser effect distance restraints, obtained due to the presence of a substantial closed state population, were used to derive the solution-state structure of the ecDERA closed state. Real-time NMR hydrogen/deuterium exchange experiments reveal that Y259 is required for efficiency of the proton abstraction step of the catalytic reaction. Phosphate titration experiments show that, in addition to the phosphate-binding residues located near the active site, as observed in the available crystal structures, ecDERA contains previously unknown auxiliary phosphate-binding residues on the C-terminal tail which could facilitate in orienting Y259 in an optimal position for catalysis. Thus, we present significant insights into the structural and mechanistic importance of the ecDERA C-terminal tail and illustrate the role of conformational sampling in enzyme catalysis.

1H, 13C, and 15N backbone and sidechain resonance assignments of a monomeric variant of E. coli deoxyribose-5-phosphate aldolase

Deoxyribose-5-phosphate aldolase (DERA) catalyses the reversible conversion of 2-deoxyribose-5-phosphate (dR5P) into glyceraldehyde-3-phosphate (G3P) and acetaldehyde. For industrial applications, this enzyme is used in organic synthesis for aldol reactions between acetaldehyde as a donor and a wide range of aldehydes as acceptors. Here, we present a near complete set of sequence-specific ¹H, ¹³C and ¹⁵N resonance assignments of a 28 kDa monomeric variant of the Escherichia coli DERA. These assignments provide the basis for ongoing structural and dynamic analysis of DERA substrate specificity.

Structural insight for substrate tolerance to 2-deoxyribose-5-phosphate aldolase from the pathogen Streptococcus suis

2-deoxyribose-5-phosphate aldolase (DERA) is a class I aldolase that catalyzes aldol condensation of two aldehydes in the active site, which is particularly germane in drug manufacture. Structural and biochemical studies have shown that the active site of DERA is typically loosely packed and displays broader substrate specificity despite sharing conserved folding architecture with other aldolases. The most distinctive structural feature of DERA compared to other aldolases is short and flexible C-terminal region. This region is also responsible for substrate recognition. Therefore, substrate tolerance may be related to the C-terminal structural features of DERA. Here, we determined the crystal structures of full length and C-terminal truncated DERA from Streptococcus suis (SsDERA). In common, both contained the typical (α/β)8 TIM-barrel fold of class I aldolases. Surprisingly, C-terminal truncation resulting in missing the last α9 and β8 secondary elements, allowed DERA to maintain activity comparable to the fulllength enzyme. Specifically, Arg186 and Ser205 residues at the C-terminus appeared mutually supplemental or less indispensible for substrate phosphate moiety recognition. Our results suggest that DERA might adopt a shorter C-terminal region than conventional aldolases during evolution pathway, resulting in a broader range of substrate tolerance through active site flexibility.

Benchmarking pKa Prediction Methods for Residues in Proteins

Methods for estimation of pKa values of residues in proteins were tested on a set of benchmark proteins with experimentally known pKa values. The benchmark set includes 80 different residues (20 each for Asp, Glu, Lys, and His), half of which consists of significantly variant cases (ΔpKa ≥ 1 pKa unit from the amino acid in solution). The method introduced by Case and co-workers [J. Am. Chem. Soc. 2004, 126, 4167-4180], referred to as the molecular dynamics/generalized-Born/thermodynamic integration (MD/GB/TI) technique, gives a root-mean-square deviation (rmsd) of 1.4 pKa units on the benchmark set. The use of explicit waters in the immediate region surrounding the residue was shown to generally reduce high errors for this method. Longer simulation time was also shown to increase the accuracy of this method. The empirical approach developed by Jensen and co-workers [Proteins 2005, 61, 704-721], PROPKA, also gives an overall rmsd of 1.4 pKa units and is more or less accurate based on residue type-the method does very well for Lys and Glu, but less so for Asp and His. Likewise, the absolute deviation is quite similar for the two methods-5.2 for PROPKA and 5.1 for MD/GB/TI. A comparison of these results with several prediction methods from the literature is presented. The error in pKa prediction is analyzed as a function of variation of the pKa from that in water and the solvent accessible surface area (SASA) of the residue. A case study of the catalytic lysine residue in 2-deoxyribose-5-phosphate aldolase (DERA) is also presented.

Linking coupled motions and entropic effects to the catalytic activity of 2-deoxyribose-5-phosphate aldolase (DERA)

Local mutations in the phosphate binding group of DERA alter global conformation dynamics, catalytic activities and reaction entropies.

DERA, 2-deoxyribose-5-phosphate aldolase, catalyzes the retro-aldol cleavage of 2-deoxy-ribose-5-phosphate (dR5P) into glyceraldehyde-3-phosphate (G3P) and acetaldehyde in a branch of the pentose phosphate pathway. In addition to the physiological reaction, DERA also catalyzes the reverse addition reaction and, hence, is an interesting candidate for bio-catalysis of carbo-ligation reactions, which are central to synthetic chemistry. An obstacle to overcome for this enzyme to become a truly useful biocatalyst, however, is to relax the very strict dependency of this enzyme on phosphorylated substrates. We have studied herein the role of the non-canonical phosphate-binding site of this enzyme, consisting of Ser²³⁸ and Ser²³⁹, by site-directed and site-saturation mutagenesis, coupled to kinetic analysis of mutants. In addition, we have performed molecular dynamics simulations on the wild-type and four mutant enzymes, to analyse how mutations at this phosphate-binding site may affect the protein structure and dynamics. Further examination of the S239P mutant revealed that this variant increases the enthalpy change at the transition state, relative to the wild-type enzyme, but concomitant loss in entropy causes an overall relative loss in the TS free energy change. This entropy loss, as measured by the temperature dependence of catalysed rates, was mirrored in both a drastic loss in dynamics of the enzyme, which contributes to phosphate binding, as well as an overall loss in anti-correlated motions distributed over the entire protein. Our combined data suggests that the degree of anticorrelated motions within the DERA structure is coupled to catalytic efficiency in the DERA-catalyzed retro-aldol cleavage reaction, and can be manipulated for engineering purposes.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

Structural and functional divergence of the aldolase fold in Toxoplasma gondii

Parasites of the phylum Apicomplexa are highly successful pathogens of humans and animals worldwide. As obligate intracellular parasites, they have significant energy requirements for invasion and gliding motility that are supplied by various metabolic pathways. Aldolases have emerged as key enzymes involved in these pathways, and all apicomplexans express one or both of fructose 1,6-bisphosphate (F16BP) aldolase and 2-deoxyribose 5-phosphate (dR5P) aldolase (DERA). Intriguingly, Toxoplasma gondii, a highly successful apicomplexan parasite, expresses F16BP aldolase (TgALD1), d5RP aldolase (TgDERA), and a divergent dR5P aldolase-like protein (TgDPA) exclusively in the latent bradyzoite stage. While the importance of TgALD1 in glycolysis is well established and TgDERA is also likely to be involved in parasite metabolism, the detailed function of TgDPA remains elusive. To gain mechanistic insight into the function of different T. gondii aldolases, we first determined the crystal structures of TgALD1 and TgDPA. Structural analysis revealed that both aldolases adopt a TIM barrel fold accessorized with divergent secondary structure elements. Structural comparison of TgALD1 and TgDPA with members of their respective enzyme families revealed that, while the active-site residues are conserved in TgALD1, key catalytic residues are absent in TgDPA. Consistent with this observation, biochemical assays showed that, while TgALD1 was active on F16BP, TgDPA was inactive on dR5P. Intriguingly, both aldolases are competent to bind polymerized actin in vitro. Altogether, structural and biochemical analyses of T. gondii aldolase and aldolase-like proteins reveal diverse functionalization of the classic TIM barrel aldolase fold.

Small molecule probes to quantify the functional fraction of a specific protein in a cell with minimal folding equilibrium shifts

Although much is known about protein folding in buffers, it remains unclear how the cellular protein homeostasis network functions as a system to partition client proteins between folded and functional, soluble and misfolded, and aggregated conformations. Herein, we develop small molecule folding probes that specifically react with the folded and functional fraction of the protein of interest, enabling fluorescence-based quantification of this fraction in cell lysate at a time point of interest. Importantly, these probes minimally perturb a protein's folding equilibria within cells during and after cell lysis, because sufficient cellular chaperone/chaperonin holdase activity is created by rapid ATP depletion during cell lysis. The folding probe strategy and the faithful quantification of a particular protein's functional fraction are exemplified with retroaldolase, a de novo designed enzyme, and transthyretin, a nonenzyme protein. Our findings challenge the often invoked assumption that the soluble fraction of a client protein is fully folded in the cell. Moreover, our results reveal that the partitioning of destabilized retroaldolase and transthyretin mutants between the aforementioned conformational states is strongly influenced by cytosolic proteostasis network perturbations. Overall, our results suggest that applying a chemical folding probe strategy to other client proteins offers opportunities to reveal how the proteostasis network functions as a system to regulate the folding and function of individual client proteins in vivo.

Exploration of alternate catalytic mechanisms and optimization strategies for retroaldolase design

Designed retroaldolases have utilized a nucleophilic lysine to promote carbon-carbon bond cleavage of β-hydroxy-ketones via a covalent Schiff base intermediate. Previous computational designs have incorporated a water molecule to facilitate formation and breakdown of the carbinolamine intermediate to give the Schiff base and to function as a general acid/base. Here we investigate an alternative active-site design in which the catalytic water molecule was replaced by the side chain of a glutamic acid. Five out of seven designs expressed solubly and exhibited catalytic efficiencies similar to previously designed retroaldolases for the conversion of 4-hydroxy-4-(6-methoxy-2-naphthyl)-2-butanone to 6-methoxy-2-naphthaldehyde and acetone. After one round of site-directed saturation mutagenesis, improved variants of the two best designs, RA114 and RA117, exhibited among the highest kcat (>10(-3)s(-1)) and kcat/KM (11-25M(-1)s(-1)) values observed for retroaldolase designs prior to comprehensive directed evolution. In both cases, the >10(5)-fold rate accelerations that were achieved are within 1-3 orders of magnitude of the rate enhancements reported for the best catalysts for related reactions, including catalytic antibodies (kcat/kuncat=10(6) to 10(8)) and an extensively evolved computational design (kcat/kuncat>10(7)). The catalytic sites, revealed by X-ray structures of optimized versions of the two active designs, are in close agreement with the design models except for the catalytic lysine in RA114. We further improved the variants by computational remodeling of the loops and yeast display selection for reactivity of the catalytic lysine with a diketone probe, obtaining an additional order of magnitude enhancement in activity with both approaches.

Aggregation and cooperative effects in the aldol reactions of lithium enolates

Density functional theory and Car-Parrinello molecular dynamics simulations have been carried out for model aldol reactions involving aggregates of lithium enolates derived from acetaldehyde and acetone. Formaldehyde and acetone have been used as electrophiles. It is found that the geometries of the enolate aggregates are in general determined by the most favorable arrangements of the point charges within the respective Lin On clusters. The reactivity of the enolates follows the sequence monomer≫dimer>tetramer. In lithium aggregates, the initially formed aldol adducts must rearrange to form more stable structures in which the enolate and alkoxide oxygen atoms are within the respective Lin On clusters. Positive cooperative effects, similar to allosteric effects found in several proteins, are found for the successive aldol reactions in aggregates. The corresponding transition structures show in general sofa geometries.

Protein identification in two phases of 1,3-propanediol production by proteomic analysis

Proteomic analysis by two-dimensional electrophoresis (2D)-mass spectrometry was used to identify differentially expressed proteins in the Clostridium sp. native strain (IBUN 158B) in two phases of the 1,3-propanediol (1,3-PD) production (lag phase and exponential growth phase). Intracellular protein fraction extraction conditions were standardised, as well as the 2D electrophoresis. Differences were found between both of the growth phases evaluated here. Thirty-two of the differentially expressed proteins were chosen to be identified by tandem mass spectrometry (MALDI TOF/TOF). The presence of four enzymes implicated in the 1,3-PD metabolic pathway was recorded: one from the reductive route (1,3-propanediol dehydrogenase) and three from the oxidative route (3-hydroxybutyryl-CoA dehydrogenase, NADPH-dependent butanol dehydrogenase and phosphate butyryl transferase). The following enzymes which have not been previously reported for Clostridium sp., were also identified: phosphoglycerate kinase, glucose 6-phosphate isomerase, deoxyribose phosphate aldolase, transketolase, cysteine synthetase, O-acetylhomoserine sulphhydrylase, glycyl-tRNA ligase, aspartate-β-semialdehyde dehydrogenase, inosine-5-monophosphate dehydrogenase, aconitate hydratase and the PrsA protein. The foregoing provides a novel contribution towards knowledge of the native strain for the purpose of designing genetic manipulation strategies to obtain strains with high production of 1,3-PD.

BIOLOGICAL SIGNIFICANCE

The article "Protein identification in two phases of 1,3-propanediol production by proteomic analysis" provides a novel contribution towards knowledge regarding the Colombian Clostridium sp. native strain (IBUN 158B) because this is a new approximation in comparative proteomics in two phases of the bacterial growth and 1,3-propanediol (1,3-PD) production conditions. The proteomic studies are very important to identify the enzymes that are expressed at different stages of production and therefore genes of interest in the genetic manipulation strategies; the results can be taken into account in future studies in metabolic engineering when optimising 1,3-PD production, in a cost-effective process having direct industrial applications.

A highly productive, whole-cell DERA chemoenzymatic process for production of key lactonized side-chain intermediates in statin synthesis

Employing DERA (2-deoxyribose-5-phosphate aldolase), we developed the first whole-cell biotransformation process for production of chiral lactol intermediates useful for synthesis of optically pure super-statins such as rosuvastatin and pitavastatin. Herein, we report the development of a fed-batch, high-density fermentation with Escherichia coli BL21 (DE3) overexpressing the native E. coli deoC gene. High activity of this biomass allows direct utilization of the fermentation broth as a whole-cell DERA biocatalyst. We further show a highly productive bioconversion processes with this biocatalyst for conversion of 2-substituted acetaldehydes to the corresponding lactols. The process is evaluated in detail for conversion of acetyloxy-acetaldehyde with the first insight into the dynamics of reaction intermediates, side products and enzyme activity, allowing optimization of the feeding strategy of the aldehyde substrates for improved productivities, yields and purities. The resulting process for production of ((2S,4R)-4,6-dihydroxytetrahydro-2H-pyran-2-yl)methyl acetate (acetyloxymethylene-lactol) has a volumetric productivity exceeding 40 g L⁻¹ h⁻¹ (up to 50 g L⁻¹ h⁻¹) with >80% yield and >80% chromatographic purity with titers reaching 100 g L⁻¹. Stereochemical selectivity of DERA allows excellent enantiomeric purities (ee >99.9%), which were demonstrated on downstream advanced intermediates. The presented process is highly cost effective and environmentally friendly. To our knowledge, this is the first asymmetric aldol condensation process achieved with whole-cell DERA catalysis and it simplifies and extends previously developed DERA-catalyzed approaches based on the isolated enzyme. Finally, applicability of the presented process is demonstrated by efficient preparation of a key lactol precursor, which fits directly into the lactone pathway to optically pure super-statins.

Characterization and application of a newly synthesized 2-deoxyribose-5-phosphate aldolase

A codon-optimized 2-deoxyribose-5-phosphate aldolase (DERA) gene was newly synthesized and expressed in Escherichia coli to investigate its biochemical properties and applications in synthesis of statin intermediates. The expressed DERA was purified and characterized using 2-deoxyribose-5-phosphate as the substrate. The specific activity of recombinant DERA was 1.8 U/mg. The optimum pH and temperature for DERA activity were pH 7.0 and 35 °C, respectively. The recombinant DERA was stable at pH 4.0-7.0 and at temperatures below 50 °C. The enzyme activity was inhibited by 1 mM of Ni(2+), Ba(2+) and Fe(2+). The apparent K (m) and V (max) values of purified enzyme for 2-deoxyribose-5-phosphate were 0.038 mM and 2.9 μmol min(-1) mg(-1), for 2-deoxyribose were 0.033 mM and 2.59 μmol min(-1) mg(-1), respectively, which revealed that the enzyme had similar catalytic efficiency towards phosphorylated and non-phosphorylated substrates. To synthesize statin intermediates, the bioconversion process for production of (3R, 5S)-6-chloro-2,4,6-trideoxyhexose from chloroacetaldehyde and acetaldehyde by the recombinant DERA was developed and a conversion of 94.4 % was achieved. This recombinant DERA could be a potential candidate for application in production of (3R, 5S)-6-chloro-2,4,6-trideoxyhexose.

Using catalytic atom maps to predict the catalytic functions present in enzyme active sites

Catalytic atom maps (CAMs) are minimal models of enzyme active sites. The structures in the Protein Data Bank (PDB) were examined to determine if proteins with CAM-like geometries in their active sites all share the same catalytic function. We combined the CAM-based search protocol with a filter based on the weighted contact number (WCN) of the catalytic residues, a measure of the "crowdedness" of the microenvironment around a protein residue. Using this technique, a CAM based on the Ser-His-Asp catalytic triad of trypsin was able to correctly identify catalytic triads in other enzymes within 0.5 Å rmsd of the CAM with 96% accuracy. A CAM based on the Cys-Arg-(Asp/Glu) active site residues from the tyrosine phosphatase active site achieved 89% accuracy in identifying this type of catalytic functionality. Both of these CAMs were able to identify active sites across different fold types. Finally, the PDB was searched to locate proteins with catalytic functionality similar to that present in the active site of orotidine 5'-monophosphate decarboxylase (ODCase), whose mechanism is not known with certainty. A CAM, based on the conserved Lys-Asp-Lys-Asp tetrad in the ODCase active site, was used to search the PDB for enzymes with similar active sites. The ODCase active site has a geometry similar to that of Schiff base-forming Class I aldolases, with lowest aldolase rmsd to the ODCase CAM at 0.48 Å. The similarity between this CAM and the aldolase active site suggests that ODCase has the correct catalytic functionality present in its active site for the generation of a nucleophilic lysine.

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.

Comparison of two metal-dependent pyruvate aldolases related by convergent evolution: substrate specificity, kinetic mechanism, and substrate channeling

HpaI and BphI are two pyruvate class II aldolases found in aromatic meta-cleavage degradation pathways that catalyze similar reactions but are not related in sequence. Steady-state kinetic analysis of the aldol addition reactions and product inhibition assays showed that HpaI exhibits a rapid equilibrium random order mechanism while BphI exhibits a compulsory order mechanism, with pyruvate binding first. Both aldolases are able to utilize aldehyde acceptors two to five carbons in length; however, HpaI showed broader specificity and had a preference for aldehydes containing longer linear alkyl chains or C2-OH substitutions. Both enzymes were able to bind 2-keto acids larger than pyruvate, but only HpaI was able to utilize both pyruvate and 2-ketobutanoate as carbonyl donors in the aldol addition reaction. HpaI lacks stereospecific control producing racemic mixtures of 4-hydroxy-2-oxopentanoate (HOPA) from pyruvate and acetaldehyde while BphI synthesizes only (4S)-HOPA. BphI is also able to utilize acetaldehyde produced by the reduction of acetyl-CoA catalyzed by the associated aldehyde dehydrogenase, BphJ. This aldehyde was directly channeled from the dehydrogenase to the aldolase active sites, with an efficiency of 84%. Furthermore, the BphJ reductive deacylation reaction increased 4-fold when BphI was catalyzing the aldol addition reaction. Therefore, the BphI-BphJ enzyme complex exhibits unique bidirectionality in substrate channeling and allosteric activation.

Cloning, expression, and characterization of a new deoxyribose 5-phosphate aldolase from Yersinia sp. EA015

A new deoC gene encoding deoxyribose 5-phosphate aldolase (DERA) was identified in Yersinia sp. EA015 isolated from soil. The DERA gene had an open reading frame (ORF) of 672 base pairs encoding 223 amino acids to yield a protein of molecular mass 24.8 kDa. The amino acid sequence was 94% identical to that of DERA from Yersinia intermedia ATCC 29909. DERA was over-expressed in Escherichia coli and purified using Ni-NTA affinity chromatography. The specific activity was 137 micromol/min/mg. The Michaelis constant (k(m) value) of DERA was 9.1 mM. DERA was optimally active at pH 6.0 and 50 degrees C. DERA was tolerant to a high concentration (300 mM) of acetaldehyde.

Catalytic cycle of human glutathione reductase near 1 A resolution

Efficient enzyme catalysis depends on exquisite details of structure beyond those resolvable in typical medium- and high-resolution crystallographic analyses. Here we report synchrotron-based cryocrystallographic studies of natural substrate complexes of the flavoenzyme human glutathione reductase (GR) at nominal resolutions between 1.1 and 0.95 A that reveal new aspects of its mechanism. Compression in the active site causes overlapping van der Waals radii and distortion in the nicotinamide ring of the NADPH substrate, which enhances catalysis via stereoelectronic effects. The bound NADPH and redox-active disulfide are positioned optimally on opposite sides of the flavin for a 1,2-addition across a flavin double bond. The new structures extend earlier observations to reveal that the redox-active disulfide loop in GR is an extreme case of sequential peptide bonds systematically deviating from planarity--a net deviation of 53 degrees across five residues. But this apparent strain is not a factor in catalysis, as it is present in both oxidized and reduced structures. Intriguingly, the flavin bond lengths in oxidized GR are intermediate between those expected for oxidized and reduced flavin, but we present evidence that this may not be due to the protein environment but instead due to partial synchrotron reduction of the flavin by the synchrotron beam. Finally, of more general relevance, we present evidence that the structures of synchrotron-reduced disulfide bonds cannot generally be used as reliable models for naturally reduced disulfide bonds.

Lysine acetylation is a highly abundant and evolutionarily conserved modification in Escherichia coli

Lysine acetylation and its regulatory enzymes are known to have pivotal roles in mammalian cellular physiology. However, the extent and function of this modification in prokaryotic cells remain largely unexplored, thereby presenting a hurdle to further functional study of this modification in prokaryotic systems. Here we report the first global screening of lysine acetylation, identifying 138 modification sites in 91 proteins from Escherichia coli. None of the proteins has been previously associated with this modification. Among the identified proteins are transcriptional regulators, as well as others with diverse functions. Interestingly, more than 70% of the acetylated proteins are metabolic enzymes and translation regulators, suggesting an intimate link of this modification to energy metabolism. The new dataset suggests that lysine acetylation could be abundant in prokaryotic cells. In addition, these results also imply that functions of lysine acetylation beyond regulation of gene expression are evolutionarily conserved from bacteria to mammals. Furthermore, we demonstrate that bacterial lysine acetylation is regulated in response to stress stimuli.

Characterization and crystal structure of Escherichia coli KDPGal aldolase

2-Keto-3-deoxy-6-phosphogluconate (KDPG) and 2-keto-3-deoxy-6-phosphogalactonate (KDPGal) aldolases catalyze an identical reaction differing in substrate specificity in only the configuration of a single stereocenter. However, the proteins show little sequence homology at the amino acid level. Here we investigate the determinants of substrate selectivity of these enzymes. The Escherichia coli KDPGal aldolase gene, cloned into a T7 expression vector and overexpressed in E. coli, catalyzes retro-aldol cleavage of the natural substrate, KDPGal, with values of k(cat)/K(M) and k(cat) of 1.9x10(4)M(-1)s(-1) and 4s(-1), respectively. In the synthetic direction, KDPGal aldolase efficiently catalyzes an aldol addition using a limited number of aldehyde substrates, including d-glyceraldehyde-3-phosphate (natural substrate), d-glyceraldehyde, glycolaldehyde, and 2-pyridinecarboxaldehyde. A preparative scale reaction between 2-pyridinecarboxaldehyde and pyruvate catalyzed by KDPGal aldolase produced the aldol adduct of the R stereochemistry in >99.7% ee, a result complementary to that observed using the related KDPG aldolase. The native crystal structure has been solved to a resolution of 2.4A and displays the same (alpha/beta)(8) topology, as KDPG aldolase. We have also determined a 2.1A structure of a Schiff base complex between the enzyme and its substrate. This model predicts that a single amino acid change, T161 in KDPG aldolase to V154 in KDPGal aldolase, plays an important role in determining the stereochemical course of enzyme catalysis and this prediction was borne out by site-directed mutagenesis studies. However, additional changes in the enzyme sequence are required to prepare an enzyme with both high catalytic efficiency and altered stereochemistry.

Carbonylation of epoxides to substituted 3-hydroxy-delta-lactones

Substituted 3-hydroxy-delta-lactones (3HLs) are valuable intermediates in the synthesis of pharmaceuticals and other biologically active natural products. Herein we report the first example of the catalytic carbonylation of substituted homoglycidols to 3HLs using HCo(CO)4. Upon optimization of the catalyst and reaction conditions, a functionally diverse set of 3HLs was prepared. Mechanistic insight was gained by observation of the carbonylation reaction using in situ IR spectroscopy, and we propose a mechanism that is consistent with previously studied epoxide carbonylation systems.

Sequential aldol condensation catalyzed by hyperthermophilic 2-deoxy-d-ribose-5-phosphate aldolase

Genes encoding 2-deoxy-d-ribose-5-phosphate aldolase (DERA) homologues from two hyperthermophiles, the archaeon Pyrobaculum aerophilum and the bacterium Thermotoga maritima, were expressed individually in Escherichia coli, after which the structures and activities of the enzymes produced were characterized and compared with those of E. coli DERA. To our surprise, the two hyperthermophilic DERAs showed much greater catalysis of sequential aldol condensation using three acetaldehydes as substrates than the E. coli enzyme, even at a low temperature (25 degrees C), although both enzymes showed much less 2-deoxy-d-ribose-5-phosphate synthetic activity. Both the enzymes were highly resistant to high concentrations of acetaldehyde and retained about 50% of their initial activities after a 20-h exposure to 300 mM acetaldehyde at 25 degrees C, whereas the E. coli DERA was almost completely inactivated after a 2-h exposure under the same conditions. The structure of the P. aerophilum DERA was determined by X-ray crystallography to a resolution of 2.0 A. The main chain coordinate of the P. aerophilum enzyme monomer was quite similar to those of the T. maritima and E. coli enzymes, whose crystal structures have already been solved. However, the quaternary structure of the hyperthermophilic enzymes was totally different from that of the E. coli DERA. The areas of the subunit-subunit interface in the dimer of the hyperthermophilic enzymes are much larger than that of the E. coli enzyme. This promotes the formation of the unique dimeric structure and strengthens the hydrophobic intersubunit interactions. These structural features are considered responsible for the extremely high stability of the hyperthermophilic DERAs.

Directed evolution of an industrial biocatalyst: 2-deoxy-D-ribose 5-phosphate aldolase

Aldolases are emerging as powerful and cost efficient tools for the industrial synthesis of chiral molecules. They catalyze enantioselective carbon-carbon bond formations, generating up to two chiral centers under mild reaction conditions. Despite their versatility, narrow substrate ranges and enzyme inactivation under synthesis conditions represented major obstacles for large-scale applications of aldolases. In this study we applied directed evolution to optimize Escherichia coli 2-deoxy-D-ribose 5-phosphate aldolase (DERA) as biocatalyst for the industrial synthesis of (3R,5S)-6-chloro-2,4,6-trideoxyhexapyranoside. This versatile chiral precursor for vastatin drugs like Lipitor (atorvastatin) is synthesized by DERA in a tandem-aldol reaction from chloroacetaldehyde and two acetaldehyde equivalents. However, E. coli DERA shows low affinity to chloroacetaldehyde and is rapidly inactivated at aldehyde concentrations useful for biocatalysis. Using high-throughput screenings for chloroacetaldehyde resistance and for higher productivity, several improved variants have been identified. By combination of the most beneficial mutations we obtained a tenfold improved variant compared to wild-type DERA with regard to (3R,5S)-6-chloro-2,4,6-trideoxyhexapyranoside synthesis, under industrially relevant conditions.

Mechanism of the Class I KDPG aldolase

In vivo, 2-keto-3-deoxy-6-phosphogluconate (KDPG) aldolase catalyzes the reversible, stereospecific retro-aldol cleavage of KDPG to pyruvate and D-glyceraldehyde-3-phosphate. The enzyme is a lysine-dependent (Class I) aldolase that functions through the intermediacy of a Schiff base. Here, we propose a mechanism for this enzyme based on crystallographic studies of wild-type and mutant aldolases. The three dimensional structure of KDPG aldolase from the thermophile Thermotoga maritima was determined to 1.9A. The structure is the standard alpha/beta barrel observed for all Class I aldolases. At the active site Lys we observe clear density for a pyruvate Schiff base. Density for a sulfate ion bound in a conserved cluster of residues close to the Schiff base is also observed. We have also determined the structure of a mutant of Escherichia coli KDPG aldolase in which the proposed general acid/base catalyst has been removed (E45N). One subunit of the trimer contains density suggesting a trapped pyruvate carbinolamine intermediate. All three subunits contain a phosphate ion bound in a location effectively identical to that of the sulfate ion bound in the T. maritima enzyme. The sulfate and phosphate ions experimentally locate the putative phosphate binding site of the aldolase and, together with the position of the bound pyruvate, facilitate construction of a model for the full-length KDPG substrate complex. The model requires only minimal positional adjustments of the experimentally determined covalent intermediate and bound anion to accommodate full-length substrate. The model identifies the key catalytic residues of the protein and suggests important roles for two observable water molecules. The first water molecule remains bound to the enzyme during the entire catalytic cycle, shuttling protons between the catalytic glutamate and the substrate. The second water molecule arises from dehydration of the carbinolamine and serves as the nucleophilic water during hydrolysis of the enzyme-product Schiff base. The second water molecule may also mediate the base-catalyzed enolization required to form the carbon nucleophile, again bridging to the catalytic glutamate. Many aspects of this mechanism are observed in other Class I aldolases and suggest a mechanistically and, perhaps, evolutionarily related family of aldolases distinct from the N-acetylneuraminate lyase (NAL) family.