O2 reactivity of flavoproteins: dynamic access of dioxygen to the active site and role of a H+ relay system in D-amino acid oxidase

Development of growth selection system and pocket engineering of d-amino acid oxidase to enhance selective deamination activity toward d-phosphinothricin

D-amino acid oxidase (DAAO)-catalyzed selective oxidative deamination is a very promising process for synthesizing l-amino acids including l-phosphinothricin ( l-PPT, a high-efficiency and broad-spectrum herbicide). However, the wild-type DAAO's low activity toward unnatural substrates like d-phosphinothricin ( d-PPT) hampers its application. Herein, a DAAO from Caenorhabditis elegans (CeDAAO) was screened and engineered to improve the catalytic potential on d-PPT. First, we designed a novel growth selection system, taking into account the intricate relationship between the growth of Escherichia coli (E. coli) and the catalytic mechanism of DAAO. The developed system was used for high-throughput screening of gene libraries, resulting in the discovery of a variant (M6) with significantly increased catalytic activity against d-PPT. The variant displays different catalytic properties on substrates with varying hydrophobicity and hydrophilicity. Analysis using Alphafold2 modeling and molecular dynamic simulations showed that the reason for the enhanced activity was the substrate-binding pocket with enlarged size and suitable charge distribution. Further QM/MM calculations revealed that the crucial factor for enhancing activity lies in reducing the initial energy barrier of the reductive half reaction. Finally, a comprehensive binding-model index to predict the enhanced activity of DAAO toward d-PPT, and an enzymatic deracemization approach was developed, enabling the efficient synthesis of l-PPT with remarkable efficiency.

Directed evolution unlocks oxygen reactivity for a nicotine-degrading flavoenzyme

The flavoenzyme nicotine oxidoreductase (NicA2) is a promising injectable treatment to aid in the cessation of smoking, a behavior responsible for one in ten deaths worldwide. NicA2 acts by degrading nicotine in the bloodstream before it reaches the brain. Clinical use of NicA2 is limited by its poor catalytic activity in the absence of its natural electron acceptor CycN. Without CycN, NicA2 is instead oxidized slowly by dioxygen (O2), necessitating unfeasibly large doses in a therapeutic setting. Here, we report a genetic selection strategy that directly links CycN-independent activity of NicA2 to growth of Pseudomonas putida S16. This selection enabled us to evolve NicA2 variants with substantial improvement in their rate of oxidation by O2. The encoded mutations cluster around a putative O2 tunnel, increasing flexibility and accessibility to O2 in this region. These mutations further confer desirable clinical properties. A variant form of NicA2 is tenfold more effective than the wild type at degrading nicotine in the bloodstream of rats.

D-Amino Acids and Cancer: Friends or Foes?

α-amino acids exist in two configurations, named D-(dextro) and L-(levo) enantiomers. L-amino acids are used in protein synthesis and play a central role in cell metabolism. The effects of the L-amino acid composition of foods and the dietary modifications of this composition on the efficacy of cancer therapies have been widely investigated in relation to the growth and reproduction of cancerous cells. However, less is known about the involvement of D-amino acids. In recent decades, D-amino acids have been identified as natural biomolecules that play interesting and specific roles as common components of the human diet. Here, we focus on recent investigations showing altered D-amino acid levels in specific cancer types and on the various roles proposed for these biomolecules related to cancer cell proliferation, cell protection during therapy, and as putative, innovative biomarkers. Notwithstanding recent progress, the relationship between the presence of D-amino acids, their nutritional value, and cancer cell proliferation and survival represents an underrated scientific issue. Few studies on human samples have been reported to date, suggesting a need for routine analysis of D-amino acid content and an evaluation of the enzymes involved in regulating their levels in clinical samples in the near future.

Development of a Versatile Method to Construct Direct Electron Transfer-Type Enzyme Complexes Employing SpyCatcher/SpyTag System

The electrochemical enzyme sensors based on direct electron transfer (DET)-type oxidoreductase-based enzymes are ideal for continuous and in vivo monitoring. However, the number and types of DET-type oxidoreductases are limited. The aim of this research is the development of a versatile method to create a DET-type oxidoreductase complex based on the SpyCatcher/SpyTag technique by preparing SpyCatcher-fused heme c and SpyTag-fused non-DET-type oxidoreductases, and by the in vitro formation of DET-type oxidoreductase complexes. A heme c containing an electron transfer protein derived from Rhizobium radiobacter (CYTc) was selected to prepare SpyCatcher-fused heme c. Three non-DET-type oxidoreductases were selected as candidates for the SpyTag-fused enzyme: fungi-derived flavin adenine dinucleotide (FAD)-dependent glucose dehydrogenase (GDH), an engineered FAD-dependent d-amino acid oxidase (DAAOx), and an engineered FMN-dependent l-lactate oxidase (LOx). CYTc-SpyCatcher (CYTc-SC) and SpyTag-Enzymes (ST-GDH, ST-DAAOx, ST-LOx) were prepared as soluble molecules while maintaining their redox properties and catalytic activities, respectively. CYTc-SC/ST-Enzyme complexes were formed by mixing CYTc-SpyCatcher and SpyTag-Enzymes, and the complexes retained their original enzymatic activity. Remarkably, the heme domain served as an electron acceptor from complexed enzymes by intramolecular electron transfer; consequently, all constructed CYTc-SC/ST-Enzyme complexes showed DET ability to the electrode, demonstrating the versatility of this method.

Transient potentiometry based d-serine sensor using engineered d-amino acid oxidase showing quasi-direct electron transfer property

