Reaction mechanism of N-acetylneuraminic acid lyase revealed by a combination of crystallography, QM/MM simulation, and mutagenesis

Biochemical and structural basis of sialic acid utilization by gut microbes

The human gastrointestinal (GI) tract harbors diverse microbial communities collectively known as the gut microbiota that exert a profound impact on human health and disease. The repartition and availability of sialic acid derivatives in the gut have a significant impact on the modulation of gut microbes and host susceptibility to infection and inflammation. Although N-acetylneuraminic acid (Neu5Ac) is the main form of sialic acids in humans, the sialic acid family regroups more than 50 structurally and chemically distinct modified derivatives. In the GI tract, sialic acids are found in the terminal location of mucin glycan chains constituting the mucus layer and also come from human milk oligosaccharides in the infant gut or from meat-based foods in adults. The repartition of sialic acid in the GI tract influences the gut microbiota composition and pathogen colonization. In this review, we provide an update on the mechanisms underpinning sialic acid utilization by gut microbes, focusing on sialidases, transporters, and metabolic enzymes.

One-pot biosynthesis of N-acetylneuraminic acid from chitin via combination of chitin-degrading enzymes, N-acetylglucosamine-2-epimerase, and N-neuraminic acid aldolase

N-acetylneuraminic acid (Neu5Ac) possesses the ability to promote mental health and enhance immunity and is widely used in both medicine and food fields as a supplement. Enzymatic production of Neu5Ac using N-acetyl-D-glucosamine (GlcNAc) as substrate was significant. However, the high-cost GlcNAc limited its development. In this study, an in vitro multi-enzyme catalysis was built to produce Neu5Ac using affordable chitin as substrate. Firstly, exochitinase SmChiA from Serratia proteamaculans and N-acetylglucosaminosidase CmNAGase from Chitinolyticbacter meiyuanensis SYBC-H1 were screened and combined to produce GlcNAc, effectively. Then, the chitinase was cascaded with N-acetylglucosamine-2-epimerase (AGE) and N-neuraminic acid aldolase (NanA) to produce Neu5Ac; the optimal conditions of the multi-enzyme catalysis system were 37°C and pH 8.5, the ratio of AGE to NanA (1:4) and addition of pyruvate (70 mM), respectively. Finally, 9.2 g/L Neu5Ac could be obtained from 20 g/L chitin within 24 h along with two supplementations with pyruvate. This work will lay a good foundation for the production of Neu5Ac from cheap chitin resources.

A mutagenic analysis of NahE, a hydratase-aldolase in the naphthalene degradative pathway

NahE is a hydratase-aldolase that converts o-substituted trans-benzylidenepyruvates (H, OH, or CO2-) to benzaldehyde, salicylaldehyde, or 2-carboxybenzaldehyde, respectively, and pyruvate. The enzyme is in a bacterial degradative pathway for naphthalene, which is a toxic and persistent environmental contaminant. Sequence, crystallographic, and mutagenic analysis identified the enzyme as a member of the N-acetylneuraminate lyase (NAL) subgroup in the aldolase superfamily. As such, it has a conserved lysine (Lys183) and tyrosine (Tyr155), for Schiff base formation, as well as a GXXGE motif for binding of the pyruvoyl carboxylate group. A crystal structure of the selenomethionine derivative of NahE shows these active site elements along with nearby residues that might be involved in the mechanism and/or specificity. Mutations of five active site amino acids (Thr65, Trp128, Tyr155, Asn157, and Asn281) were constructed and kinetic parameters measured in order to assess the effect(s) on catalysis. The results show that the two Trp128 mutants (Phe and Tyr) have the least effect on catalysis, whereas amino acids with bulky side chains at Thr65 (Val) and Asn281 (Leu) have the greatest effect. Changing Tyr155 to Phe and Asn157 to Ala also hinders catalysis, and the effects fall in between these extremes. These observations are put into a structural context using a crystal structure of the Schiff base of the reaction intermediate. Trapping experiments with substrate, Na(CN)BH3, and wild type enzyme and selected mutants mostly paralleled the kinetic analysis, and identified two salicylaldehyde-modified lysines: the active site lysine (Lys183) and one outside the active site (Lys279). The latter could be responsible for the observed inhibition of NahE by salicylaldehyde. Together, the results provide new insights into the NahE-catalyzed reaction.

