The Medium-Chain Dehydrogenase/reductase Engineering Database: a systematic analysis of a diverse protein family to understand sequence-structure-function relationship

Stairway to Stereoisomers: Engineering Short- and Medium-Chain Ketoreductases To Produce Chiral Alcohols

The short- and medium-chain dehydrogenase/reductase superfamilies are responsible for most chiral alcohol production in laboratories and industries. In nature, they participate in diverse roles such as detoxification, housekeeping, secondary metabolite production, and catalysis of several chemicals with commercial and environmental significance. As a result, they are used in industries to create biopolymers, active pharmaceutical intermediates (APIs), and are also used as components of modular enzymes like polyketide synthases for fabricating bioactive molecules. Consequently, random, semi-rational and rational engineering have helped transform these enzymes into product-oriented efficient catalysts. The rise of newer synthetic chemicals and their enantiopure counterparts has proved challenging, and engineering them has been the subject of numerous studies. However, they are frequently limited to the synthesis of a single chiral alcohol. The study attempts to defragment and describe hotspots of engineering short- and medium-chain dehydrogenases/reductases for the production of chiral synthons.

Utilizing Alcohol for Alkane Biosynthesis by Introducing a Fatty Alcohol Dehydrogenase

Alkanes produced by microorganisms are expected to be an alternative to fossil fuels as an energy source. Microbial synthesis of alkanes involves the formation of fatty aldehydes via fatty acyl coenzyme A (acyl-CoA) intermediates derived from fatty acid metabolism, followed by aldehyde decarbonylation to generate alkanes. Advancements in metabolic engineering have enabled the construction of such pathways in various microorganisms, including Escherichia coli. However, endogenous aldehyde reductases in the host microorganisms are highly active in converting fatty aldehydes to fatty alcohols, limiting the substrate pool for alkane production. To reuse the alcohol by-product, a screening of fatty alcohol-assimilating microorganisms was conducted, and a bacterial strain, Pantoea sp. strain 7-4, was found to convert 1-tetradecanol to tetradecanal. From this strain, an alcohol dehydrogenase, PsADH, was purified and found to be involved in 1-tetradecanol-oxidizing reaction. Subsequent heterologous expression of the PsADH gene in E. coli was conducted, and recombinant PsADH was purified for a series of biochemical characterizations, including cofactors, optimal reaction conditions, and kinetic parameters. Furthermore, direct alkane production from alcohol was achieved in E. coli by coexpressing PsADH with a cyanobacterial aldehyde-deformylating oxygenase and a reducing system, including ferredoxin and ferredoxin reductase, from Nostoc punctiforme PCC73102. The alcohol-aldehyde-alkane synthetic route established in this study will provide a new approach to utilizing fatty alcohols for the production of alkane biofuel. IMPORTANCE Alcohol dehydrogenases are a group of enzymes found in many organisms. Unfortunately, studies on these enzymes mainly focus on their activities toward short-chain alcohols. In this study, we discovered an alcohol dehydrogenase, PsADH, from the bacterium Pantoea sp. 7-4, which can oxidize 1-tetradecanol to tetradecanal. The medium-chain aldehyde products generated by this enzyme can serve as the substrate of aldehyde-deformylating oxygenase to produce alkanes. The enzyme found in this study can be applied to the biosynthetic pathway involving the formation of medium-chain aldehydes to produce alkanes and other valuable compounds.

Advanced Insights into Catalytic and Structural Features of the Zinc-Dependent Alcohol Dehydrogenase from Thauera aromatica

The asymmetric reduction of ketones to chiral hydroxyl compounds by alcohol dehydrogenases (ADHs) is an established strategy for the provision of valuable precursors for fine chemicals and pharmaceutics. However, most ADHs favor linear aliphatic and aromatic carbonyl compounds, and suitable biocatalysts with preference for cyclic ketones and diketones are still scarce. Among the few candidates, the alcohol dehydrogenase from Thauera aromatica (ThaADH) stands out with a high activity for the reduction of the cyclic α‐diketone 1,2‐cyclohexanedione to the corresponding α‐hydroxy ketone. This study elucidates catalytic and structural features of the enzyme. ThaADH showed a remarkable thermal and pH stability as well as stability in the presence of polar solvents. A thorough description of the substrate scope combined with the resolution and description of the crystal structure, demonstrated a strong preference of ThaADH for cyclic α‐substituted cyclohexanones, and indicated structural determinants responsible for the unique substrate acceptance.