d-Serine biosensing has been extensively reported based on enzyme sensors using flavin adenine dinucleotide (FAD) -dependent d-amino acid oxidase (DAAOx), based on the monitoring of hydrogen peroxide generated by the enzymatic reaction, which is affected by dissolved oxygen concentration in the measurement environment in in vivo use. Here we report a novel sensing principle for d-serine, transient potentiometry based d-serine sensor using engineered DAAOx showing quasi-direct electron transfer (DET) property. DAAOx Gly52Val mutant, revealed to possess dye-mediated dehydrogenase activity using artificial synthetic electron acceptors, while its oxidase activity was negligible. The enzyme was immobilized on electrode and was modified with amine-reactive phenazine ethosulfate, resulted an enzyme electrode showing quasi-DET type response. Although OCP based monitoring took more than several minutes to obtain steady state OCP value, the time dependent OCP change monitoring, transient potentiometry, provided rapid and sensitive sensor signals. While dOCP/dt based monitoring was suitable for sensing with longer than 5 s time resolution with d-serine concentration range between 0.5 mM and 5 mM, dOCP/d t based monitoring is suitable for d-serine monitoring with much shorter time resolution (less than 1 s) with high sensitivity with wider dynamic range (20 μM-30 mM). The maximum dOCP/d t was -39.2 ± 2.0 mV/s^1/2, the K_m(app) was 1.9 mM, and the lower limit of detection was 20 μM. In addition, d-serine monitoring was also possible in the artificial cerebrospinal fluid. The transient potentiometry based sensing reported in this study will be further utilized to realize miniaturized, continuous, real-time, in vivo sensor for d-serine monitoring.

Reactive oxygen species as a double-edged sword: The role of oxidative enzymes in antitumor therapy

A number of approaches have been developed over the years to manage cancer, such as chemotherapy using low-molecular-mass molecules and radiotherapy. Here, enzymes can also find useful applications. Among them, oxidases have attracted attention because of their ability to produce reactive oxygen species (ROS, especially hydrogen peroxide) in tumors and potentially modulate the production of this cytotoxic compound when enzymes active on substrates present in low amounts are used, such as the d-amino acid oxidase and d-amino acid couple system. These treatments have been also developed for additional cancer treatment approaches, such as phototherapy, nutrient starvation, and metal-induced hydroxyl radical production. In addition, to improve tumor specificity and decrease undesired side effects, oxidases have been targeted by means of nanotechnologies and protein engineering (i.e., by designing chimeric proteins able to accumulate in the tumor). The most recent advances obtained by using six different oxidases (i.e., the FAD-containing enzymes glucose oxidase, d- and l-amino acid oxidases, cholesterol oxidase and xanthine oxidase, and the copper-containing amine oxidase) have been reported. Anticancer therapy based on oxidase-based ROS production has now reached maturity and can be applied in the clinic.

Catalytic mechanism of ancestral L-lysine oxidase assigned by sequence data mining

A large number of protein sequences are registered in public databases such as PubMed. Functionally uncharacterized enzymes are included in these databases, some of which likely have potential for industrial applications. However, assignment of the enzymes remained difficult tasks for now. In this study, we assigned a total of 28 original sequences to uncharacterized enzymes in the FAD-dependent oxidase family expressed in some species of bacteria including Chryseobacterium, Flavobacterium, and Pedobactor. Progenitor sequence of the assigned 28 sequences was generated by ancestral sequence reconstruction, and the generated sequence exhibited L-lysine oxidase activity; thus, we named the enzyme AncLLysO. Crystal structures of ligand-free and ligand-bound forms of AncLLysO were determined, indicating that the enzyme recognizes L-Lys by hydrogen bond formation with R76 and E383. The binding of L-Lys to AncLLysO induced dynamic structural change at a plug loop formed by residues 251 to 254. Biochemical assays of AncLLysO variants revealed the functional importance of these substrate recognition residues and the plug loop. R76A and E383D variants were also observed to lose their activity, and the kcat/Km value of G251P and Y253A mutations were approximately 800- to 1800-fold lower than that of AncLLysO, despite the indirect interaction of the substrates with the mutated residues. Taken together, our data demonstrate that combinational approaches to sequence classification from database and ancestral sequence reconstruction may be effective not only to find new enzymes using databases of unknown sequences but also to elucidate their functions.

Human D-aspartate Oxidase: A Key Player in D-aspartate Metabolism

In recent years, the D-enantiomers of amino acids have been recognized as natural molecules present in all kingdoms, playing a variety of biological roles. In humans, d-serine and d-aspartate attracted attention for their presence in the central nervous system. Here, we focus on d-aspartate, which is involved in glutamatergic neurotransmission and the synthesis of various hormones. The biosynthesis of d-aspartate is still obscure, while its degradation is due to the peroxisomal flavin adenine dinucleotide (FAD)-containing enzyme d-aspartate oxidase. d-Aspartate emergence is strictly controlled: levels decrease in brain within the first days of life while increasing in endocrine glands postnatally and through adulthood. The human d-aspartate oxidase (hDASPO) belongs to the d-amino acid oxidase-like family: its tertiary structure closely resembles that of human d-amino acid oxidase (hDAAO), the enzyme that degrades neutral and basic d-amino acids. The structure-function relationships of the physiological isoform of hDASPO (named hDASPO_341) and the regulation of gene expression and distribution and properties of the longer isoform hDASPO_369 have all been recently elucidated. Beyond the substrate preference, hDASPO and hDAAO also differ in kinetic efficiency, FAD-binding affinity, pH profile, and oligomeric state. Such differences suggest that evolution diverged to create two different ways to modulate d-aspartate and d-serine levels in the human brain. Current knowledge about hDASPO is shedding light on the molecular mechanisms underlying the modulation of d-aspartate levels in human tissues and is pushing novel, targeted therapeutic strategies. Now, it has been proposed that dysfunction in NMDA receptor-mediated neurotransmission is caused by disrupted d-aspartate metabolism in the nervous system during the onset of various disorders (such as schizophrenia): the design of suitable hDASPO inhibitors aimed at increasing d-aspartate levels thus represents a novel and useful form of therapy.

The family of sarcosine oxidases: Same reaction, different products

The subfamily of sarcosine oxidase is a set of enzymes within the larger family of amine oxidases. It is ubiquitously distributed among different kingdoms of life. The member enzymes catalyze the oxidization of an N-methyl amine bond of amino acids to yield unstable imine species that undergo subsequent spontaneous non-enzymatic reactions, forming an array of different products. These products range from demethylated simple species to complex alkaloids. The enzymes belonging to the sarcosine oxidase family, namely, monomeric and heterotetrameric sarcosine oxidase, l-pipecolate oxidase, N-methyltryptophan oxidase, NikD, l-proline dehydrogenase, FsqB, fructosamine oxidase and saccharopine oxidase have unique features differentiating them from other amine oxidases. This review highlights the key attributes of the sarcosine oxidase family enzymes, in terms of their substrate binding motif, type of oxidation reaction mediated and FAD regeneration, to define the boundaries of this group and demarcate these enzymes from other amine oxidase families.

