Beyond active site residues: overall structural dynamics control catalysis in flavin-containing and heme-containing monooxygenases

Evolution of the catalytic mechanism at the dawn of the Baeyer-Villiger monooxygenases

Enzymes are crucial for the emergence and sustenance of life on earth. How they became catalytically active during their evolution is still an open question. Two opposite explanations are plausible: acquiring a mechanism in a series of discrete steps or all at once in a single evolutionary event. Here, we use molecular phylogeny, ancestral sequence reconstruction, and biochemical characterization to follow the evolution of a specialized group of flavoprotein monooxygenases, the bacterial Baeyer-Villiger monooxygenases (BVMOs). These enzymes catalyze an intricate chemical reaction relying on three different elements: a reduced nicotinamide cofactor, dioxygen, and a substrate. Characterization of ancestral BVMOs shows that the catalytic mechanism evolved in a series of steps starting from a FAD-binding protein and further acquiring reactivity and specificity toward each of the elements participating in the reaction. Together, the results of our work portray how an intrinsically complex catalytic mechanism emerged during evolution.

A Conserved Allosteric Site on Drug-Metabolizing CYPs: A Systematic Computational Assessment

Cytochrome P450 enzymes (CYPs) are the largest group of enzymes involved in human drug metabolism. Ligand tunnels connect their active site buried at the core of the membrane-anchored protein to the surrounding solvent environment. Recently, evidence of a superficial allosteric site, here denoted as hotspot 1 (H1), involved in the regulation of ligand access in a soluble prokaryotic CYP emerged. Here, we applied multi-scale computational modeling techniques to study the conservation and functionality of this allosteric site in the nine most relevant mammalian CYPs responsible for approximately 70% of drug metabolism. In total, we systematically analyzed over 44 μs of trajectories from conventional MD, cosolvent MD, and metadynamics simulations. Our bioinformatic analysis and simulations with organic probe molecules revealed the site to be well conserved in the CYP2 family with the exception of CYP2E1. In the presence of a ligand bound to the H1 site, we could observe an enlargement of a ligand tunnel in several members of the CYP2 family. Further, we could detect the facilitation of ligand translocation by H1 interactions with statistical significance in CYP2C8 and CYP2D6, even though all other enzymes except for CYP2C19, CYP2E1, and CYP3A4 presented a similar trend. As the detailed comprehension of ligand access and egress phenomena remains one of the most relevant challenges in the field, this work contributes to its elucidation and ultimately helps in estimating the selectivity of metabolic transformations using computational techniques.

Characterization of a Nitro-Forming Enzyme Involved in Fosfazinomycin Biosynthesis

N-hydroxylating monooxygenases (NMOs) are a subclass of flavin-dependent enzymes that hydroxylate nitrogen atoms. Recently, unique NMOs that perform multiple reactions on one substrate molecule have been identified. Fosfazinomycin M (FzmM) is one such NMO, forming nitrosuccinate from aspartate (Asp) in the fosfazinomycin biosynthetic pathway in some Streptomyces sp. This work details the biochemical and kinetic analysis of FzmM. Steady-state kinetic investigation shows that FzmM performs a coupled reaction with Asp (k_cat, 3.0 ± 0.01 s^-1) forming nitrosuccinate, which can be converted to fumarate and nitrite by the action of FzmL. FzmM displays a 70-fold higher k_cat/K_M value for NADPH compared to NADH and has a narrow optimal pH range (7.5-8.0). Contrary to other NMOs where the k_red is rate-limiting, FzmM exhibits a very fast k_red (50 ± 0.01 s^-1 at 4 °C) with NADPH. NADPH binds at a K_D value of ∼400 μM, and hydride transfer occurs with pro-R stereochemistry. Oxidation of FzmM in the absence of Asp exhibits a spectrum with a shoulder at ∼370 nm, consistent with the formation of a C(4a)-hydroperoxyflavin intermediate, which decays into oxidized flavin and hydrogen peroxide at a rate 100-fold slower than the k_cat. This reaction is enhanced in the presence of Asp with a slightly faster k_ox than the k_cat, suggesting that flavin dehydration or Asp oxidation is partially rate limiting. Multiple sequence analyses of FzmM to NMOs identified conserved residues involved in flavin binding but not for NADPH. Additional sequence analysis to related monooxygenases suggests that FzmM shares sequence motifs absent in other NMOs.