Structurally Informed Mutagenesis of a Stereochemically Promiscuous Aldolase Produces Mutants That Catalyze the Diastereoselective Syntheses of All Four Stereoisomers of 3-Deoxy-hexulosonic Acid

A 2-keto-3-deoxygluconate aldolase from the hyperthermophile Sulfolobus solfataricus catalyzes the nonstereoselective aldol reaction of pyruvate and d-glyceraldehyde to produce 2-keto-3-deoxygluconate (d-KDGlc) and 2-keto-3-deoxy-d-galactonate (d-KDGal). Previous investigations into curing the stereochemical promiscuity of this hyperstable aldolase used high-resolution structures of the aldolase bound to d-KDGlc or d-KDGal to identify critical amino acids involved in substrate binding for mutation. This structure-guided approach enabled mutant variants to be created that could stereoselectively catalyze the aldol reaction of pyruvate and natural d-glyceraldehyde to selectively afford d-KDGlc or d-KDGal. Here we describe the creation of two further mutants of this Sulfolobus aldolase that can be used to catalyze aldol reactions between pyruvate and non-natural l-glyceraldehyde to enable the diastereoselective synthesis of l-KDGlc and l-KDGal. High-resolution crystal structures of all four variant aldolases have been determined (both unliganded and liganded), including Variant 1 with d-KDGlc, Variant 2 with pyruvate, Variant 3 with l-KDGlc, and Variant 4 with l-KDGal. These structures have enabled us to rationalize the observed changes in diastereoselectivities in these variant-catalyzed aldol reactions at a molecular level. Interestingly, the active site of Variant 4 was found to be sufficiently flexible to enable catalytically important amino acids to be replaced while still retaining sufficient enzymic activity to enable production of l-KDGal.

Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory

An improved understanding of the P450 structure is relevant to the development of biomimetic catalysts and inhibitors for controlled CH-bond activation, an outstanding challenge of synthetic chemistry. Motivated by the experimental findings of an unusually short Fe-S bond of 2.18 Å for the wild-type (WT) OleT P450 decarboxylase relative to a cysteine pocket mutant form (A369P), a computational model that captures the effect of the thiolate axial ligand on the iron-sulfur distance is presented. With the computational efficiency and streamlined analysis in mind, this model combines a cluster representation of the enzyme─40-110 atoms, depending on the heme and ligand truncation level─with a density functional theory (DFT) description of the electronic structure (ES) and is calibrated against the experimental data. The optimized Fe-S distances show a difference of 0.25 Å between the low and high spin states, in agreement with the crystallographic structures of the OleT WT and mutant forms. We speculate that this difference is attributable to the packing of the ligand; the mutant is bulkier due to an alanine-to-proline replacement, meaning that it is excluded from the energetically favored low-spin minimum because of steric constraints. The presence of pure spin-state pairs and the intersection of the low/high spin states for the enzyme model is indicative of the limitations of single-reference ES methods in such systems and emphasizes the significance of using the proper state when modeling the hydrogen atom transfer (HAT) reaction catalyzed by OleT. At the same time, the correct characterization of both the short and long Fe-S bonds within a small DFT-based model of 42 atoms paves the way for quantum dynamics modeling of the HAT step, which initiates the OleT decarboxylation reaction.

Recombinant snail sialic acid aldolase is promiscuous towards aliphatic aldehydes

Aldolases are enzymes that reversibly catalyze the cleavage of carbon-carbon bonds. Here we describe a recombinant sialic acid aldolase originating from the freshwater snail Biomphalaria glabrata (sNPL), and compare its substrate spectrum with a sialic acid aldolase originating from chicken (chNPL). In contrast to vertebrate animals which can synthesize, degrade, and incorporate sialic acids on glycoconjugate ubiquitously, snails (as all mollusks) cannot synthesize sialic acids endogenously, and therefore the biological function and substrate scope of sNPL ought to differ significantly from vertebrate sialic aldolases such as chNPL. sNPL was active towards a series of sialic acid derivatives but was in contrast to chNPL unable to catalyze the cleavage of N-acetylneuraminic acid into N-acetylmannosamine and pyruvate. Interestingly, chNPL and sNPL showed contrasting C4 (R) / (S) diastereoselectivity towards the substrates d-mannose and d-galactose in the presence of pyruvate. In addition, sNPL was also able to synthesize a series of 4-hydroxy-2-oxoates using the corresponding aliphatic aldehyde substrates in the presence of pyruvate, which could be not achieved by chNPL.

