Unveiling the Pathways of Dioxygen Through the C2 Component of the Environmentally Relevant Monooxygenase p‐Hydroxyphenylacetate Hydroxylase from Acinetobacter baumannii: A Molecular Dynamics Investigation

In this work, models of the homotetrameric C2 component of the monooxygenase p‐hydroxyphenylacetate hydroxylase from Acinetobacter baumannii, in complex with dioxygen (O₂) and, or not, the substrate p‐hydroxyphenylacetate (HPA) were built. Both models proved to be amenable to random‐acceleration molecular dynamics (RAMD) simulations, whereby a tiny randomly oriented external force, acting on O₂ at the active site in front of flavin mononucleotide (FMNH^?), accelerated displacement of O₂ toward the bulk solvent. This allowed us to carry out a sufficiently large number of RAMD simulations to be of statistical significance. The two systems behaved very similarly under RAMD, except for O₂ leaving the active site more easily in the absence of HPA, but then finding similar obstacles in getting to the gate as when the active site was sheltered by HPA. This challenges previous conclusions that HPA can only reach the active center after that the C4aOOH derivative of FMNH^? is formed, requiring uptake of O₂ at the active site before HPA. According to these RAMD simulations, O₂ could well get to FMNH^? also in the presence of the substrate at the active site.     

Structural and biochemical characterization of EDTA monooxygenase and its physical interaction with a partner flavin reductase

Ethylenediaminetetraacetate (EDTA) is currently the most abundant organic pollutant due to its recalcitrance and extensive use. Only a few bacteria can degrade it, using EDTA monooxygenase (EmoA) to initiate the degradation. EmoA is an FMNH₂‐dependent monooxygenase that requires an NADH:FMN oxidoreductase (EmoB) to provide FMNH₂ as a cosubstrate. Although EmoA has been identified from Chelativorans (ex. Mesorhizobium) sp. BNC1, its catalytic mechanism is unknown. Crystal structures of EmoA revealed a domain‐like insertion into a TIM‐barrel, which might serve as a flexible lid for the active site. Docking of MgEDTA²⁻ into EmoA identified an intricate hydrogen bond network connected to Tyr⁷¹, which should potentially lower its pKa. Tyr⁷¹, along with nearby Glu⁷⁰ and a peroxy flavin, facilitates a keto‐enol transition of the leaving acetyl group of EDTA. Further, for the first time, the physical interaction between EmoA and EmoB was observed by ITC, molecular docking and enzyme kinetic assay, which enhanced both EmoA and EmoB activities probably through coupled channelling of FMNH₂.     

Purification and Characterization of EDTA Monooxygenase from the EDTA-Degrading Bacterium BNC1

The synthetic chelating agent EDTA can mobilize radionuclides and heavy metals in the environment. Biodegradation of EDTA should reduce this mobilization. Although several bacteria have been reported to mineralize EDTA, little is known about the biochemistry of EDTA degradation. Understanding the biochemistry will facilitate the removal of EDTA from the environment. EDTA-degrading activities were detected in cell extracts of bacterium BNC1 when flavin mononucleotide (FMN), NADH, and O₂ were present. The degradative enzyme system was separated into two different enzymes, EDTA monooxygenase and an FMN reductase. EDTA monooxygenase oxidized EDTA to glyoxylate and ethylenediaminetriacetate (ED3A), with the coconsumption of FMNH₂ and O₂. The FMN reductase provided EDTA monooxygenase with FMNH₂ by reducing FMN with NADH. The FMN reductase was successfully substituted in the assay mixture by other FMN reductases. EDTA monooxygenase was purified to greater than 95% homogeneity and had a single polypeptide with a molecular weight of 45,000. The enzyme oxidized both EDTA complexed with various metal ions and uncomplexed EDTA. The optimal conditions for activity were pH 7.8 and 35°C. K_ms were 34.1 μM for uncomplexed EDTA and 8.5 μM for MgEDTA²⁻; this difference in K_m indicates that the enzyme has greater affinity for MgEDTA²⁻. The enzyme also catalyzed the release of glyoxylate from nitrilotriacetate and diethylenetriaminepentaacetate. EDTA monooxygenase belongs to a small group of FMNH₂-utilizing monooxygenases that attack carbon-nitrogen, carbon-sulfur, and carbon-carbon double bonds.     