Engineered Corynebacterium glutamicum as the Platform for the Production of Aromatic Aldehydes

Aromatic aldehydes, including 4-hydroxybenzaldehyde (4-HB aldehyde), protocatechuic (PC) aldehyde, and vanillin, are used as important flavors, fragrances, and pharmaceutical precursors and have several biological and therapeutic effects. Production of aromatic aldehydes in microbial hosts poses a challenge due to its rapid and endogenous reduction to alcohols. To address this hurdle, prospecting of the genome of Corynebacterium glutamicum yielded 27 candidate proteins that were used in comprehensive screening with a 4-hydroxybenzyl (4-HB) alcohol–producing strain. We identified that NCgl0324 has aromatic aldehyde reductase activity and contributed to 4-HB aldehyde reduction in vivo since the NCgl0324 deletion strain HB-Δ0324 produced 1.36 g/L of 4-HB aldehyde, that is, about 188% more than its parental strain. To demonstrate that NCgl0324 knockout can also improve production of PC aldehyde and vanillin, first, a basal MA303 strain that produces protocatechuate was engineered from 4-hydroxybenzoate-synthesizing C. glutamicum APS963, followed by deletion of NCgl0324 to generate PV-Δ0324. The PC aldehyde/alcohol or vanillin/vanillyl alcohol biosynthetic pathways, respectively, were able to be expanded from protocatechuate upon introduction of carboxylic acid reductase (CAR) and catechol O-methyltransferase encoded by a mutated comt^m gene. In shake flask culture, the resulting NCgl0324 deletion strains PV-IΔ0324 and PV-IYΔ0324 were shown to produce 1.18 g/L PC aldehyde and 0.31 g/L vanillin, respectively. Thus, modulation of the identified NCgl0324 gene was shown to have the potential to boost production of valuable aromatic aldehydes and alcohols.

Microbial Hydroxysteroid Dehydrogenases: From Alpha to Omega

Bile acids (BAs) and glucocorticoids are steroid hormones derived from cholesterol that are important signaling molecules in humans and other vertebrates. Hydroxysteroid dehydrogenases (HSDHs) are encoded both by the host and by their resident gut microbiota, and they reversibly convert steroid hydroxyl groups to keto groups. Pairs of HSDHs can reversibly epimerize steroids from α-hydroxy conformations to β-hydroxy, or β-hydroxy to ω-hydroxy in the case of ω-muricholic acid. These reactions often result in products with drastically different physicochemical properties than their precursors, which can result in steroids being activators or inhibitors of host receptors, can affect solubility in fecal water, and can modulate toxicity. Microbial HSDHs modulate sterols associated with diseases such as colorectal cancer, liver cancer, prostate cancer, and polycystic ovary syndrome. Although the role of microbial HSDHs is not yet fully elucidated, they may have therapeutic potential as steroid pool modulators or druggable targets in the future. In this review, we explore metabolism of BAs and glucocorticoids with a focus on biotransformation by microbial HSDHs.

Phylogenetics-based identification and characterization of a superior 2,3-butanediol dehydrogenase for Zymomonas mobilis expression

BACKGROUND

Zymomonas mobilis has recently been shown to be capable of producing the valuable platform biochemical, 2,3-butanediol (2,3-BDO). Despite this capability, the production of high titers of 2,3-BDO is restricted by several physiological parameters. One such bottleneck involves the conversion of acetoin to 2,3-BDO, a step catalyzed by 2,3-butanediol dehydrogenase (Bdh). Several Bdh enzymes have been successfully expressed in Z. mobilis, although a highly active enzyme is yet to be identified for expression in this host. Here, we report the application of a phylogenetic approach to identify and characterize a superior Bdh, followed by validation of its structural attributes using a mutagenesis approach.

RESULTS