Varying the Directionality of Protein Catalysts for Aldol and Retro-Aldol Reactions

Natural aldolase enzymes and created retro-aldolase protein catalysts often catalyze both aldol and retro-aldol reactions depending on the concentrations of the reactants and the products. Here, we report that the directionality of protein catalysts can be altered by replacing one amino acid. The protein catalyst derived from a scaffold of a previously reported retro-aldolase catalyst, catalyzed aldol reactions more efficiently than the previously reported retro-aldolase catalyst. The retro-aldolase catalyst efficiently catalyzed the retro-aldol reaction but was less efficient in catalyzing the aldol reaction. The results indicate that protein catalysts with varying levels of directionality in usually reversibly catalyzed aldol and retro-aldol reactions can be generated from the same protein scaffold.

Recent trends in the stereoselective synthesis of (poly)-substituted 2-oxo acids by biocatalyzed aldol reaction

Recently, an increased interest toward enzymatic carboligation was observed, as biocatalytic carbon-carbon bond formation is a common obstacle in retrosynthetic planning. The construction of extended 2-oxoacid frameworks by 2-oxoacid aldolases and enzymes acting as aldolases is a potent tool for synthetic chemists since a broad spectrum of downstream reactions through functional group interconversions gives access to a plethora of compound classes. In the search for selective biocatalysts, successful protein engineering efforts and high throughput screenings from biodiversity expand the structural diversity of nucleophile and electrophile substrates. Several successful examples with an emphasis on reactions catalyzed by class II aldolases and enzymes acting as class II aldolases are highlighted, including reactions in which both enantiomeric products and in selected cases even diastereomeric products are accessed.

Demystifying the catalytic pathway of Mycobacterium tuberculosis isocitrate lyase

Pulmonary tuberculosis, caused by Mycobacterium tuberculosis, is one of the most persistent diseases leading to death in humans. As one of the key targets during the latent/dormant stage of M. tuberculosis, isocitrate lyase (ICL) has been a subject of interest for new tuberculosis therapeutics. In this work, the cleavage of the isocitrate by M. tuberculosis ICL was studied using quantum mechanics/molecular mechanics method at M06-2X/6-31+G(d,p): AMBER level of theory. The electronic embedding approach was applied to provide a better depiction of electrostatic interactions between MM and QM regions. Two possible pathways (pathway I that involves Asp108 and pathway II that involves Glu182) that could lead to the metabolism of isocitrate was studied in this study. The results suggested that the core residues involved in isocitrate catalytic cleavage mechanism are Asp108, Cys191 and Arg228. A water molecule bonded to Mg²⁺ acts as the catalytic base for the deprotonation of isocitrate C(2)–OH group, while Cys191 acts as the catalytic acid. Our observation suggests that the shuttle proton from isocitrate hydroxyl group C(2) atom is favourably transferred to Asp108 instead of Glu182 with a lower activation energy of 6.2 kcal/mol. Natural bond analysis also demonstrated that pathway I involving the transfer of proton to Asp108 has a higher intermolecular interaction and charge transfer that were associated with higher stabilization energy. The QM/MM transition state stepwise catalytic mechanism of ICL agrees with the in vitro enzymatic assay whereby Asp108Ala and Cys191Ser ICL mutants lost their isocitrate cleavage activities.

Biofilm formation and extracellular microvesicles-The way of foodborne pathogens toward resistance

Almost all known foodborne pathogens are able to form biofilms as one of the strategies for survival under harsh living conditions, to ward off the inhibition and the disinfection during food production, transport and storage, as well as during cleaning and sanitation of corresponding facilities. Biofilms are communities where microbial cells live under constant intracellular interaction and communication. Members of the biofilm community are embedded into extracellular matrix that contains polysaccharides, DNA, lipids, proteins, and small molecules that protect microorganisms and enable their intercellular communication under stress conditions. Membrane vesicles (MVs) are produced by both Gram positive and Gram negative bacteria. These lipid membrane-enveloped nanoparticles play an important role in biofilm genesis and in communication between different biofilm members. Furthermore, MVs are involved in other important steps of bacterial life like cell wall modeling, cellular division, and intercellular communication. They also carry toxins and virulence factors, as well as nucleic acids and different metabolites, and play a key role in host infections. After entering host cells, MVs can start many pathologic processes and cause serious harm and cell death. Prevention and inhibition of both biofilm formation and shedding of MVs by foodborne pathogens has a very important role in food production, storage, and food safety in general. Better knowledge of biofilm formation and maintaining, as well as the role of microbial vesicles in this process and in the process of host cells' infection is essential for food safety and prevention of both food spoilage and host infection.