Catalytic role of a conserved cysteine residue in the desulfonation reaction by the alkanesulfonate monooxygenase enzyme

Detailed kinetic studies were performed in order to determine the role of the single cysteine residue in the desulfonation reaction catalyzed by SsuD. Mutation of the conserved cysteine at position 54 in SsuD to either serine or alanine had little effect on FMNH₂ binding. The k_cat/K_m value for the C54S SsuD variant increased 3-fold, whereas the k_cat/K_m value for C54A SsuD decreased 6-fold relative to wild-type SsuD. An initial fast phase was observed in kinetic traces obtained for the oxidation of flavin at 370 nm when FMNH₂ was mixed against C54S SsuD (k_obs, 111 s^− 1) in oxygenated buffer that was 10-fold faster than wild-type SsuD (k_obs, 12.9 s^− 1). However, there was no initial fast phase observed in similar kinetic traces obtained for C54A SsuD. This initial fast phase was previously assigned to the formation of the C4a-(hydro)peroxyflavin in studies with wild-type SsuD. There was no evidence for the formation of the C4a-(hydro)peroxyflavin with either SsuD variant when octanesulfonate was included in rapid reaction kinetic studies, even at low octanesulfonate concentrations. The absence of any C4a-(hydro)peroxyflavin accumulation correlates with the increased catalytic activity of C54S SsuD. These results suggest that the conservative serine substitution is able to effectively take the place of cysteine in catalysis. Conversely, decreased accumulation of the C4a-(hydro)peroxyflavin intermediate with the C54A SsuD variant may be due to decreased activity. The data described suggest that Cys54 in SsuD may be either directly or indirectly involved in stabilizing the C4a-(hydro)peroxyflavin intermediate formed during catalysis through hydrogen bonding interactions.     

Both FMNH2 and FADH2 can be utilized by the dibenzothiophene monooxygenase from a desulfurizing bacterium Mycobacterium goodii X7B

To investigate the flavin utilization by dibenzothiophene monooxygenase (DszC), DszC of a desulfurizing bacterium Mycobacterium goodii X7B was purified from the recombinant Escherichia coli. It was shown to be able to utilize either FMNH₂ or FADH₂ when coupled with a flavin reductase that reduces either FMN or FAD. Sequence analysis indicated that DszC was similar to the C₂ component of p-hydroxyphenylacetate hydroxylase from Acinetobacter baumannii, which can use both FADH₂ and FMNH₂ as substrates. Both flavins at high concentrations could inhibit the activity of DszC due to autocatalytic oxidation of reduced flavins. The results suggest that DszC should be reclassified as an FMNH₂ and FADH₂ both-utilizing monooxygenase component and the flavins should be controlled at properly reduced levels to obtain optimal biodesulfurization results.     

Altered mechanism of the alkanesulfonate FMN reductase with the monooxygenase enzyme

The two-component alkanesulfonate monooxygenase system from Escherichia coli is comprised of an FMN reductase (SsuE) and a monooxygenase enzyme (SsuD) that together catalyze the oxidation of alkanesulfonate to the corresponding aldehyde and sulfite products. To determine the effects of protein interactions on catalysis, the steady-state kinetic parameters for SsuE were determined in single-enzyme assays and in the presence of the monooxygenase enzyme and alkanesulfonate substrate. In single-enzyme kinetic assays, SsuE followed an ordered sequential mechanism, with NADPH as the first substrate to bind and NADP+ as the last product to dissociate. However, in the presence of SsuD and octanesulfonate the kinetic mechanism of SsuE is altered to a rapid equilibrium ordered mechanism, and the Km value for FMN is increased 10-fold. These results suggest that both the SsuD enzyme and alkanesulfonate substrate are required to ensure that the FMN reductase reaction proceeds to form the ternary complex with the subsequent generation of reduced flavin transfer.     