Transient pockets as mediators of gas molecules routes inside proteins: The case study of dioxygen pathway in homogentisate 1,2-dioxygenase and its implication in Alkaptonuria development

Alkaptonuria (AKU) is an ultra-rare disease caused by mutations in homogentisate 1,2-dioxygenase (HGD) enzyme, characterized by the loss of enzymatic activity and the accumulation of its substrate, homogentisic acid (HGA) in different tissues, leading to ochronosis and organ degeneration. Although the pathological effects of HGD mutations are largely studied, less is known about the structure of the enzyme, in particular the pathways for dioxygen diffusion to the active site, required for the enzymatic reaction, are still uninvestigated. In the present project, the combination of two in silico techniques, Molecular Dynamics (MD) simulation and Implicit Ligand Sampling (ILS), was used to delineate gas diffusion routes in HGD enzyme. A route from the central opening of the hexameric structure of the enzyme to the back of the active site trough the protein moiety was identified as the path for dioxygen diffusion, also overlapping with a transient pocket, which then assumes an important role in dioxygen diffusion. Along the route the sequence location of the missense variant E401Q, responsible for AKU development, was also found, suggesting such mutation to be conducive of enzymatic activity loss by altering the flow dynamics of dioxygen. Our in silico approach allowed also to delineate the route of HGA substrate to the active site, until now only supposed.

Biosensors for D-Amino Acids: Detection Methods and Applications

D-enantiomers of amino acids (D-AAs) are only present in low amounts in nature, frequently at trace levels, and for this reason, their biological function was undervalued for a long time. In the past 25 years, the improvements in analytical methods, such as gas chromatography, HPLC, and capillary electrophoresis, allowed to detect D-AAs in foodstuffs and biological samples and to attribute them specific biological functions in mammals. These methods are time-consuming, expensive, and not suitable for online application; however, life science investigations and industrial applications require rapid and selective determination of D-AAs, as only biosensors can offer. In the present review, we provide a status update concerning biosensors for detecting and quantifying D-AAs and their applications for safety and quality of foods, human health, and neurological research. The review reports the main challenges in the field, such as selectivity, in order to distinguish the different D-AAs present in a solution, the simultaneous assay of both L- and D-AAs, the production of implantable devices, and surface-scanning biosensors. These innovative tools will push future research aimed at investigating the neurological role of D-AAs, a vibrant field that is growing at an accelerating pace.

Alteration of Electron Acceptor Preferences in the Oxidative Half-Reaction of Flavin-Dependent Oxidases and Dehydrogenases

In this review, recent progress in the engineering of the oxidative half-reaction of flavin-dependent oxidases and dehydrogenases is discussed, considering their current and future applications in bioelectrochemical studies, such as for the development of biosensors and biofuel cells. There have been two approaches in the studies of oxidative half-reaction: engineering of the oxidative half-reaction with oxygen, and engineering of the preference for artificial electron acceptors. The challenges for engineering oxidative half-reactions with oxygen are further categorized into the following approaches: (1) mutation to the putative residues that compose the cavity where oxygen may be located, (2) investigation of the vicinities where the reaction with oxygen may take place, and (3) investigation of possible oxygen access routes to the isoalloxazine ring. Among these approaches, introducing a mutation at the oxygen access route to the isoalloxazine ring represents the most versatile and effective strategy. Studies to engineer the preference of artificial electron acceptors are categorized into three different approaches: (1) engineering of the charge at the residues around the substrate entrance, (2) engineering of a cavity in the vicinity of flavin, and (3) decreasing the glycosylation degree of enzymes. Among these approaches, altering the charge in the vicinity where the electron acceptor may be accessed will be most relevant.

A detailed mechanism of the oxidative half-reaction of d-amino acid oxidase: another route for flavin oxidation

d-Amino acid oxidase (DAAO) is a flavoenzyme whose inhibition is expected to have therapeutic potential in schizophrenia. DAAO catalyses hydride transfer from the substrate to the flavin in the reductive half-reaction, and the flavin is reoxidized by O2 in the oxidative half-reaction. Quantum mechanical/molecular mechanical calculations were performed and their results together with available experimental information were used to elucidate the detailed mechanism of the oxidative half-reaction. The reaction starts with a single electron transfer from FAD to O2, followed by triplet-singlet transition. FAD oxidation is completed by a proton coupled electron transfer to the oxygen species and the reaction terminates with H2O2 formation by proton transfer from the oxidized substrate to the oxygen species via a chain of water molecules. The substrate plays a double role by facilitating the first electron transfer and by providing a proton in the last step. The mechanism differs from the oxidative half-reaction of other oxidases.

Following the Evolutionary Track of a Highly Specific l-Arginine Oxidase by Reconstruction and Biochemical Analysis of Ancestral and Native Enzymes