Of the 11 distinct bdh genes that were expressed in Z. mobilis, crude extracts expressing Serratia marcescens Bdh (SmBdh) were found to have the highest activity (8.89 µmol/min/mg), when compared to other Bdh enzymes (0.34-2.87 µmol/min/mg). The SmBdh crystal structure was determined through crystallization with cofactor (NAD⁺) and substrate (acetoin) molecules bound in the active site. Active SmBdh was shown to be a tetramer with the active site populated by a Gln247 residue contributed by the diagonally opposite subunit. SmBdh showed a more extensive supporting hydrogen-bond network in comparison to the other well-studied Bdh enzymes, which enables improved substrate positioning and substrate specificity. This protein also contains a short α6 helix, which provides more efficient entry and exit of molecules from the active site, thereby contributing to enhanced substrate turnover. Extending the α6 helix to mimic the lower activity Enterobacter cloacae (EcBdh) enzyme resulted in reduction of SmBdh function to nearly 3% of the total activity. In great contrast, reduction of the corresponding α6 helix of the EcBdh to mimic the SmBdh structure resulted in ~ 70% increase in its activity.

CONCLUSIONS

This study has demonstrated that SmBdh is superior to other Bdhs for expression in Z. mobilis for 2,3-BDO production. SmBdh possesses unique structural features that confer biochemical advantage to this protein. While coordinated active site formation is a unique structural characteristic of this tetrameric complex, the smaller α6 helix and extended hydrogen network contribute towards improved activity and substrate promiscuity of the enzyme.

The influence of polycyclic aromatic hydrocarbons in protein profile of Medicago sativa L