Prediction of the enantioselectivity of lipases and esterases by molecular docking method with modified force field parameters

In the work, molecular docking method was applied to extensively predict the enantioselectivity of lipases and esterases. A ligand library consisted of 69 chiral substrates was docked to four lipases and two esterases to set up the prediction model. During the docking process, necessary modification was carried out on van de Waals and hydrogen bond parameters of enzyme/substrate pair so that the ligands were able to adopt productive geometry in the enzymes. The docking results correctly indicated the enantiopreference for 91% (63/69) of docking pairs and the docking energy difference between substrate enantiomers (ΔΔG_docking) was significantly (correlation coefficient = 0.72, P < 0.05) correlated with the activation free energy difference (ΔΔG^≠) that was quantitatively correlated with enantioselectivity of the enzymes. The prediction method was further validated by docking with another 12 enzyme/substrate pairs. Moreover, the prediction error was susceptible to the size of groups bonded to substrate's chiral center and expected ΔΔG^≠ values but was not related to the substrate type and reaction medium. The possible reasons of observed error were discussed. It is demonstrated that the docking method has great application potential in high performance prediction of enzyme enantioselectivity. Biotechnol. Bioeng. 2010. 105: 687–696. © 2009 Wiley Periodicals, Inc.     

The acetylene reduction technique as an assay for nitrogenase activity in the methane oxidizing bacterium Methylococcus capsulatus strain bath

The use of acetylene as a convenient assay substrate for nitrogenase in methane oxidising bacteria is complicated by the observation that it is a potent inhibitor of the methane monooxygenase enzyme in both whole cells and cell-free extracts. If the cells were provided with alternative oxidisable carbon substrates other than methane then nitrogen fixing cells would reduce acetylene to ethylene. Hydrogen gas also served as an oxidisable substrate in the assay. Nitrous oxide, which is reduced by nitrogenase to N₂ and H₂O, was not an inhibitor of methane monooxygenase function and could be used as a convenient assay substrate for nitrogenase. Reduction of both substrates by whole cells showed similar response to oxygen in the assay system and in this respect Methylococcus resembles other free living nitrogen fixing aerobes.     

Computational design of a lipase for catalysis of the Diels-Alder reaction

Combined molecular docking, molecular dynamics (MD) and density functional theory (DFT) studies have been employed to study catalysis of the Diels-Alder reaction by a modified lipase. Six variants of the versatile enzyme Candida Antarctica lipase B (CALB) have been rationally engineered in silico based on the specific characteristics of the pericyclic addition. A kinetic analysis reveals that hydrogen bond stabilization of the transition state and substrate binding are key components of the catalytic process. In the case of substrate binding, which has the greater potential for optimization, both binding strength and positioning of the substrates are important for catalytic efficiency. The binding strength is determined by hydrophobic interactions and can be tuned by careful selection of solvent and substrates. The MD simulations show that substrate positioning is sensitive to cavity shape and size, and can be controlled by a few rational mutations. The well-documented S105A mutation is essential to enable sufficient space in the vicinity of the oxyanion hole. Moreover, bulky residues on the edge of the active site hinders the formation of a sandwich-like nearattack conformer (NAC), and the I189A mutation is needed to obtain enough space above the face of the α,β-double bond on the dienophile. The double mutant S105A/I189A performs quite well for two of three dienophiles. Based on binding constants and NAC energies obtained from MD simulations combined with activation energies from DFT computations, relative catalytic rates (v _cat /v _uncat ) of up to 103 are predicted.     

Analysis of the bacterial luciferase mobile loop by replica-exchange molecular dynamics