Following the evolutionary track of enzymes can help elucidate how enzymes attain their characteristic functions, such as thermal adaptation and substrate selectivity, during the evolutionary process. Ancestral sequence reconstruction (ASR) is effective for following evolutionary processes if sufficient sequence data are available. Selecting sequences from the data to generate a curated sequence library is necessary for the successful design of artificial proteins by ASR. In this study, we tried to follow the evolutionary track of l-arginine oxidase (AROD), a flavin adenine dinucleotide (FAD)-dependent amino acid oxidase (LAAO) that exhibits high specificity for l-arginine. The library was generated by selecting sequences in which the 15th, 50th, 332nd, and 580th residues are Gly, Ser, Trp, and Thr, respectively. We excluded sequences that are either extremely short or long and those with a low degree of sequence identity. Three ancestral ARODs (AncARODn0, AncARODn1, and AncARODn2) were designed using the library. Subsequently, we expressed the ancestral ARODs as well as native Oceanobacter kriegii AROD (OkAROD) in bacteria. AncARODn0 is phylogenetically most remote from OkAROD, whereas AncARODn2 is most similar to OkAROD. Thermal stability was gradually increased by extending AROD sequences back to the progenitor, while the temperature at which the residual activity is half of the maximum measured activity (T _1/2) of AncARODn0 was >20°C higher than that of OkAROD. Remarkably, only AncARODn0 exhibited broad substrate selectivity similar to that of conventional promiscuous LAAO. Taken together, our findings led us to infer that AROD may have evolved from a highly thermostable and promiscuous LAAO.IMPORTANCE In this study, we attempted to infer the molecular evolution of a recently isolated FAD-dependent l-arginine oxidase (AROD) that oxidizes l-arginine to 2-ketoarginine. Utilizing 10 candidate AROD sequences, we obtained a total of three ancestral ARODs. In addition, one native AROD was obtained by cloning one of the candidate ARODs. The candidate sequences were selected utilizing a curation method defined in this study. All the ARODs were successfully expressed in Escherichia coli for analysis of their biochemical functions. The catalytic activity of our bacterially expressed ancestral ARODs suggests that our ASR was successful. The ancestral AROD that is phylogenetically most remote from a native AROD has the highest thermal stability and substrate promiscuity. Our findings led us to infer that AROD evolved from a highly thermostable and promiscuous LAAO. As an application, we can design artificial ARODs with improved functions compared with those of native ones.

Oxygen Activation Switch in the Copper Amine Oxidase of Escherichia coli

Copper amine oxidases (CuAOs) are metalloenzymes that reduce molecular oxygen to hydrogen peroxide during catalytic turnover of primary amines. In addition to Cu²⁺ in the active site, two peripheral calcium sites, ∼32 Å from the active site, have roles in Escherichia coli amine oxidase (ECAO). The buried Ca²⁺ (Asp533, Leu534, Asp535, Asp678, and Ala679) is essential for full-length protein production, while the surface Ca²⁺ (Glu573, Tyr667, Asp670, and Glu672) modulates biogenesis of the 2,4,5-trihydroxyphenylalanine quinone (TPQ) cofactor. The E573Q mutation at the surface site prevents calcium binding and TPQ biogenesis. However, TPQ biogenesis can be restored by a suppressor mutation (I342F) in the proposed oxygen delivery channel to the active site. While supporting TPQ biogenesis (∼60% WTECAO TPQ), I342F/E573Q has almost no amine oxidase activity (∼4.6% WTECAO activity). To understand how these long-range mutations have major effects on TPQ biogenesis and catalysis, we employed ultraviolet-visible spectroscopy, steady-state kinetics, inhibition assays, and X-ray crystallography. We show that the surface metal site controls the equilibrium (disproportionation) of the Cu²⁺-substrate reduced TPQ (TPQ_AMQ) Cu⁺-TPQ semiquinone (TPQ_SQ) couple. Removal of the calcium ion from this site by chelation or mutagenesis shifts the equilibrium to Cu²⁺-TPQ_AMQ or destabilizes Cu⁺-TPQ_SQ. Crystal structure analysis shows that TPQ biogenesis is stalled at deprotonation in the Cu²⁺-tyrosinate state. Our findings support WTECAO using the inner sphere electron transfer mechanism for oxygen reduction during catalysis, and while a Cu⁺-tyrosyl radical intermediate is not essential for TPQ biogenesis, it is required for efficient biogenesis.

Modulating D-amino acid oxidase (DAAO) substrate specificity through facilitated solvent access

D-amino acid oxidase (DAAO) degrades D-amino acids to produce α-ketoacids, hydrogen peroxide and ammonia. DAAO has often been investigated and engineered for industrial and clinical applications. We combined information from literature with a detailed analysis of the structure to engineer mammalian DAAOs. The structural analysis was complemented with molecular dynamics simulations to characterize solvent accessibility and product release mechanisms. We identified non-obvious residues located on the loops on the border between the active site and the secondary binding pocket essential for pig and human DAAO substrate specificity and activity. We engineered DAAOs by mutating such critical residues and characterised the biochemical activity of the resulting variants. The results highlight the importance of the selected residues in modulating substrate specificity, product egress and enzyme activity, suggesting further steps of DAAO re-engineering towards desired clinical and industrial applications.

The ins and outs of vanillyl alcohol oxidase: Identification of ligand migration paths

Vanillyl alcohol oxidase (VAO) is a homo-octameric flavoenzyme belonging to the VAO/PCMH family. Each VAO subunit consists of two domains, the FAD-binding and the cap domain. VAO catalyses, among other reactions, the two-step conversion of p-creosol (2-methoxy-4-methylphenol) to vanillin (4-hydroxy-3-methoxybenzaldehyde). To elucidate how different ligands enter and exit the secluded active site, Monte Carlo based simulations have been performed. One entry/exit path via the subunit interface and two additional exit paths have been identified for phenolic ligands, all leading to the si side of FAD. We argue that the entry/exit path is the most probable route for these ligands. A fourth path leading to the re side of FAD has been found for the co-ligands dioxygen and hydrogen peroxide. Based on binding energies and on the behaviour of ligands in these four paths, we propose a sequence of events for ligand and co-ligand migration during catalysis. We have also identified two residues, His466 and Tyr503, which could act as concierges of the active site for phenolic ligands, as well as two other residues, Tyr51 and Tyr408, which could act as a gateway to the re side of FAD for dioxygen. Most of the residues in the four paths are also present in VAO’s closest relatives, eugenol oxidase and p-cresol methylhydroxylase. Key path residues show movements in our simulations that correspond well to conformations observed in crystal structures of these enzymes. Preservation of other path residues can be linked to the electron acceptor specificity and oligomerisation state of the three enzymes. This study is the first comprehensive overview of ligand and co-ligand migration in a member of the VAO/PCMH family, and provides a proof of concept for the use of an unbiased method to sample this process.