Benchmarking quantum mechanical methods for calculating reaction energies of reactions catalyzed by enzymes

To assess the accuracy of different quantum mechanical methods for biochemical modeling, the reaction energies of 20 small model reactions (chosen to represent chemical steps catalyzed by commonly studied enzymes) were calculated. The methods tested included several popular Density Functional Theory (DFT) functionals, second-order Møller Plesset perturbation theory (MP2) and its spin-component scaled variant (SCS-MP2), and coupled cluster singles and doubles and perturbative triples (CCSD(T)). Different basis sets were tested. CCSD(T)/aug-cc-pVTZ results for all 20 reactions were used to benchmark the other methods. It was found that MP2 and SCS-MP2 reaction energy calculation results are similar in quality to CCSD(T) (mean absolute error (MAE) of 1.2 and 1.3 kcal mol−1, respectively). MP2 calculations gave a large error in one case, and are more subject to basis set effects, so in general SCS-MP2 calculations are a good choice when CCSD(T) calculations are not feasible. Results with different DFT functionals were of reasonably good quality (MAEs of 2.5–5.1 kcal mol−1), whereas popular semi-empirical methods (AM1, PM3, SCC-DFTB) gave much larger errors (MAEs of 11.6–14.6 kcal mol−1). These results should be useful in guiding methodological choices and assessing the accuracy of QM/MM calculations on enzyme-catalyzed reactions.

Getting zeolite catalysts to play your tune: methyl carboxylate esters as switchable promoters for methanol dehydration to DME

Methyl carboxylate esters can be used as additives to promote the zeolite catalysed formation of dimethly ether (DME) from methanol. By taking advantage of the well-known confinement effect in combination with further functionalisation the potency of the promoter can be markedly enhanced, giving significant increases in DME yield at promoter concentrations as low as 10 ppm relative to methanol. The promotion is readily reversible and promoter concentration can be used to fine tune the zeolite productivity to DME.

Elucidation of a sialic acid metabolism pathway in mucus-foraging Ruminococcus gnavus unravels mechanisms of bacterial adaptation to the gut

Sialic acid (Neu5Ac) is commonly found in terminal location of colonic mucins glycans where it is a much-coveted nutrient for gut bacteria including Ruminococcus gnavus. R. gnavus is part of the healthy gut microbiota in humans but shows a disproportionate representation in diseases. There is therefore a need in understanding the molecular mechanisms underpinning its adaptation to the gut. Previous in vitro work demonstrated that R. gnavus mucin glycan-foraging strategy is strain-dependent and associated with the expression of an intramolecular trans-sialidase releasing 2,7-anhydro-Neu5Ac instead of Neu5Ac from mucins. Here, we have unravelled the metabolism pathway of 2,7-anhydro-Neu5Ac in R. gnavus which is underpinned by the exquisite specificity of the sialic transporter for 2,7-anhydro-Neu5Ac, and by the action of an oxidoreductase converting 2,7-anhydro-Neu5Ac into Neu5Ac which then becomes substrate of a Neu5Ac-specific aldolase. Having generated a R. gnavus nan cluster deletion mutant that lost the ability to grow on sialylated substrates, we showed that in gnotobiotic mice colonised with R. gnavus wild-type and mutant strains, the fitness of the nan mutant was significantly impaired with a reduced ability to colonise the mucus layer. Overall, our study revealed a unique sialic acid pathway in bacteria, with significant implications for the spatial adaptation of mucin-foraging gut symbionts in health and disease.

The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study

Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.

Features and structure of a cold active N-acetylneuraminate lyase