Mycobacterial crypto-AcpM as a tool to investigate the consequence of drug binding on its key FAS II partner enzyme HadAB

Background Mycobacterial FASII pathway is governed by the Protein-Protein Interaction mediated dynamics existent between Acyl Carrier Protein and its partner enzymes. The dehydratase HadAB, involved in the third step of FASII synthesis has remained a key target of drugs like Thiacetazone (TAC) and its consequence on AcpM binding is yet to be deciphered. Owing to the transient nature of these interactions, analysing their implications as a drug target has been exhausting. Methods In this context, we have developed an in vitro method to study the effect of thiocarbamide-containing compounds, TAC and SPA0355 (a thiourea analogue) against mycobacterial HadAB. Additionally, by utilizing crypto-ACP (NBD-tagged Acyl Carrier Protein) as a tool of our choice, we attempted at exploring the effect of TAC and SPA0355 on mycobacterial HadAB. Results SPA0355 behaves at par with TAC and undergoes activation in the presence of monooxygenase EthA thus, bringing about a covalent modification in HadA subunit of HadAB. The crypto-ACP method provides insights into the altered substrate housing capability in HadAB associated with the impediment of its AcpM mediated functionality; an outcome attributed to the repercussions associated with the binding of the aforementioned thiourea compounds. Conclusion This investigation has assisted in unveiling a two-step mechanism undertaken by AcpM for interacting with its corresponding partner protein during acyl chain transfer. General significance This study highlights the alterations brought about by drug binding in the interplay between ACP and HadAB. Additionally, this work for the first time establishes the role of SPA0355 as a promising drug candidate against dehydratase HadAB.

Flavoprotein monooxygenases: Versatile biocatalysts

Flavoprotein monooxygenases (FPMOs) are single- or two-component enzymes that catalyze a diverse set of chemo-, regio- and enantioselective oxyfunctionalization reactions. In this review, we describe how FPMOs have evolved from model enzymes in mechanistic flavoprotein research to biotechnologically relevant catalysts that can be applied for the sustainable production of valuable chemicals. After a historical account of the development of the FPMO field, we explain the FPMO classification system, which is primarily based on protein structural properties and electron donor specificities. We then summarize the most appealing reactions catalyzed by each group with a focus on the different types of oxygenation chemistries. Wherever relevant, we report engineering strategies that have been used to improve the robustness and applicability of FPMOs.

Computational design of enzymes for biotechnological applications

Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.

New frontiers in flavin-dependent monooxygenases

Flavin-dependent monooxygenases catalyze a wide variety of redox reactions in important biological processes and are responsible for the synthesis of highly complex natural products. Although much has been learned about FMO chemistry in the last ~80 years of research, several aspects of the reactions catalyzed by these enzymes remain unknown. In this review, we summarize recent advancements in the flavin-dependent monooxygenase field including aspects of flavin dynamics, formation and stabilization of reactive species, and the hydroxylation mechanism. Novel catalysis of flavin-dependent N-oxidases involving consecutive oxidations of amines to generate oximes or nitrones is presented and the biological relevance of the products is discussed. In addition, the activity of some FMOs have been shown to be essential for the virulence of several human pathogens. We also discuss the biomedical relevance of FMOs in antibiotic resistance and the efforts to identify inhibitors against some members of this important and growing family enzymes.

The devil is in the details: The chemical basis and mechanistic versatility of flavoprotein monooxygenases