Enzymes are bionanomachines, which speed up chemical reactions in organisms. To understand how they achieve that, we need to study their mechanisms. Computational enzymology can show us what happens in the enzyme’s active site during a reaction. But molecules need first to reach the active site before a reaction can start. The process of substrate entry and product exit to the active site is often neglected when studying enzymes. However, these two events are of fundamental importance to the proper functioning of any enzyme. We are interested in these dynamic processes to complete our understanding of the mode of action of enzymes. In our work, we have studied substrate and product migration in vanillyl alcohol oxidase. This enzyme can produce the flavour vanillin and enantiopure alcohols, but also catalyses other reactions. The named products are of interest to the flavour- and fine-chemical industries.

On Atomistic Models for Molecular Oxygen

Molecular oxygen (O₂) is key to all life on earth, as it is constantly cycled via photosynthesis and cellular respiration. Substantial scientific effort has been devoted to understanding every part of this cycle. Classical molecular dynamics (MD) simulations have been used to study some of the key processes involved in cellular respiration: O₂ permeation through alveolar monolayers and cellular membranes, its binding to hemoglobin during transport in the bloodstream, as well as its transport along optimal pathways toward its reduction sites in proteins. Moreover, MD simulations can help interpret the results of several imaging techniques in which O₂ is used because of its paramagnetic nature. However, despite the widespread use of computational models for the O₂ molecule, their performances have never been systematically evaluated. In this paper, we assess the performances of 14 different models of O₂ available in the literature by calculating four thermodynamic properties: density, heat of vaporization, free energy of hydration, and free energy of solvation in hexadecane. For each property, reliable experimental data are available. Most models perform reasonably well in predicting the correct trends, but they fail to reproduce the experimental data quantitatively. We then develop new models for O₂, with and without a quadrupole moment, and compare their behavior with the behavior of previously published models. The new models show significant improvement in terms of density, heat of vaporization, and free energy of hydration. However, quantitative agreement with water-oil partitioning is not reached due to discrepancies between the calculated and measured free energies of solvation in hexadecane. We suggest that classical pairwise-additive models may be inadequate to properly describe the thermodynamics of solvation of apolar species, such as O₂, in apolar solvents.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

Control of the Position of Oxygen Delivery in Soybean Lipoxygenase-1 by Amino Acid Side Chains within a Gas Migration Channel

Understanding gas migration pathways is critical to unraveling structure-function relationships in enzymes that process gaseous substrates such as O2, H2, and N2 This work investigates the role of a defined pathway for O2 in regulating the peroxidation of linoleic acid by soybean lipoxygenase 1. Computational and mutagenesis studies provide strong support for a dominant delivery channel that shuttles molecular oxygen to a specific region of the active site, thereby ensuring the regio- and stereospecificity of product. Analysis of reaction kinetics and product distribution in channel mutants also reveals a plasticity to the gas migration pathway. The findings show that a single site mutation (I553W) limits oxygen accessibility to the active site, greatly increasing the fraction of substrate that reacts with oxygen free in solution. They also show how a neighboring site mutation (L496W) can result in a redirection of oxygen toward an alternate position of the substrate, changing the regio- and stereospecificity of peroxidation. The present data indicate that modest changes in a protein scaffold may modulate the access of small gaseous molecules to enzyme-bound substrates.

Structure-Function Relationships in l-Amino Acid Deaminase, a Flavoprotein Belonging to a Novel Class of Biotechnologically Relevant Enzymes

l-Amino acid deaminase from Proteus myxofaciens (PmaLAAD) is a membrane flavoenzyme that catalyzes the deamination of neutral and aromatic l-amino acids into α-keto acids and ammonia. PmaLAAD does not use dioxygen to re-oxidize reduced FADH2 and thus does not produce hydrogen peroxide; instead, it uses a cytochrome b-like protein as an electron acceptor. Although the overall fold of this enzyme resembles that of known amine or amino acid oxidases, it shows the following specific structural features: an additional novel α+β subdomain placed close to the putative transmembrane α-helix and to the active-site entrance; an FAD isoalloxazine ring exposed to solvent; and a large and accessible active site suitable to bind large hydrophobic substrates. In addition, PmaLAAD requires substrate-induced conformational changes of part of the active site, particularly in Arg-316 and Phe-318, to achieve the correct geometry for catalysis. These studies are expected to pave the way for rationally improving the versatility of this flavoenzyme, which is critical for biocatalysis of enantiomerically pure amino acids.

Role of F357 as an Oxygen Gate in the Oxidative Half-Reaction of Choline Oxidase

Choline oxidase from Arthrobacter globiformis catalyzes the oxidation of choline to glycine betaine by using oxygen as an electron acceptor. A partially rate limiting isomerization of the reduced wild-type enzyme during the reaction with oxygen was previously detected using solvent viscosity effects. In this study, we hypothesized that the side chains of M62 and F357, located at the entrance to the active site of choline oxidase, may be related to the slow isomerization detected. We engineered a double-variant enzyme M62A/F357A. The kinetic characterization of the double-variant enzyme showed a lack of the isomerization detected in wild-type choline oxidase, and a lack of saturation with an oxygen concentration as high as 1 mM, while most other kinetic parameters were similar to those of wild-type choline oxidase. The kinetic characterization of the single-variant enzymes established that only the side chain of F357 plays a role in the isomerization of choline oxidase in the oxidative half-reaction. Molecular dynamics studies suggest that the slow isomerization related to F357 is possibly due to the participation of the phenyl ring in a newly proposed gating mechanism for a narrow tunnel, assumed to regulate the access of oxygen to the reduced cofactor.

All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

The N-terminal pre-A region of Mycobacterium tuberculosis 2/2HbN promotes NO-dioxygenase activity