N-acetylneuraminate lyases (NALs) are enzymes that catalyze the reversible cleavage and synthesis of sialic acids. They are therefore commonly used for the production of these high-value sugars. This study presents the recombinant production, together with biochemical and structural data, of the NAL from the psychrophilic bacterium Aliivibrio salmonicida LFI1238 (AsNAL). Our characterization shows that AsNAL possesses high activity and stability at alkaline pH. We confirm that these properties allow for the use in a one-pot reaction at alkaline pH for the synthesis of N-acetylneuraminic acid (Neu5Ac, the most common sialic acid) from the inexpensive precursor N-acetylglucosamine. We also show that the enzyme has a cold active nature with an optimum temperature for Neu5Ac synthesis at 20°C. The equilibrium constant for the reaction was calculated at different temperatures, and the formation of Neu5Ac acid is favored at low temperatures, making the cold active enzyme a well-suited candidate for use in such exothermic reactions. The specific activity is high compared to the homologue from Escherichia coli at three tested temperatures, and the enzyme shows a higher catalytic efficiency and turnover number for cleavage at 37°C. Mutational studies reveal that amino acid residue Asn 168 is important for the high k_cat. The crystal structure of AsNAL was solved to 1.65 Å resolution and reveals a compact, tetrameric protein similar to other NAL structures. The data presented provides a framework to guide further optimization of its application in sialic acid production and opens the possibility for further design of the enzyme.

Enzymatic Synthesis of Trideuterated Sialosides

Sialic acids are a family of acidic monosaccharides often found on the termini of cell surface proteins or lipid glycoconjugates of higher animals. Herein we describe the enzymatic synthesis of the two isotopically labeled sialic acid derivatives d₃-X-Gal-α-2,3-Neu5Ac and d₃-X-Gal-α-2,3-Neu5Gc. Using deuterium oxide as the reaction solvent, deuterium atoms could be successfully introduced during the enzymatic epimerization and aldol addition reactions when the sialosides were generated. NMR and mass spectrometric analyses confirmed that the resulting sialosides were indeed tri-deuterated. These compounds may be of interest as internal standards in liquid chromatography/mass spectrometric assays for biochemical or clinical studies of sialic acids. This was further exemplified by the use of this tri-deuterated sialosides as internal standards for the quantification of sialic acids in meat and egg samples.

Exploration of the Sialic Acid World

Sialic acids are cytoprotectors, mainly localized on the surface of cell membranes with multiple and outstanding cell biological functions. The history of their structural analysis, occurrence, and functions is fascinating and described in this review. Reports from different researchers on apparently similar substances from a variety of biological materials led to the identification of a 9-carbon monosaccharide, which in 1957 was designated "sialic acid." The most frequently occurring member of the sialic acid family is N-acetylneuraminic acid, followed by N-glycolylneuraminic acid and O-acetylated derivatives, and up to now over about 80 neuraminic acid derivatives have been described. They appeared first in the animal kingdom, ranging from echinoderms up to higher animals, in many microorganisms, and are also expressed in insects, but are absent in higher plants. Sialic acids are masks and ligands and play as such dual roles in biology. Their involvement in immunology and tumor biology, as well as in hereditary diseases, cannot be underestimated. N-Glycolylneuraminic acid is very special, as this sugar cannot be expressed by humans, but is a xenoantigen with pathogenetic potential. Sialidases (neuraminidases), which liberate sialic acids from cellular compounds, had been known from very early on from studies with influenza viruses. Sialyltransferases, which are responsible for the sialylation of glycans and elongation of polysialic acids, are studied because of their significance in development and, for instance, in cancer. As more information about the functions in health and disease is acquired, the use of sialic acids in the treatment of diseases is also envisaged.

Pathological macromolecular crystallographic data affected by twinning, partial-disorder and exhibiting multiple lattices for testing of data processing and refinement tools

Twinning is a crystal growth anomaly, which has posed a challenge in macromolecular crystallography (MX) since the earliest days. Many approaches have been used to treat twinned data in order to extract structural information. However, in most cases it is usually simpler to rescreen for new crystallization conditions that yield an untwinned crystal form or, if possible, collect data from non-twinned parts of the crystal. Here, we report 11 structures of engineered variants of the E. coli enzyme N-acetyl-neuraminic lyase which, despite twinning and incommensurate modulation, have been successfully indexed, solved and deposited. These structures span a resolution range of 1.45-2.30 Å, which is unusually high for datasets presenting such lattice disorders in MX and therefore these data provide an excellent test set for improving and challenging MX data processing programs.

Dihydrodipicolinate synthase is absent in fungi