The ubiquitous flavoenzymes commonly catalyze redox chemistry such as the monooxygenation of organic substrates and are both widely utilized in nature (e.g., in primary and secondary metabolism) and of significant industrial interest. In this work, we highlight the structural and mechanistic characteristics of the distinct types of flavoprotein monooxygenases (FPMOs). We thereby illustrate the chemical basis of FPMO catalysis, which enables reactions such as (aromatic) hydroxylation, epoxidation, (de)halogenation, heteroatom oxygenation, Baeyer-Villiger oxidation, α-hydroxylation of ketones, or non-oxidative carbon-hetero bond cleavage. This seemingly unmatched versatility in oxygenation chemistry results from extensive fine-tuning and regiospecific functionalization of the flavin cofactor that is tightly controlled by the surrounding protein matrix. Accordingly, FPMOs steer the formation of covalent flavin-oxygen adducts for oxygen transfer in the form of the classical flavin-C4a-(hydro)peroxide or the recently discovered N5-functionalized flavins (i.e. the flavin-N5-oxide and the flavin-N5-peroxide), while in rare cases covalent oxygen adduct formation may be foregone entirely. Finally, we speculate about hitherto undiscovered flavin-mediated oxygenation reactions and compare FPMOs to cytochrome P450 monooxygenases, before addressing open questions and challenges for the future investigation of FPMOs.

Riboflavin Is Directly Involved in N-Dealkylation Catalyzed by Bacterial Cytochrome P450 Monooxygenases

Like a vast number of enzymes in nature, bacterial cytochrome P450 monooxygenases require an activated form of flavin as a cofactor for catalytic activity. Riboflavin is the precursor of FAD and FMN that serves as indispensable cofactor for flavoenzymes. In contrast to previous notions, herein we describe the identification of an electron-transfer process that is directly mediated by riboflavin for N-dealkylation by bacterial P450 monooxygenases. The electron relay from NADPH to riboflavin and then via activated oxygen to heme was proposed based on a combination of X-ray crystallography, molecular modeling and molecular dynamics simulation, site-directed mutagenesis and biochemical analysis of representative bacterial P450 monooxygenases. This study provides new insights into the electron transfer mechanism in bacterial P450 enzyme catalysis and likely in yeasts, fungi, plants and mammals.

Using in vivo multiphoton fluorescence lifetime imaging to unravel disease-specific changes in the liver redox state

Multiphoton fluorescence lifetime microscopy has revolutionized studies of pathophysiological and xenobiotic dynamics, enabling the spatial and temporal quantification of these processes in intact organs in vivo. We have previously used multiphoton fluorescence lifetime microscopy to characterise the morphology and amplitude weighted mean fluorescence lifetime of the endogenous fluorescent metabolic cofactor nicotinamide adenine dinucleotide (phosphate) (NAD(P)H) of mouse livers in vivo following induction of various disease states. Here, we extend the characterisation of liver disease models by using nonlinear regression to estimate the unbound, bound fluorescence lifetimes for NAD(P)H, flavin adenine dinucleotide (FAD), along with metabolic ratios and examine the impact of using multiple segmentation methods. We found that NAD(P)H amplitude ratio, and fluorescence lifetime redox ratio can be used as discriminators of diseased liver from normal liver. The redox ratio provided a sensitive measure of the changes in hepatic fibrosis and biliary fibrosis. Hepatocellular carcinoma was associated with an increase in spatial heterogeneity and redox ratio coupled with a decrease in mean fluorescence lifetime. We conclude that multiphoton fluorescence lifetime microscopy parameters and metabolic ratios provided insights into the in vivo redox state of diseased compared to normal liver that were not apparent from a global, mean fluorescence lifetime measurement alone.

Insights into the Mechanism of Ethionamide Resistance in Mycobacterium tuberculosis through an in silico Structural Evaluation of EthA and Mutants Identified in Clinical Isolates

Mutation in the ethionamide (ETH) activating enzyme, EthA, is the main factor determining resistance to this drug, used to treat TB patients infected with MDR and XDR Mycobacterium tuberculosis isolates. Many mutations in EthA of ETH resistant (ETH-R) isolates have been described but their roles in resistance remain uncharacterized, partly because structural studies on the enzyme are lacking. Thus, we took a two-tier approach to evaluate two mutations (Y50C and T453I) found in ETH-R clinical isolates. First, we used a combination of comparative modeling, molecular docking, and molecular dynamics to build an EthA model in complex with ETH that has hallmark features of structurally characterized homologs. Second, we used free energy computational calculations for the reliable prediction of relative free energies between the wild type and mutant enzymes. The ΔΔG values for Y50C and T453I mutant enzymes in complex with FADH2-NADP-ETH were 3.34 (+/−0.55) and 8.11 (+/−0.51) kcal/mol, respectively, compared to the wild type complex. The positive ΔΔG values indicate that the wild type complex is more stable than the mutants, with the T453I complex being the least stable. These are the first results shedding light on the molecular basis of ETH resistance, namely reduced complex stability of mutant EthA.