A unique defense mechanisms by which Mycobacterium tuberculosis protects itself from nitrosative stress is based on the O2 -dependent NO-dioxygenase (NOD) activity of truncated hemoglobin 2/2HbN (Mt2/2HbN). The NOD activity largely depends on the efficiency of ligand migration to the heme cavity through a two-tunnel (long and short) system; recently, it was also correlated with the presence at the Mt2/2HbN N-terminus of a short pre-A region, not conserved in most 2/2HbNs, whose deletion results in a drastic reduction of NO scavenging. In the present study, we report the crystal structure of Mt2/2HbN-ΔpreA, lacking the pre-A region, at a resolution of 1.53 Å. We show that removal of the pre-A region results in long range effects on the protein C-terminus, promoting the assembly of a stable dimer, both in the crystals and in solution. In the Mt2/2HbN-ΔpreA dimer, access of heme ligands to the short tunnel is hindered. Molecular dynamics simulations show that the long tunnel branch is the only accessible pathway for O2 -ligand migration to/from the heme, and that the gating residue Phe(62)E15 partly restricts the diameter of the tunnel. Accordingly, kinetic measurements indicate that the kon value for peroxynitrite isomerization by Mt2/2HbN-ΔpreA-Fe(III) is four-fold lower relative to the full-length protein, and that NO scavenging by Mt2/2HbN-ΔpreA-Fe(II)-O2 is reduced by 35-fold. Therefore, we speculate that Mt2/2HbN evolved to host the pre-A region as a mechanism for preventing dimerization, thus reinforcing the survival of the microorganism against the reactive nitrosative stress in macrophages.

DATABASE

Coordinates and structure factors have been deposited in the Protein Data Bank under accession number 5AB8.

Oxygen diffusion pathways in a cofactor-independent dioxygenase

A combination of experimental and computational studies reveals the main O₂ diffusion pathways, providing insight into how cofactor-independent oxygenases control stereospecificity and prevent oxidative inactivation.