The class I aldolase dihydrodipicolinate synthase (DHDPS) catalyzes the first committed step of the diaminopimelate (DAP) lysine biosynthesis pathway in bacteria, archaea and plants. Despite the existence, in databases, of numerous fungal sequences annotated as DHDPS, its presence in fungi has been the subject of contradictory claims. We report the characterization of DHDPS candidates from fungi. Firstly, the putative DHDPS from Coccidioides immitis (PDB ID: 3QFE) was shown to have negligible enzyme activity. Sequence analysis of 3QFE showed that three out of the seven amino acid residues critical for DHDPS activity are absent; however, exact matches to catalytic residues from two other class I aldolases, 2-keto-3-deoxygluconate aldolase (KDGA), and 4-hydroxy-2-oxoglutarate aldolase (HOGA), were identified. The presence of both KDGA and HOGA activity in 3QFE was confirmed in vitro using enzyme assays, the first report of such dual activity. Subsequent analyses of all publically available fungal sequences revealed that no entry contains all seven residues important for DHDPS function. The candidate with the highest number of identities (6 of 7), KIW77228 from Fonsecaea pedrosoi, was shown to have trace DHDPS activity in vitro, partially restored by substitution of the seventh critical residue, and to be incapable of complementing DHDPS-deficient E. coli cells. Combined with the presence of all seven sequences for the alternative α-aminoadipate (AAA) lysine biosynthesis pathway in C. immitis and F. pedrosoi, we believe that DHDPS and the DAP pathway are absent in fungi, and further, that robust informed methods for annotating genes need to be implemented.

Structural Characterization of the Hydratase-Aldolases, NahE and PhdJ: Implications for the Specificity, Catalysis, and N-Acetylneuraminate Lyase Subgroup of the Aldolase Superfamily

NahE and PhdJ are bifunctional hydratase-aldolases in bacterial catabolic pathways for naphthalene and phenanthrene, respectively. Bacterial species with these pathways can use polycyclic aromatic hydrocarbons (PAHs) as sole sources of carbon and energy. Because of the harmful properties of PAHs and their widespread distribution and persistence in the environment, there is great interest in understanding these degradative pathways, including the mechanisms and specificities of the enzymes found in the pathways. This knowledge can be used to develop and optimize bioremediation techniques. Although hydratase-aldolases catalyze a major step in the PAH degradative pathways, their mechanisms are poorly understood. Sequence analysis identified NahE and PhdJ as members of the N-acetylneuraminate lyase (NAL) subgroup in the aldolase superfamily. Both have a conserved lysine and tyrosine (for Schiff base formation) as well as a GXXGE motif (to bind the pyruvoyl carboxylate group). Herein, we report the structures of NahE, PhdJ, and PhdJ covalently bound to substrate via a Schiff base. Structural analysis and dynamic light scattering experiments show that both enzymes are tetramers. A hydrophobic helix insert, present in the active sites of NahE and PhdJ, might differentiate them from other NAL subgroup members. The individual specificities of NahE and PhdJ are governed by Asn-281/Glu-285 and Ser-278/Asp-282, respectively. Finally, the PhdJ complex structure suggests a potential mechanism for hydration of substrate and subsequent retro-aldol fission. The combined findings fill a gap in our mechanistic understanding of these enzymes and their place in the NAL subgroup.

Multiscale analysis of enantioselectivity in enzyme-catalysed 'lethal synthesis' using projector-based embedding

The action of fluoroacetate as a broad-spectrum mammalian pesticide depends on the 'lethal synthesis' of fluorocitrate by citrate synthase, through a subtle enantioselective enolization of fluoroacetyl-coenzyme A. In this work, we demonstrate how a projection-based embedding method can be applied to calculate coupled cluster (CCSD(T)) reaction profiles from quantum mechanics/molecular mechanics optimized pathways for this enzyme reaction. Comparison of pro-R and pro-S proton abstraction in citrate synthase at the CCSD(T)-in-DFT//MM level yields the correct enantioselectivity. We thus demonstrate the potential of projection-based embedding for determining stereoselectivity in enzymatic systems. We further show that the method is simple to apply, eliminates variability due to the choice of density functional theory functional and allows the efficient calculation of CCSD(T) quality enzyme reaction barriers.

Extending enzyme molecular recognition with an expanded amino acid alphabet

