Mutations of an NAD(P)H-dependent flavoprotein monooxygenase that influence cofactor promiscuity and enantioselectivity

Genome analysis for the identification of genes involved in phenanthrene biodegradation pathway in Stenotrophomonas indicatrix CPHE1. Phenanthrene mineralization in soils assisted by integrated approaches

Phenanthrene (PHE) is a highly toxic compound, widely present in soils. For this reason, it is essential to remove PHE from the environment. Stenotrophomonas indicatrix CPHE1 was isolated from an industrial soil contaminated by polycyclic aromatic hydrocarbons (PAHs) and was sequenced to identify the PHE degrading genes. Dioxygenase, monooxygenase, and dehydrogenase gene products annotated in S. indicatrix CPHE1 genome were clustered into different trees with reference proteins. Moreover, S. indicatrix CPHE1 whole-genome sequences were compared to genes of PAHs-degrading bacteria retrieved from databases and literature. On these basis, reverse transcriptase-polymerase chain reaction (RT-PCR) analysis pointed out that cysteine dioxygenase (cysDO), biphenyl-2,3-diol 1,2-dioxygenase (bphC), and aldolase hydratase (phdG) were expressed only in the presence of PHE. Therefore, different techniques have been designed to improve the PHE mineralization process in five PHE artificially contaminated soils (50 mg kg^-1), including biostimulation, adding a nutrient solution (NS), bioaugmentation, inoculating S. indicatrix CPHE1 which was selected for its PHE-degrading genes, and the use of 2-hydroxypropyl-β-cyclodextrin (HPBCD) as a bioavailability enhancer. High percentages of PHE mineralization were achieved for the studied soils. Depending on the soil, different treatments resulted to be successful; in the case of a clay loam soil, the best strategy was the inoculation of S. indicatrix CPHE1 and NS (59.9% mineralized after 120 days). In sandy soils (CR and R soils) the highest percentage of mineralization was achieved in presence of HPBCD and NS (87.3% and 61.3%, respectively). However, the combination of CPHE1 strain, HPBCD, and NS showed to be the most efficient strategy for sandy and sandy loam soils (LL and ALC soils showed 35% and 74.6%, respectively). The results indicated a high degree of correlation between gene expression and the rates of mineralization.

Evolution of enzyme functionality in the flavin-containing monooxygenases

Among the molecular mechanisms of adaptation in biology, enzyme functional diversification is indispensable. By allowing organisms to expand their catalytic repertoires and adopt fundamentally different chemistries, animals can harness or eliminate new-found substances and xenobiotics that they are exposed to in new environments. Here, we explore the flavin-containing monooxygenases (FMOs) that are essential for xenobiotic detoxification. Employing a paleobiochemistry approach in combination with enzymology techniques we disclose the set of historical substitutions responsible for the family's functional diversification in tetrapods. Remarkably, a few amino acid replacements differentiate an ancestral multi-tasking FMO into a more specialized monooxygenase by modulating the oxygenating flavin intermediate. Our findings substantiate an ongoing premise that enzymatic function hinges on a subset of residues that is not limited to the active site core.

Rossmann-toolbox: a deep learning-based protocol for the prediction and design of cofactor specificity in Rossmann fold proteins

The Rossmann fold enzymes are involved in essential biochemical pathways such as nucleotide and amino acid metabolism. Their functioning relies on interaction with cofactors, small nucleoside-based compounds specifically recognized by a conserved βαβ motif shared by all Rossmann fold proteins. While Rossmann methyltransferases recognize only a single cofactor type, the S-adenosylmethionine, the oxidoreductases, depending on the family, bind nicotinamide (nicotinamide adenine dinucleotide, nicotinamide adenine dinucleotide phosphate) or flavin-based (flavin adenine dinucleotide) cofactors. In this study, we showed that despite its short length, the βαβ motif unambiguously defines the specificity towards the cofactor. Following this observation, we trained two complementary deep learning models for the prediction of the cofactor specificity based on the sequence and structural features of the βαβ motif. A benchmark on two independent test sets, one containing βαβ motifs bearing no resemblance to those of the training set, and the other comprising 38 experimentally confirmed cases of rational design of the cofactor specificity, revealed the nearly perfect performance of the two methods. The Rossmann-toolbox protocols can be accessed via the webserver at https://lbs.cent.uw.edu.pl/rossmann-toolbox and are available as a Python package at https://github.com/labstructbioinf/rossmann-toolbox.