Molecular oxygen plays an important role in a wide variety of enzymatic reactions. Through recent research efforts combining computational and experimental methods a new view of O₂ diffusion is emerging, where specific channels guide O₂ to the active site. The focus of this work is DpgC, a cofactor-independent oxygenase. Molecular dynamics simulations, together with mutagenesis experiments and xenon-binding data, reveal that O₂ reaches the active site of this enzyme using three main pathways and four different access points. These pathways connect a series of dynamic hydrophobic pockets, concentrating O₂ at a specific face of the enzyme substrate. Extensive molecular dynamics simulations provide information about which pathways are more frequently used. This data is consistent with the results of kinetic measurements on mutants and is difficult to obtain using computational cavity-location methods. Taken together, our results reveal that although DpgC is rare in its ability of activating O₂ in the absence of cofactors or metals, the way O₂ reaches the active site is similar to that reported for other O₂-using proteins: multiple access channels are available, and the architecture of the pathway network can provide regio- and stereoselectivity. Our results point to the existence of common themes in O₂ access that are conserved among very different types of proteins.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O₂) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O₂ to travel from the solvent/membrane to the heme a₃-Cu_B binuclear center (BNC). Moreover, the identification of all possible O₂ channels as well as the atomic details of O₂ diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O₂ distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O₂ channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O₂ molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O₂ inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the BNC. Both alternative O₂ channels start in the membrane spanning region and terminate close to Y288_I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O₂ to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O₂ to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O₂ delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O₂ can diffuse via three well-defined channels that start in membrane region (where O₂ solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the binuclear center.

Engineering pyranose 2-oxidase for modified oxygen reactivity

Pyranose 2-oxidase (POx), a member of the GMC family of flavoproteins, catalyzes the regioselective oxidation of aldopyranoses at position C2 to the corresponding 2-ketoaldoses. During the first half-reaction, FAD is reduced to FADH2 and reoxidized in the second half-reaction by reducing molecular oxygen to H2O2. Alternative electron acceptors including quinones, radicals or chelated metal ions show significant and in some cases even higher activity. While oxygen as cheap and abundantly available electron acceptor is favored for many processes, reduced oxygen reactivity is desirable for some applications such as in biosensors/biofuel cells because of reduced oxidative damages to the biocatalyst from concomitant H2O2 production as well as reduced electron "leakage" to oxygen. The reactivity of flavoproteins with oxygen is of considerable scientific interest, and the determinants of oxygen activation and reactivity are the subject of numerous studies. We applied site-saturation mutagenesis on a set of eleven amino acids around the active site based on the crystal structure of the enzyme. Using microtiter plate screening assays with peroxidase/2,2'-azino-bis(3-ethylbenzthiazoline-6-sulphonic acid) and 2,6-dichlorophenolindophenol, variants of POx with decreased oxidase activity and maintained dehydrogenase activity were identified. Variants T166R, Q448H, L545C, L547R and N593C were characterized with respect to their apparent steady-state constants with oxygen and the alternative electron acceptors DCPIP, 1,4-benzoquinone and ferricenium ion, and the effect of the mutations was rationalized based on structural properties.

Aldonolactone oxidoreductases

Vitamin C is a widely used vitamin. Here we review the occurrence and properties of aldonolactone oxidoreductases, an important group of flavoenzymes responsible for the ultimate production of vitamin C and its analogs in animals, plants, and single-cell organisms.

L-amino acid oxidase as biocatalyst: a dream too far?

L-amino acid oxidase (LAAO) is a flavoenzyme containing non-covalently bound flavin adenine dinucleotide, which catalyzes the stereospecific oxidative deamination of l-amino acids to α-keto acids and also produces ammonia and hydrogen peroxide via an imino acid intermediate. LAAOs purified from snake venoms are the best-studied members of this family of enzymes, although a number of LAAOs from bacterial and fungal sources have been also reported. From a biochemical point of view, LAAOs from different sources are distinguished by molecular mass, substrate specificity, post-translational modifications and regulation. In analogy to the well-known biotechnological applications of d-amino acid oxidase, important results are expected from the availability of suitable LAAOs; however, these expectations have not been fulfilled yet because none of the "true" LAAOs has successfully been expressed as a recombinant protein in prokaryotic hosts, such as Escherichia coli. In enzyme biotechnology, recombinant production of a protein is mandatory both for the production of large amounts of the catalyst and to improve its biochemical properties by protein engineering. As an alternative, flavoenzymes active on specific l-amino acids have been identified, e.g., l-aspartate oxidase, l-lysine oxidase, l-phenylalanine oxidase, etc. According to presently available information, amino acid oxidases with "narrow" or "strict" substrate specificity represent as good candidates to obtain an enzyme more suitable for biotechnological applications by enlarging their substrate specificity by means of protein engineering.

Engineering of pyranose dehydrogenase for increased oxygen reactivity

Pyranose dehydrogenase (PDH), a member of the GMC family of flavoproteins, shows a very broad sugar substrate specificity but is limited to a narrow range of electron acceptors and reacts extremely slowly with dioxygen as acceptor. The use of substituted quinones or (organo)metals as electron acceptors is undesirable for many production processes, especially of food ingredients. To improve the oxygen reactivity, site-saturation mutagenesis libraries of twelve amino acids around the active site of Agaricus meleagris PDH were expressed in Saccharomyces cerevisiae. We established high-throughput screening assays for oxygen reactivity and standard dehydrogenase activity using an indirect Amplex Red/horseradish peroxidase and a DCIP/D-glucose based approach. The low number of active clones confirmed the catalytic role of H512 and H556. Only one position was found to display increased oxygen reactivity. Histidine 103, carrying the covalently linked FAD cofactor in the wild-type, was substituted by tyrosine, phenylalanine, tryptophan and methionine. Variant H103Y was produced in Pichia pastoris and characterized and revealed a five-fold increase of the oxygen reactivity.

Molecular dynamics simulations reveal highly permeable oxygen exit channels shared with water uptake channels in photosystem II

Photosystem II (PSII) catalyzes the oxidation of water in the conversion of light energy into chemical energy in photosynthesis. Water delivery and oxygen removal from the oxygen evolving complex (OEC), buried deep within PSII, are critical requirements to facilitate the reaction and minimize reactive oxygen damage. It has often been assumed that water and oxygen travel through separate channels within PSII, as demonstrated in cytochrome c oxidase. This study describes all-atom molecular dynamics simulations of PSII designed to investigate channels by fully characterizing the distribution and permeation of both water and oxygen. Interestingly, most channels found in PSII were permeable to both oxygen and water, however individual channels exhibited different energetic barriers for the two solutes. Several routes for oxygen diffusion within PSII with low energy permeation barriers were found, ensuring its fast removal from the OEC. In contrast, all routes for water showed significant energy barriers, corresponding to a much slower permeation rate for water through PSII. Two major factors were responsible for this selectivity: (1) hydrogen bonds between water and channel amino acids, and (2) steric restraints. Our results reveal the presence of a shared network of channels in PSII optimized to both facilitate the quick removal of oxygen and effectively restrict the water supply to the OEC to help stabilize and protect it from small water soluble inhibitors.

Monoamine oxidases in development

Monoamine oxidases (MAOs) are flavoproteins of the outer mitochondrial membrane that catalyze the oxidative deamination of biogenic and xenobiotic amines. In mammals there are two isoforms (MAO-A and MAO-B) that can be distinguished on the basis of their substrate specificity and their sensitivity towards specific inhibitors. Both isoforms are expressed in most tissues, but their expression in the central nervous system and their ability to metabolize monoaminergic neurotransmitters have focused MAO research on the functionality of the mature brain. MAO activities have been related to neurodegenerative diseases as well as to neurological and psychiatric disorders. More recently evidence has been accumulating indicating that MAO isoforms are expressed not only in adult mammals, but also before birth, and that defective MAO expression induces developmental abnormalities in particular of the brain. This review is aimed at summarizing and critically evaluating the new findings on the developmental functions of MAO isoforms during embryogenesis.

Fluorescent protein barrel fluctuations and oxygen diffusion pathways in mCherry

Fluorescent proteins (FPs) are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. In this paper, we use molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variant. We employ implicit ligand sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. We also show that these pathways can be blocked or altered and barrel fluctuations can be reduced by strategic amino acid substitutions.

Flavin-linked Erv-family sulfhydryl oxidases release superoxide anion during catalytic turnover

Typically, simple flavoprotein oxidases couple the oxidation of their substrates with the formation of hydrogen peroxide without release of significant levels of the superoxide ion. However, two evolutionarily related single-domain sulfhydryl oxidases (Erv2p; a yeast endoplasmic reticulum resident protein and augmenter of liver regeneration, ALR, an enzyme predominantly found in the mitochondrial intermembrane) release up to ~30% of the oxygen they reduce as the superoxide ion. Both enzymes oxidize dithiol substrates via a redox-active disulfide adjacent to the flavin cofactor within the helix-rich Erv domain. Subsequent reduction of the flavin is followed by transfer of reducing equivalents to molecular oxygen. Superoxide release was initially detected using tris(3-hydroxypropyl)phosphine (THP) as an alternative reducing substrate to dithiothreitol (DTT). THP, and other phosphines, showed anomalously high turnover numbers with Erv2p and ALR in the oxygen electrode, but oxygen consumption was drastically suppressed upon the addition of superoxide dismutase. The superoxide ion initiates a radical chain reaction promoting the aerobic oxidation of phosphines with the formation of hydrogen peroxide. Use of a known flux of superoxide generated by the xanthine/xanthine oxidase system showed that one superoxide ion stimulates the reduction of 27 and 4.5 molecules of oxygen using THP and tris(2-carboxyethyl)phosphine (TCEP), respectively. This superoxide-dependent amplification of oxygen consumption by phosphines provides a new kinetic method for the detection of superoxide. Superoxide release was also observed by a standard chemiluminescence method using a luciferin analogue (MCLA) when 2 mM DTT was employed as a substrate of Erv2p and ALR. The percentage of superoxide released from Erv2p increased to ~65% when monomeric mutants of the normally homodimeric enzyme were used. In contrast, monomeric multidomain quiescin sulfhydryl oxidase enzymes that also contain an Erv FAD-binding fold release only 1-5% of their total reduced oxygen species as the superoxide ion. Aspects of the mechanism and possible physiological significance of superoxide release from these Erv-domain flavoproteins are discussed.

Oxygen reactivity in flavoenzymes: context matters

Many flavoenzymes--oxidases and monooxygenases--react faster with oxygen than free flavins do. There are many ideas on how enzymes cause this. Recent work has focused on the importance of a positive charge near N5 of the reduced flavin. Fructosamine oxidase has a lysine near N5 of its flavin. We measured a rate constant of 1.6 × 10(5) M(-1) s(-1) for its reaction with oxygen. The Lys276Met mutant reacted with a rate constant of 291 M(-1) s(-1), suggesting an important role for this lysine in oxygen activation. The dihydroorotate dehydrogenases from E. coli and L. lactis also have a lysine near N5 of the flavin. They react with O(2) with rate constants of 6.2 × 10(4) and 3.0 × 10(3) M(-1) s(-1), respectively. The Lys66Met and Lys43Met mutant enzymes react with rate constants that are nearly the same as those for the wild-type enzymes, demonstrating that simply placing a positive charge near N5 of the flavin does not guarantee increased oxygen reactivity. Our results show that the lysine near N5 does not exert an effect without an appropriate context; evolution did not find only one mechanism for activating the reaction of flavins with O(2).

Modulating O2 reactivity in a fungal flavoenzyme: involvement of aryl-alcohol oxidase Phe-501 contiguous to catalytic histidine

Aryl-alcohol oxidase (AAO) is a flavoenzyme responsible for activation of O(2) to H(2)O(2) in fungal degradation of lignin. The AAO crystal structure shows a buried active site connected to the solvent by a hydrophobic funnel-shaped channel, with Phe-501 and two other aromatic residues forming a narrow bottleneck that prevents the direct access of alcohol substrates. However, ligand diffusion simulations show O(2) access to the active site following this channel. Site-directed mutagenesis of Phe-501 yielded a F501A variant with strongly reduced O(2) reactivity. However, a variant with increased reactivity, as shown by kinetic constants and steady-state oxidation degree, was obtained by substitution of Phe-501 with tryptophan. The high oxygen catalytic efficiency of F501W, ∼2-fold that of native AAO and ∼120-fold that of F501A, seems related to a higher O(2) availability because the turnover number was slightly decreased with respect to the native enzyme. Free diffusion simulations of O(2) inside the active-site cavity of AAO (and several in silico Phe-501 variants) yielded >60% O(2) population at 3-4 Å from flavin C4a in F501W compared with 44% in AAO and only 14% in F501A. Paradoxically, the O(2) reactivity of AAO decreased when the access channel was enlarged and increased when it was constricted by introducing a tryptophan residue. This is because the side chain of Phe-501, contiguous to the catalytic histidine (His-502 in AAO), helps to position O(2) at an adequate distance from flavin C4a (and His-502 Nε). Phe-501 substitution with a bulkier tryptophan residue resulted in an increase in the O(2) reactivity of this flavoenzyme.

Molecular mimicry and ligand recognition in binding and catalysis by the histone demethylase LSD1-CoREST complex

Histone demethylases LSD1 and LSD2 (KDM1A/B) catalyze the oxidative demethylation of Lys4 of histone H3. We used molecular dynamics simulations to probe the diffusion of the oxygen substrate. Oxygen can reach the catalytic center independently from the presence of a bound histone peptide, implying that LSD1 can complete subsequent demethylation cycles without detaching from the nucleosomal particle. The simulations highlight the role of a strictly conserved active-site Lys residue providing general insight into the enzymatic mechanism of oxygen-reacting flavoenzymes. The crystal structure of LSD1-CoREST bound to a peptide of the transcription factor SNAIL1 unravels a fascinating example of molecular mimicry. The SNAIL1 N-terminal residues bind to the enzyme active-site cleft, effectively mimicking the H3 tail. This finding predicts that other members of the SNAIL/Scratch transcription factor family might associate to LSD1/2. The combination of selective histone-modifying activity with the distinct recognition mechanisms underlies the biological complexity of LSD1/2.

Pleiotropic impact of a single lysine mutation on biosynthesis of and catalysis by N-methyltryptophan oxidase

N-Methyltryptophan oxidase (MTOX) contains covalently bound FAD. N-Methyltryptophan binds in a cavity above the re face of the flavin ring. Lys259 is located above the opposite, si face. Replacement of Lys259 with Gln, Ala, or Met blocks (>95%) covalent flavin incorporation in vivo. The mutant apoproteins can be reconstituted with FAD. Apparent turnover rates (k(cat,app)) of the reconstituted enzymes are ~2500-fold slower than those of wild-type MTOX. Wild-type MTOX forms a charge-transfer E(ox)·S complex with the redox-active anionic form of NMT. The E(ox)·S complex formed with Lys259Gln does not exhibit a charge-transfer band and is converted to a reduced enzyme·imine complex (EH(2)·P) at a rate 60-fold slower than that of wild-type MTOX. The mutant EH(2)·P complex contains the imine zwitterion and exhibits a charge-transfer band, a feature not observed with the wild-type EH(2)·P complex. Reaction of reduced Lys259Gln with oxygen is 2500-fold slower than that of reduced wild-type MTOX. The latter reaction is unaffected by the presence of bound product. Dissociation of the wild-type EH(2)·P complex is 80-fold slower than k(cat). The mutant EH(2)·P complex dissociates 15-fold faster than k(cat,app). Consequently, EH(2)·P and free EH(2) are the species that react with oxygen during turnover of the wild-type and mutant enzyme, respectively. The results show that (i) Lys259 is the site of oxygen activation in MTOX and also plays a role in holoenzyme biosynthesis and N-methyltryptophan oxidation and (ii) MTOX contains separate active sites for N-methyltryptophan oxidation and oxygen reduction on opposite faces of the flavin ring.

New biotech applications from evolved D-amino acid oxidases

D-Amino acid oxidase (DAAO) is a well-known flavoenzyme that catalyzes the oxygen-dependent oxidative deamination of amino acid D-isomers with absolute stereospecificity, which results in α-keto acids, ammonia and hydrogen peroxide. Recently, the extraordinary functional plasticity of DAAO has become evident; in turn, boosting research on this flavoprotein. Protein engineering has allowed for a redesign of DAAO substrate specificity, oxygen affinity, cofactor binding, stability, and oligomeric state. We review recent developments in utilizing DAAO, including as a biocatalyst for resolving racemic amino acid mixtures, as a tool for biosensing, and as a new mechanism of herbicide resistance. Perspectives for future biotechnological applications of this oxidative biocatalyst are also outlined.