Natural Variation in the ‘Control Loop’ of BVMOAFL210 and Its Influence on Regioselectivity and Sulfoxidation

Baeyer-Villiger monooxygenases (BVMOs) are flavin-dependent enzymes that primarily convert ketones to esters, but can also catalyze heteroatom oxidation. Several structural studies have highlighted the importance of the ‘control loop’ in BVMOs, which adopts different conformations during catalysis. Central to the ‘control loop’ is a conserved tryptophan that has been implicated in NADP(H) binding. BVMOAFL210 from Aspergillus flavus, however, contains a threonine in the equivalent position. Here, we report the structure of BVMOAFL210 in complex with NADP+ in both the ‘open’ and ‘closed’ conformations. In neither conformation does Thr513 contact the NADP+. Although mutagenesis of Thr513 did not significantly alter the substrate scope, changes in peroxyflavin stability and reaction rates were observed. Mutation of this position also brought about changes in the regio- and enantioselectivity of the enzyme. Moreover, lower rates of overoxidation during sulfoxidation of thioanisole were also observed.

Ancestral-sequence reconstruction unveils the structural basis of function in mammalian FMOs

Flavin-containing monooxygenases (FMOs) are ubiquitous in all domains of life and metabolize a myriad of xenobiotics, including toxins, pesticides and drugs. However, despite their pharmacological importance, structural information remains bereft. To further our understanding behind their biochemistry and diversity, we used ancestral-sequence reconstruction, kinetic and crystallographic techniques to scrutinize three ancient mammalian FMOs: AncFMO2, AncFMO3-6 and AncFMO5. Remarkably, all AncFMOs could be crystallized and were structurally resolved between 2.7- and 3.2-Å resolution. These crystal structures depict the unprecedented topology of mammalian FMOs. Each employs extensive membrane-binding features and intricate substrate-profiling tunnel networks through a conspicuous membrane-adhering insertion. Furthermore, a glutamate-histidine switch is speculated to induce the distinctive Baeyer-Villiger oxidation activity of FMO5. The AncFMOs exhibited catalysis akin to human FMOs and, with sequence identities between 82% and 92%, represent excellent models. Our study demonstrates the power of ancestral-sequence reconstruction as a strategy for the crystallization of proteins.

Stabilization of cyclohexanone monooxygenase by computational and experimental library design

Enzymes often by far exceed the activity, selectivity, and sustainability achieved with chemical catalysts. One of the main reasons for the lack of biocatalysis in the chemical industry is the poor stability exhibited by many enzymes when exposed to process conditions. This dilemma is exemplified in the usually very temperature‐sensitive enzymes catalyzing the Baeyer–Villiger reaction, which display excellent stereo‐ and regioselectivity and offer a green alternative to the commonly used, explosive peracids. Here we describe a protein engineering approach applied to cyclohexanone monooxygenase from Rhodococcus sp. HI‐31, a substrate‐promiscuous enzyme that efficiently catalyzes the production of the nylon‐6 precursor ε‐caprolactone. We used a framework for rapid enzyme stabilization by computational libraries (FRESCO), which predicts protein‐stabilizing mutations. From 128 screened point mutants, approximately half had a stabilizing effect, albeit mostly to a small degree. To overcome incompatibility effects observed upon combining the best hits, an easy shuffled library design strategy was devised. The most stable and highly active mutant displayed an increase in unfolding temperature of 13°C and an approximately 33x increase in half‐life at 30°C. In contrast to the wild‐type enzyme, this thermostable 8x mutant is an attractive biocatalyst for biotechnological applications.