Bacterial luciferase contains an extended 29-residue mobile loop. Movements of this loop are governed by binding of either flavin mononucleotide (FMNH₂) or polyvalent anions. To understand this process, loop dynamics were investigated using replica-exchange molecular dynamics that yielded conformational ensembles in either the presence or absence of FMNH₂. The resulting data were analyzed using clustering and network analysis. We observed the closed conformations that are visited only in the simulations with the ligand. Yet the mobile loop is intrinsically flexible, and FMNH₂ binding modifies the relative populations of conformations. This model provides unique information regarding the function of a crystallographically disordered segment of the loop near the binding site. Structures at or near the fringe of this network were compatible with flavin binding or release. Finally, we demonstrate that the crystallographically observed conformation of the mobile loop bound to oxidized flavin was influenced by crystal packing. Thus, our study has revealed what we believe are novel conformations of the mobile loop and additional context for experimentally determined structures.     

Detection of protein-protein interactions in the alkanesulfonate monooxygenase system from Escherichia coli

The two-component alkanesulfonate monooxygenase system utilizes reduced flavin as a substrate to catalyze a unique desulfonation reaction during times of sulfur starvation. The importance of protein-protein interactions in the mechanism of flavin transfer was analyzed in these studies. The results from affinity chromatography and cross-linking experiments support the formation of a stable complex between the flavin mononucleotide (FMN) reductase (SsuE) and monooxygenase (SsuD). Interactions between the two proteins do not lead to overall conformational changes in protein structure, as indicated by the results from circular dichroism spectroscopy in the far-UV region. However, subtle changes in the flavin environment of FMN-bound SsuE that occur in the presence of SsuD were identified by circular dichroism spectroscopy in the visible region. These data are supported by the results from fluorescent spectroscopy experiments, where a dissociation constant of 0.0022 +/- 0.0010 muM was obtained for the binding of SsuE to SsuD. Based on these studies, the stoichiometry for protein-protein interactions is proposed to involve a 1:1 monomeric association of SsuE with SsuD.     

Catalytic importance of the substrate binding order for the FMNH2-dependent alkanesulfonate monooxygenase enzyme

The two-component alkanesulfonate monooxygenase system from Escherichia coli includes an FMN reductase (SsuE) and an FMNH2-dependent alkanesulfonate monooxygenase (SsuD) involved in the acquisition of sulfur from alkanesulfonates during sulfur starvation. The SsuD enzyme directly catalyzes the oxidation of alkanesulfonate to aldehyde and sulfite in the presence of O2 and FMNH2. The goal of these studies was to investigate the kinetic mechanism of SsuD through rapid reaction kinetics and substrate binding studies. The SsuD enzyme shows a clear preference for FMNH2 (Kd, 0.32 +/- 0.15 microM) compared to FMN (Kd, 10.2 +/- 0.4 microM) with a 1:1 binding stoichiometry for each form of the flavin. The kinetic trace of premixed SsuD and FMNH2 mixed with oxygenated buffer was best fit to a double exponential with no observed formation of the C4a-(hydro)peroxyflavin. However, when FMNH2 was mixed with SsuD and oxygenated buffer an initial fast phase (kobs, 12.9 s-1) was observed, suggesting that the mixing order is critical for the accumulation of the C4a-(hydro)peroxyflavin. Results from fluorimetric titrations with octanesulfonate imply that reduced flavin must bind first to promote octanesulfonate binding. When octanesulfonate was included in the kinetic studies the C4a-(hydro)peroxyflavin was observed at 370 nm when FMNH2 was not premixed with SsuD, which correlated with an increase in octanal product. There was a clear hyperbolic dependence on octanesulfonate binding, indicating that octanesulfonate binds in rapid equilibrium, and further results indicated there was a second isomerization step following binding. These results suggest that an ordered substrate binding mechanism is important in the desulfonation reaction by SsuD with reduced flavin binding first followed by either O2 or octanesulfonate.     

Deletional studies to investigate the functional role of a dynamic loop region of alkanesulfonate monooxygenase