The adaptive landscape of a metallo-enzyme is shaped by environment-dependent epistasis

Enzymes can evolve new catalytic activity when environmental changes present them with novel substrates. Despite this seemingly straightforward relationship, factors other than the direct catalytic target can also impact adaptation. Here, we characterize the catalytic activity of a recently evolved bacterial methyl-parathion hydrolase for all possible combinations of the five functionally relevant mutations under eight different laboratory conditions (in which an alternative divalent metal is supplemented). The resultant adaptive landscapes across this historical evolutionary transition vary in terms of both the number of “fitness peaks” as well as the genotype(s) at which they are found as a result of genotype-by-environment interactions and environment-dependent epistasis. This suggests that adaptive landscapes may be fluid and molecular adaptation is highly contingent not only on obvious factors (such as catalytic targets), but also on less obvious secondary environmental factors that can direct it towards distinct outcomes.

The metaphor of an adaptive landscape is presented quantitatively by looking at molecular adaptations and their catalytic consequences in a recently evolved bacterial enzyme. The study identifies both genotype-by-environment interactions and environment-dependent epistasis as factors that can alter the fitness of functional mutations.

Overview of structurally homologous flavoprotein oxidoreductases containing the low Mr thioredoxin reductase-like fold - A functionally diverse group

Structural studies show that enzymes have a limited number of unique folds, although structurally related enzymes have evolved to perform a large variety of functions. In this review, we have focused on enzymes containing the low molecular weight thioredoxin reductase (low M_r TrxR) fold. This fold consists of two domains, both containing a three-layer ββα sandwich Rossmann-like fold, serving as flavin adenine dinucleotide (FAD) and, in most cases, pyridine nucleotide (NAD(P)H) binding-domains. Based on a search of the Protein Data Bank for all published structures containing the low M_r TrxR-like fold, we here present a comprehensive overview of enzymes with this structural architecture. These range from TrxR-like ferredoxin/flavodoxin NAD(P)⁺ oxidoreductases, through glutathione reductase, to NADH peroxidase. Some enzymes are solely composed of the low M_r TrxR-like fold, while others contain one or two additional domains. In this review, we give a detailed description of selected enzymes containing only the low M_r TrxR-like fold, however, catalyzing a diversity of chemical reactions. Our overview of this structurally similar, yet functionally distinct group of flavoprotein oxidoreductases highlights the fascinating and increasing number of studies describing the diversity among these enzymes, especially during the last decade(s).

Effect of His6-tag Position on the Expression and Properties of Phenylacetone Monooxygenase from Thermobifida fusca

Phenylacetone monooxygenase (EC 1.14.13.92, PAMО) catalyzes oxidation of ketones with molecular oxygen and NADPH with the formation of esters. PAMО is a promising enzyme for biotechnological processes. In this work, we generated genetic constructs coding for PAMO from Thermobifida fusca, containing N- or C-terminal His₆-tags (PAMO N and PAMO C, respectively), as well as PAMO L with the His₆-tag attached to the enzyme C-terminus via a 19-a.a. spacer. All PAMO variants were expressed as catalytically active proteins in Escherichia coli BL21(DE3) cells; however, the expression level of PAMO N was 3 to 5 times higher than for the other two enzymes. The catalytic constants (k_cat) of PAMO C and PAMO L were similar to that published for PAMO L produced in a different expression system; the catalytic constant for PAMO N was slightly lower (by 15%). The values of Michaelis constants with NADPH for all PAMО variants were in agreement within the published data for PAMO L (within the experimental error); however, the K_M for benzylacetone was several times higher. Thermal inactivation studies and differential scanning calorimetry demonstrated that the thermal stability of PAMO N was 3 to 4 times higher compared to that of the enzymes with the C-terminal His₆-tag.

Protein Engineering for Nicotinamide Coenzyme Specificity in Oxidoreductases: Attempts and Challenges