Medicago sativa L. (alfalfa) are studied as potential phytoremediation agents of priority pollutants like polycyclic aromatic hydrocarbons (PAH). However, elucidation of the biochemical mechanisms involved in phytoremediation is a topic to be explored with knowledge gaps. This study aims to identify and classify proteins expressed in the aerial parts of laboratory-cultivated alfalfa in the presence and absence of pyrene, anthracene, and phenanthrene. Soil samples were amended with 100 mg.kg^-1 of each PAH (total concentration of 300 ppm) and cultivated with alfalfa plants for 20 days. After this, aerial parts of cultivated plants from each condition were collected for qualitative proteomic analysis (ESI-Q/TOF). The results showed a significant increase (Student's t-test p < 0.05) of 41.7% in the concentration of proteins from plants grown in PAH-amended substrates, changes in the protein profile, with intense protein bands observed at 40-55, 34, 28, and 15 kDa when compared to the control. A total of 504 proteins were identified and classified into 12 functional categories, highlighting the identification of 11 phytoremediation-related proteins candidates in plants grown in the presence of PAH, with biological functions related to diverse metabolisms involved in the xenobiotics biodegradation (included PAH), glutathione and response to stress.

Substitution of cysteine-153 ligated to the catalytic zinc in yeast alcohol dehydrogenase with aspartic acid and analysis of mechanisms of related medium chain dehydrogenases

The catalytic zincs in complexes of horse liver and yeast alcohol dehydrogenases (ADH) with NAD⁺ and the substrate analogue, 2,2,2-trifluoroethanol, are ligated to two cysteine residues and one histidine residue from the protein and the oxygen from the alcohol. The zinc facilitates deprotonation of the alcohol and is essential for catalysis. In the yeast apoenzyme, the zinc is coordinated to a nearby glutamic acid, which is displaced by the alcohol in the complex with NAD⁺. Some homologous medium chain dehydrogenases have a cysteine replaced by aspartic or glutamic acid residues. How an aspartic acid would affect catalysis was studied by replacing Cys-153 in Saccharomyces cerevisiae ADH1 by using site-directed mutagenesis. The C153D enzyme was about as stable as the wild-type enzyme, if EDTA was not included in the buffers. The substitution increased affinity for NAD⁺ by 3-fold, but did not affect NADH binding. At pH 7.3, the turnover number for ethanol oxidation (V₁/E_t) decreased by 7-fold and catalytic efficiency decreased 18-fold (V₁/E_tK_b), but turnover for acetaldehyde reduction (V₂/E_t) was the same as for wild-type enzyme and catalytic efficiency decreased 8-fold (V₂/E_tK_p). Deuterium isotope effects of 3.0 on V₁/E_t and 3.8 on V₁/E_tK_b for ethanol oxidation suggest that hydride transfer is more rate-limiting for turnover for the C153D enzyme than by wild-type enzyme. The patterns of pH dependence for V₁/E_tK_b for ethanol oxidation were similar for both enzymes in the pH range from 7 to 9. The C153D substitution decreased binding of trifluoroethanol by 5-fold and of pyrazole by 65-fold. Substrate specificities for C153D and wild-type ADHs for primary alcohols have similar patterns. Efficiency for secondary alcohols decreased only about 4-fold, and efficiencies for 1,2-propanediol and acetone were about the same as for wild-type enzyme. The C153D substitution modestly affects catalysis by altering ligand exchange on the zinc or local structure. Structures and mechanisms for acid-base catalysis in related medium chain dehydrogenases with substitutions of the homologous cysteine are reviewed and analyzed.

Biosynthesis of ilamycins featuring unusual building blocks and engineered production of enhanced anti-tuberculosis agents

Tuberculosis remains one of the world’s deadliest communicable diseases, novel anti-tuberculosis agents are urgently needed due to severe drug resistance and the co-epidemic of tuberculosis/human immunodeficiency virus. Here, we show the isolation of six anti-mycobacterial ilamycin congeners (1–6) bearing rare L-3-nitro-tyrosine and L-2-amino-4-hexenoic acid structural units from the deep sea-derived Streptomyces atratus SCSIO ZH16. The biosynthesis of the rare L-3-nitrotyrosine and L-2-amino-4-hexenoic acid units as well as three pre-tailoring and two post-tailoring steps are probed in the ilamycin biosynthetic machinery through a series of gene inactivation, precursor chemical complementation, isotope-labeled precursor feeding experiments, as well as structural elucidation of three intermediates (6–8) from the respective mutants. Most impressively, ilamycins E₁/E₂, which are produced in high titers by a genetically engineered mutant strain, show very potent anti-tuberculosis activity with an minimum inhibitory concentration value ≈9.8 nM to Mycobacterium tuberculosis H37Rv constituting extremely potent and exciting anti-tuberculosis drug leads.

Tuberculosis (TB) remains one of the world’s deadliest communicable diseases, novel anti-TB agents are urgently needed due to severe drug resistance and the co-epidemic of TB/HIV. Here, the authors show that anti-mycobacterial ilamycin congeners bearing unusual structural units possess extremely potent anti-tuberculosis activities.

Clade classification of monolignol biosynthesis gene family members reveals target genes to decrease lignin in Lolium perenne

In monocots, lignin content has a strong impact on the digestibility of the cell wall fraction. Engineering lignin biosynthesis requires a profound knowledge of the role of paralogues in the multigene families that constitute the monolignol biosynthesis pathway. We applied a bioinformatics approach for genome-wide identification of candidate genes in Lolium perenne that are likely to be involved in the biosynthesis of monolignols. More specifically, we performed functional subtyping of phylogenetic clades in four multigene families: 4CL, COMT, CAD and CCR. Essential residues were considered for functional clade delineation within these families. This classification was complemented with previously published experimental evidence on gene expression, gene function and enzymatic activity in closely related crops and model species. This allowed us to assign functions to novel identified L. perenne genes, and to assess functional redundancy among paralogues. We found that two 4CL paralogues, two COMT paralogues, three CCR paralogues and one CAD gene are prime targets for genetic studies to engineer developmentally regulated lignin in this species. Based on the delineation of sequence conservation between paralogues and a first analysis of allelic diversity, we discuss possibilities to further study the roles of these paralogues in lignin biosynthesis, including expression analysis, reverse genetics and forward genetics, such as association mapping. We propose criteria to prioritise paralogues within multigene families and certain SNPs within these genes for developing genotyping assays or increasing power in association mapping studies. Although L. perenne was the target of the analyses presented here, this functional subtyping of phylogenetic clades represents a valuable tool for studies investigating monolignol biosynthesis genes in other monocot species.

Yeast alcohol dehydrogenase structure and catalysis

Yeast (Saccharomyces cerevisiae) alcohol dehydrogenase I (ADH1) is the constitutive enzyme that reduces acetaldehyde to ethanol during the fermentation of glucose. ADH1 is a homotetramer of subunits with 347 amino acid residues. A structure for ADH1 was determined by X-ray crystallography at 2.4 Å resolution. The asymmetric unit contains four different subunits, arranged as similar dimers named AB and CD. The unit cell contains two different tetramers made up of "back-to-back" dimers, AB:AB and CD:CD. The A and C subunits in each dimer are structurally similar, with a closed conformation, bound coenzyme, and the oxygen of 2,2,2-trifluoroethanol ligated to the catalytic zinc in the classical tetrahedral coordination with Cys-43, Cys-153, and His-66. In contrast, the B and D subunits have an open conformation with no bound coenzyme, and the catalytic zinc has an alternative, inverted coordination with Cys-43, Cys-153, His-66, and the carboxylate of Glu-67. The asymmetry in the dimeric subunits of the tetramer provides two structures that appear to be relevant for the catalytic mechanism. The alternative coordination of the zinc may represent an intermediate in the mechanism of displacement of the zinc-bound water with alcohol or aldehyde substrates. Substitution of Glu-67 with Gln-67 decreases the catalytic efficiency by 100-fold. Previous studies of structural modeling, evolutionary relationships, substrate specificity, chemical modification, and site-directed mutagenesis are interpreted more fully with the three-dimensional structure.

The triterpene cyclase protein family: a systematic analysis

Triterpene cyclases catalyze a broad range of cyclization reactions to form polycyclic triterpenes. Triterpene cyclases that convert squalene to hopene are named squalene-hopene cyclases (SHC) and triterpene cyclases that convert oxidosqualene are named oxidosqualene cyclases (OSC). Many sequences have been published, but there is only one structure available for each of SHCs and OSCs. Although they catalyze a similar reaction, the sequence similarity between SHCs and OSCs is low. A family classification based on phylogenetic analysis revealed 20 homologous families which are grouped into two superfamilies, SHCs and OSCs. Based on this family assignment, the Triterpene Cyclase Engineering Database (TTCED) was established. It integrates available information on sequence and structure of 639 triterpene cyclases as well as on structurally and functionally relevant amino acids. Family specific multiple sequence alignments were generated to identify the functionally relevant residues. Based on sequence alignments, conserved residues in SHCs and OSCs were analyzed and compared to experimentally confirmed mutational data. Functional schematic models of the central cavities of OSCs and SHCs were derived from structure comparison and sequence conservation analysis. These models demonstrate the high similarity of the substrate binding cavity of SHCs and OSCs and the equivalences of the respective residues. The TTCED is a novel source for comprehensive information on the triterpene cyclase family, including a compilation of previously described mutational data. The schematic models present the conservation analysis in a readily available fashion and facilitate the correlation of residues to a specific function or substrate interaction.

Geraniol and geranial dehydrogenases induced in anaerobic monoterpene degradation by Castellaniella defragrans

Castellaniella defragrans is a Betaproteobacterium capable of coupling the oxidation of monoterpenes with denitrification. Geraniol dehydrogenase (GeDH) activity was induced during growth with limonene in comparison to growth with acetate. The N-terminal sequence of the purified enzyme directed the cloning of the corresponding open reading frame (ORF), the first bacterial gene for a GeDH (geoA, for geraniol oxidation pathway). The C. defragrans geraniol dehydrogenase is a homodimeric enzyme that affiliates with the zinc-containing benzyl alcohol dehydrogenases in the superfamily of medium-chain-length dehydrogenases/reductases (MDR). The purified enzyme most efficiently catalyzes the oxidation of perillyl alcohol (k(cat)/K(m) = 2.02 × 10(6) M(-1) s(-1)), followed by geraniol (k(cat)/K(m) = 1.57 × 10(6) M(-1) s(-1)). Apparent K(m) values of <10 μM are consistent with an in vivo toxicity of geraniol above 5 μM. In the genetic vicinity of geoA is a putative aldehyde dehydrogenase that was named geoB and identified as a highly abundant protein during growth with phellandrene. Extracts of Escherichia coli expressing geoB demonstrated in vitro a geranial dehydrogenase (GaDH) activity. GaDH activity was independent of coenzyme A. The irreversible formation of geranic acid allows for a metabolic flux from β-myrcene via linalool, geraniol, and geranial to geranic acid.

Identification of important residues in diketoreductase from Acinetobacter baylyi by molecular modeling and site-directed mutagenesis

Diketoreductase (DKR) from Acinetobacter baylyi exhibits a unique property of double reduction of a β, δ-diketo ester with excellent stereoselectivity, which can serve as an efficient biocatalyst for the preparation of an important chiral intermediate for cholesterol lowering statin drugs. Taken the advantage of high homology between DKR and human heart 3-hydroxyacyl-CoA dehydrogenase (HAD), a molecular model was created to compare the tertiary structures of DKR and HAD. In addition to the possible participation of His-143 in the enzyme catalysis by pH profile, three key amino acid residues, Ser-122, His-143 and Glu-155, were identified and mutated to explore the possibility of involving in the catalytic process. The catalytic activities for mutants S122A/C, H143A/K and E155Q were below detectable level, while their binding affinities to the diketo ester substrate and cofactor NADH did not change obviously. The experimental results were further supported by molecular docking, suggesting that Ser-122 and His-143 were essential for the proton transfer to the carbonyl functional groups of the substrate. Moreover, Glu-155 was crucial for maintaining the proper orientation and protonation of the imidazole ring of His-143 for efficient catalysis.

The isoelectric region of proteins: a systematic analysis

Background

Binding of proteins in ion exchange chromatography is dominated by electrostatic interactions and can be tuned by adjusting pH and ionic strength of the solvent. Therefore, the isoelectric region (IER), the pH region of almost zero charge near the pI, has been used to predict the binding properties of proteins.

Principal findings

Usually the IER is small and binding and elution is carried out at pH values near to the pI. However, some proteins with an extended IER have been shown to bind and elute far away from its pI. To analyze factors that mediate the size of the IER and to identify proteins with an extended IER, two protein families consisting of more than 7000 proteins were systematically investigated. Most proteins were found to have a small IER and thus are expected to bind or elute near to their pI, while only a small fraction of less than 2% had a large IER.

Conclusions

Only four factors, the number of histidines, the pI, the number of titratable amino acids and the ratio of acidic to basic residues, are sufficient to reliably classify proteins by their IER based on their sequence only, and thus to predict their binding and elution behaviour in ion exchange chromatography.

The cytochrome P450 engineering database: Integration of biochemical properties

Background

Cytochrome P450 monooxygenases (CYPs) form a vast and diverse enzyme class of particular interest in drug development and a high biotechnological potential. Although very diverse in sequence, they share a common structural fold. For the comprehensive and systematic comparison of protein sequences and structures the Cytochrome P450 Engineering Database (CYPED) was established. It was built up based on an extensible data model that enables its functions readily enhanced.

Description

The new version of the CYPED contains information on sequences and structures of 8613 and 47 proteins, respectively, which strictly follow Nelson's classification rules for homologous families and superfamilies. To gain biochemical information on substrates and inhibitors, the CYPED was linked to the Cytochrome P450 Knowledgebase (CPK). To overcome differences in the data model and inconsistencies in the content of CYPED and CPK, a metric was established based on sequence similarity to link protein sequences as primary keys. In addition, the annotation of structurally and functionally relevant residues was extended by a reliable prediction of conserved secondary structure elements and by information on the effect of single nucleotide polymorphisms.

Conclusion

The online accessible version of the CYPED at http://www.cyped.uni-stuttgart.de provides a valuable tool for the analysis of sequences, structures and their relationships to biochemical properties.

The PHA Depolymerase Engineering Database: A systematic analysis tool for the diverse family of polyhydroxyalkanoate (PHA) depolymerases

Background

Polyhydroxyalkanoates (PHAs) can be degraded by many microorganisms using intra- or extracellular PHA depolymerases. PHA depolymerases are very diverse in sequence and substrate specificity, but share a common α/β-hydrolase fold and a catalytic triad, which is also found in other α/β-hydrolases.

Results

The PHA Depolymerase Engineering Database (DED, ) has been established as a tool for systematic analysis of this enzyme family. The DED contains sequence entries of 587 PHA depolymerases, which were assigned to 8 superfamilies and 38 homologous families based on their sequence similarity. For each family, multiple sequence alignments and profile hidden Markov models are provided, and functionally relevant residues are annotated.

Conclusion

The DED is a valuable tool which can be applied to identify new PHA depolymerase sequences from complete genomes in silico, to classify PHA depolymerases, to predict their biochemical properties, and to design enzyme variants with improved properties.