Several bacterial organisms rely on the two-component alkanesulfonate monooxygenase system for the acquisition of organosulfonate compounds when inorganic sulfur is limiting in the environment. This system is comprised of an FMN reductase (SsuE) that supplies reduced flavin to the alkanesulfonate monooxygenase (SsuD). Desulfonation of alkanesulfonates by SsuD is catalyzed through the activation of dioxygen by reduced flavin. The three-dimensional structure of SsuD exists as a TIM-barrel fold with several discrete insertion regions. An extensive insertion region near the putative active site was disordered in the SsuD structure, suggesting the importance of protein dynamics in the desulfonation mechanism. Three variants containing a partial deletion of the loop region were constructed to evaluate the functional properties of this region. There were no overall gross changes in secondary structure for the three SsuD deletion variants compared to wild-type SsuD, but each variant was found to be catalytically inactive. The deletion variants were unable to undergo the conformational changes necessary for catalysis even though they were able to bind reduced flavin. Rapid kinetic analyses monitoring the reductive and oxidative half-reactions indicated that the SsuD deletion variants failed to protect reduced flavin from unproductive oxidation. These studies define the importance of dynamic loop region for protection and stabilization of reduced flavin and reaction intermediates.     

Docking of 4‐oxalocrotonate tautomerase substrates: Implications for the catalytic mechanism

The enzyme 4‐oxalocrotonate tautomerase catalyzes the ketonization of dienols, which after further processing become intermediates in the Krebs cycle. The enzyme uses a general acid–base mechanism for proton transfer: the amino‐terminal proline has been shown to function as the catalytic base and Arg39 has been implicated as the catalytic acid. We report the results of molecular docking simulations of 4‐oxalocrotonate tautomerase with two substrates, 2‐hydroxymuconate and 5‐carboxymethyl‐2‐hydroxymuconate. pK_a calculations are also performed for the free enzyme. The predicted binding mode of 2‐hydroxymuconate is in agreement with experimental data. A model for the binding mode of 5‐carboxymethyl‐2‐hydroxymuconate is proposed which explains the lower catalytic efficiency of the enzyme toward this substrate. The pK_a predictions and docking simulations support residue Arg39 as the general acid for the enzyme catalysis. © 1999 John Wiley & Sons, Inc. Biopoly 50: 319–328, 1999     

Exploration of human serum albumin binding sites by docking and molecular dynamics flexible ligand–protein interactions

Five‐nanosecond molecular dynamics (MD) simulations were performed on human serum albumin (HSA) to study the conformational features of its primary ligand binding sites (I and II). Additionally, 11 HSA snapshots were extracted every 0.5 ns to explore the binding affinity (K_d) of 94 known HSA binding drugs using a blind docking procedure. MD simulations indicate that there is considerable flexibility for the protein, including the known sites I and II. Movements at HSA sites I and II were evidenced by structural analyses and docking simulations. The latter enabled the study and analysis of the HSA–ligand interactions of warfarin and ketoprofen (ligands binding to sites I and II, respectively) in greater detail. Our results indicate that the free energy values by docking (K_d observed) depend upon the conformations of both HSA and the ligand. The 94 HSA–ligand binding K_d values, obtained by the docking procedure, were subjected to a quantitative structure‐activity relationship (QSAR) study by multiple regression analysis. The best correlation between the observed and QSAR theoretical (K_d predicted) data was displayed at 2.5 ns. This study provides evidence that HSA binding sites I and II interact specifically with a variety of compounds through conformational adjustments of the protein structure in conjunction with ligand conformational adaptation to these sites. These results serve to explain the high ligand‐promiscuity of HSA. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 161–170, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Regiospecificity of Human Cytochrome P450 1A1-Mediated Oxidations: The Role of Steric Effects

Abstract