Oxidoreductases are ubiquitous enzymes that catalyze an extensive range of chemical reactions with great specificity, efficiency, and selectivity. Most oxidoreductases are nicotinamide cofactor-dependent enzymes with a strong preference for NADP or NAD. Because these coenzymes differ in stability, bioavailability and costs, the enzyme preference for a specific coenzyme is an important issue for practical applications. Different approaches for the manipulation of coenzyme specificity have been reported, with different degrees of success. Here we present various attempts for the switching of nicotinamide coenzyme preference in oxidoreductases by protein engineering. This review covers 103 enzyme engineering studies from 82 articles and evaluates the accomplishments in terms of coenzyme specificity and catalytic efficiency compared to wild type enzymes of different classes. We analyzed different protein engineering strategies and related them with the degree of success in inverting the cofactor specificity. In general, catalytic activity is compromised when coenzyme specificity is reversed, however when switching from NAD to NADP, better results are obtained. In most of the cases, rational strategies were used, predominantly with loop exchange generating the best results. In general, the tendency of removing acidic residues and incorporating basic residues is the strategy of choice when trying to change specificity from NAD to NADP, and vice versa. Computational strategies and algorithms are also covered as helpful tools to guide protein engineering strategies. This mini review aims to give a general introduction to the topic, giving an overview of tools and information to work in protein engineering for the reversal of coenzyme specificity.

Identification of structural determinants of NAD(P)H selectivity and lysine binding in lysine N(6)-monooxygenase

l-lysine (l-Lys) N(6)-monooxygenase (NbtG), from Nocardia farcinica, is a flavin-dependent enzyme that catalyzes the hydroxylation of l-Lys in the presence of oxygen and NAD(P)H in the biosynthetic pathway of the siderophore nocobactin. NbtG displays only a 3-fold preference for NADPH over NADH, different from well-characterized related enzymes, which are highly selective for NADPH. The structure of NbtG with bound NAD(P)(+) or l-Lys is currently not available. Herein, we present a mutagenesis study targeting M239, R301, and E216. These amino acids are conserved and located in either the NAD(P)H binding domain or the l-Lys binding pocket. M239R resulted in high production of hydrogen peroxide and little hydroxylation with no change in coenzyme selectivity. R301A caused a 300-fold decrease on kcat/Km value with NADPH but no change with NADH. E216Q increased the Km value for l-Lys by 30-fold with very little change on the kcat value or in the binding of NAD(P)H. These results suggest that R301 plays a major role in NADPH selectivity by interacting with the 2'-phosphate of the adenine-ribose moiety of NADPH, while E216 plays a role in l-Lys binding.

Exploring nicotinamide cofactor promiscuity in NAD(P)H-dependent flavin containing monooxygenases (FMOs) using natural variation within the phosphate binding loop. Structure and activity of FMOs from Cellvibrio sp. BR and Pseudomonas stutzeri NF13

•

Genes encoding ‘Type II FMOs’ CFMO and PSFMO were selected for diversity in the NADPH binding loop.

•

CFMO and PSFMO were shown to both accept either NADH or NADPH as cofactor in the reduction of FAD.

•

The activities with both NADPH and NADH have been evaluated in the oxidation of sulfide substrates.

•

Structures of CFMO and PSFMO were solved, revealing the nature of the NADPH phosphate binding loop.

Flavin-containing monooxygenases (FMOs) catalyse asymmetric oxidation reactions that have potential for preparative organic synthesis, but most use the more expensive, phosphorylated nicotinamide cofactor NADPH to reduce FAD to FADH₂ prior to formation of the (hydro)peroxy intermediate required for substrate oxygenation. A comparison of the structures of NADPH-dependent FMO from Methylophaga aminisulfidivorans (mFMO) and SMFMO from Stenotrophomonas maltophilia, which is able to use both NADPH and NADH, suggested that the promiscuity of the latter enzyme may be due in part to the substitution of an Arg–Thr couple in the NADPH phosphate recognition site in mFMO, for a Gln–His couple in SMFMO (Jensen et al., 2012, Chembiochem, 13, 872–878). Natural variation within the phosphate binding region, and its influence on nicotinamide cofactor promiscuity, was explored through the cloning, expression, characterisation and structural studies of FMOs from Cellvibrio sp. BR (CFMO) and Pseudomonas stutzeri NF13 (PSFMO), which possess Thr–Ser and Gln–Glu in the putative phosphate recognition positions, respectively. CFMO and PSFMO displayed 5- and 1.5-fold greater activity, respectively, than SMFMO for the reduction of FAD with NADH, and were also cofactor promiscuous, displaying a ratio of activity with NADH:NADPH of 1.7:1 and 1:1.3, respectively. The structures of CFMO and PSFMO revealed the context of the phosphate binding loop in each case, and also clarified the structure of the mobile helix–loop–helix motif that appears to shield the FAD-binding pocket from bulk solvent in this class of FMOs, a feature that was absent from the structure of SMFMO.