Natural enzymes are constructed from the 20 proteogenic amino acids, which may then require posttranslational modification or the recruitment of coenzymes or metal ions to achieve catalytic function. Here, we demonstrate that expansion of the alphabet of amino acids can also enable the properties of enzymes to be extended. A chemical mutagenesis strategy allowed a wide range of noncanonical amino acids to be systematically incorporated throughout an active site to alter enzymic substrate specificity. Specifically, 13 different noncanonical side chains were incorporated at 12 different positions within the active site of N-acetylneuraminic acid lyase (NAL), and the resulting chemically modified enzymes were screened for activity with a range of aldehyde substrates. A modified enzyme containing a 2,3-dihydroxypropyl cysteine at position 190 was identified that had significantly increased activity for the aldol reaction of erythrose with pyruvate compared with the wild-type enzyme. Kinetic investigation of a saturation library of the canonical amino acids at the same position showed that this increased activity was not achievable with any of the 20 proteogenic amino acids. Structural and modeling studies revealed that the unique shape and functionality of the noncanonical side chain enabled the active site to be remodeled to enable more efficient stabilization of the transition state of the reaction. The ability to exploit an expanded amino acid alphabet can thus heighten the ambitions of protein engineers wishing to develop enzymes with new catalytic properties.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

Structure and inhibition of N-acetylneuraminate lyase from methicillin-resistant Staphylococcus aureus

N-Acetylneuraminate lyase is the first committed enzyme in the degradation of sialic acid by bacterial pathogens. In this study, we analyzed the kinetic parameters of N-acetylneuraminate lyase from methicillin-resistant Staphylococcus aureus (MRSA). We determined that the enzyme has a relatively high K_M of 3.2 mm, suggesting that flux through the catabolic pathway is likely to be controlled by this enzyme. Our data indicate that sialic acid alditol, a known inhibitor of N-acetylneuraminate lyase enzymes, is a stronger inhibitor of MRSA N-acetylneuraminate lyase than of Clostridium perfringens N-acetylneuraminate lyase. Our analysis of the crystal structure of ligand-free and 2R-sialic acid alditol-bound MRSA N-acetylneuraminate lyase suggests that subtle dynamic differences in solution and/or altered binding interactions within the active site may account for species-specific inhibition.

A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase

Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.

Evaluation of fluoropyruvate as nucleophile in reactions catalysed by N-acetyl neuraminic acid lyase variants: scope, limitations and stereoselectivity

The catalysis of reactions involving fluoropyruvate as donor by N-acetyl neuraminic acid lyase (NAL) variants was investigated. Under kinetic control, the wild-type enzyme catalysed the reaction between fluoropyruvate and N-acetyl mannosamine to give a 90 : 10 ratio of the (3R,4R)- and (3S,4R)-configured products; after extended reaction times, equilibration occurred to give a 30 : 70 mixture of these products. The efficiency and stereoselectivity of reactions of a range of substrates catalysed by the E192N, E192N/T167V/S208V and E192N/T167G NAL variants were also studied. Using fluoropyruvate and (2R,3S)- or (2S,3R)-2,3-dihydroxy-4-oxo-N,N-dipropylbutanamide as substrates, it was possible to obtain three of the four possible diastereomeric products; for each product, the ratio of anomeric and pyranose/furanose forms was determined. The crystal structure of S. aureus NAL in complex with fluoropyruvate was determined, assisting rationalisation of the stereochemical outcome of C-C bond formation.

Direct synthesis of C-glycosides from unprotected 2-N-acyl-aldohexoses via aldol condensation–oxa-Michael reactions with unactivated ketones

C-glycosides were synthesized from unprotected 2-N-acyl-aldohexoses and unactivated ketones in one pot via aldol condensation–oxa-Michael reactions.

Probing the Catalytic Mechanism Involved in the Isocitrate Lyase Superfamily: Hybrid Quantum Mechanical/Molecular Mechanical Calculations on 2,3-Dimethylmalate Lyase

The isocitrate lyase (ICL) superfamily catalyzes the cleavage of the C(2)-C(3) bond of various α-hydroxy acid substrates. Members of the family are found in bacteria, fungi, and plants and include ICL itself, oxaloacetate hydrolase (OAH), 2-methylisocitrate lyase (MICL), and (2R,3S)-dimethylmalate lyase (DMML) among others. ICL and related targets have been the focus of recent studies to treat bacterial and fungal infections, including tuberculosis. The catalytic process by which this family achieves C(2)-C(3) bond breaking is still not clear. Extensive structural studies have been performed on this family, leading to a number of plausible proposals for the catalytic mechanism. In this paper, we have applied quantum mechanical/molecular mechanical (QM/MM) methods to the most recently reported family member, DMML, to assess whether any of the mechanistic proposals offers a clear energetic advantage over the others. Our results suggest that Arg161 is the general base in the reaction and Cys124 is the general acid, giving rise to a rate-determining barrier of approximately 10 kcal/mol.