Cytochrome P450 1A1 oxidizes a diverse range of substrates, including the procarcinogenic xenobiotic benzo[a]pyrene (B[a]P) and endogenous fatty acid precursors of prostaglandins, such as arachidonic acid (AA) and eicosapentaenoic acid (EA). We have investigated the extent to which enzyme-substrate interactions govern regio- and stereoselectivity of oxidation of these compounds by using docking and molecular dynamics (MD) simulations to examine the likelihood of substrate oxidation at various sites. Due to structural differences between the substrates analyzed, B[a]P and its diols (planar, rigid), and the fatty acids AA and EA (long, flexible), different docking strategies were required. B[a]P, B[a]P-7,8-diols, (+) 7S,85- and (-) 7R,8R-diols, were docked into the active site of a homology model of P450 1A1 using an automated routine. Affinity (Accelrys, San Diego, CA). AA and EA, on the other hand, required a series of restrained MD simulations to obtain a variety of productive binding modes. All complexes were evaluated by MD-based in silico site scoring to predict product profiles based on certain geometric criteria, such as angle and distance of a given substrate atom from the ferryl oxygen. For all substrates studied, the in vitro profiles were generally reflected by the in silico scores, which suggests that steric factors play a key role in determining regiospecificity in P450 1A1-mediated oxidations. We have also shown that molecular dynamics simulations may be very useful in determination of product profiles for structurally diverse substrates of P450 enzymes.     

Functional role of a conserved arginine residue located on a mobile loop of alkanesulfonate monooxygenase

The structure of the flavin-dependent alkanesulfonate monooxygenase (SsuD) exists as a TIM-barrel structure with an insertion region located over the active site that contains a conserved arginine (Arg297) residue present in all SsuD homologues. Substitution of Arg297 with alanine (R297A SsuD) or lysine (R297K SsuD) was performed to determine the functional role of this conserved residue in SsuD catalysis. While the more conservative R297K SsuD possessed a lower k(cat)/K(m) value (0.04 ± 0.01 μM(-1) min(-1)) relative to wild-type (1.17 ± 0.22 μM(-1) min(-1)), there was no activity observed with the R297A SsuD variant. Each of the arginine variants had similar K(d) values for flavin binding as wild-type SsuD (0.32 ± 0.15 μM), but there was no measurable binding of octanesulfonate. The low levels of activity for the R297A and R297K SsuD variants correlated with the absence of any detectable C4a-(peroxy)flavin formation in stopped-flow kinetic studies. Single-turnover experiments were performed in the presence of SsuE to evaluate both the reductive and oxidative half-reaction. With wild-type SsuD a lag phase is observed following the reductive half-reaction by SsuE that represents flavin transfer or conformational changes associated with the binding of substrates. Evaluation of the Arg297 SsuD variants in the presence of SsuE showed no lag phase following reduction by SsuE, and the flavin was oxidized immediately following the reductive half-reaction. These results corresponded with a lack of detectable changes in the proteolytic susceptibility of R297A and R297K SsuD in the presence of reduced flavin and/or octanesulfonate, signifying the absence of a conformational change in these variants with the substitution of Arg297.     

Relationship between mutation of serine residue at 315th position in M. tuberculosis catalase-peroxidase enzyme and Isoniazid susceptibility: An in silico analysis

Remarkable advances have been made in the drug therapy of tuberculosis. However much remains to be learned about the molecular and structural basis of drug resistance in Mycobacterium tuberculosis. It is known that, activation of Isoniazid (INH) is mediated by Mycobacterium tuberculosis catalase-peroxidase (MtBKatG) and mutation at position 315 (serine to threonine) leads to resistance. We have conducted studies on the drug resistance through docking and binding analysis supported by time-scale (∼1000 ps) and unrestrained all-atom molecular dynamics simulations of wild and mutant MtBKatG. The study showed conformational changes of binding residues. Mutant (S315T) showed high docking score and INH binding affinity as compared to wild enzyme. In molecular dynamics simulation, mutant enzyme exhibited less structure fluctuation at INH binding residues and more degree of fluctuation at C-terminal domain compared to wild enzyme. Our computational studies and data endorse that MtBKatG mutation (S315T) decrease the flexibility of binding residues and made them rigid by altering the conformational changes, in turn it hampers the INH activity. We ascertain from this work that, this study on structural mechanism of resistance development in Mycobacterium tuberculosis would lead to new therapeutics based on the result obtained in this study.