Characterization of a novel N-acetylneuraminic acid lyase favoring industrial N-acetylneuraminic acid synthesis

N-Acetylneuraminic acid lyase (NAL, E.C. number 4.1.3.3) is a Class I aldolase that catalyzes the reversible aldol cleavage of N-acetylneuraminic acid (Neu5Ac) from pyruvate and N-acetyl-D-mannosamine (ManNAc). Due to the equilibrium favoring Neu5Ac cleavage, the enzyme catalyzes the rate-limiting step of two biocatalytic reactions producing Neu5Ac in industry. We report the biochemical characterization of a novel NAL from a “GRAS” (General recognized as safe) strain C. glutamicum ATCC 13032 (CgNal). Compared to all previously reported NALs, CgNal exhibited the lowest k_cat/K_m value for Neu5Ac and highest k_cat/K_m values for ManNAc and pyruvate, which makes CgNal favor Neu5Ac synthesis the most. The recombinant CgNal reached the highest expression level (480 mg/L culture), and the highest reported yield of Neu5Ac was achieved (194 g/L, 0.63 M). All these unique properties make CgNal a promising biocatalyst for industrial Neu5Ac biosynthesis. Additionally, although showing the best Neu5Ac synthesis activity among the NAL family, CgNal is more related to dihydrodipicolinate synthase (DHDPS) by phylogenetic analysis. The activities of CgNal towards both NAL's and DHDPS' substrates are fairly high, which indicates CgNal a bi-functional enzyme. The sequence analysis suggests that CgNal might have adopted a unique set of residues for substrates recognition.

Structure and function of a decarboxylating Agrobacterium tumefaciens keto-deoxy-d-galactarate dehydratase

Agrobacterium tumefaciens (At) strain C58 contains an oxidative enzyme pathway that can function on both d-glucuronic and d-galacturonic acid. The corresponding gene coding for At keto-deoxy-d-galactarate (KDG) dehydratase is located in the same gene cluster as those coding for uronate dehydrogenase (At Udh) and galactarolactone cycloisomerase (At Gci) which we have previously characterized. Here, we present the kinetic characterization and crystal structure of At KDG dehydratase, which catalyzes the next step, the decarboxylating hydrolyase reaction of KDG to produce α-ketoglutaric semialdehyde (α-KGSA) and carbon dioxide. The crystal structures of At KDG dehydratase and its complexes with pyruvate and 2-oxoadipic acid, two substrate analogues, were determined to 1.7 Å, 1.5 Å, and 2.1 Å resolution, respectively. Furthermore, mass spectrometry was used to confirm reaction end-products. The results lead us to propose a structure-based mechanism for At KDG dehydratase, suggesting that while the enzyme belongs to the Class I aldolase protein family, it does not follow a typical retro-aldol condensation mechanism.

First functional and mutational analysis of group 3 N-acetylneuraminate lyases from Lactobacillus antri and Lactobacillus sakei 23K

N-acetyl neuraminate lyases (NALs) catalyze the reversible aldol cleavage of N-acetyl neuraminic acid (Neu5Ac) to pyruvate and N-acetyl-D-mannosamine (ManNAc). Previous phylogenetic studies divided NALs into four different groups. Groups 1 and 2 have been well characterized at both kinetic and molecular levels, but no NAL from group 3 has been studied to date. In this work, a functional characterization of two group 3 members was performed using the recombinant NALs from Lactobacillus antri and Lactobacillus sakei 23K, revealing an optimal pH of between 6.0 and 7.0, low stability at basic pHs (>8.0), low optimal temperatures and, especially, low catalytic efficiency compared with their counterparts in group 1 and 2. The mutational analysis carried out showed that a plausible molecular reason for the low activity shown by Lactobacillus antri and Lactobacillus sakei 23k NALs compared with group 1 and 2 NALs could be the relatively small sugar-binding pocket they contain. A functional divergence analysis concluding that group 3 is more closely related to group 2 than to group 1.

Reactions of pyruvates: organocatalytic synthesis of functionalized dihydropyrans in one pot and further transformations to functionalized carbocycles and heterocycles

Concise cascade reactions of pyruvates with aldehydes that generate functionalized dihydropyran derivatives in one pot have been developed. The products, dihydropyrans, were further concisely transformed to various functionalized molecules.