Glutarate Hydroxylation by the Carbon Starvation-Induced Protein D: A Computational Study into the Stereo- and Regioselectivities of the Reaction

Adaptive laboratory evolution and transcriptomics-guided engineering of Escherichia coli for increased isobutanol tolerance

As a renewable energy from biomass, isobutanol is considered as a promising alternative to fossil fuels. To biotechnologically produce isobutanol, strain development using industrial microbial hosts, such as Escherichia coli, has been conducted by introducing a heterologous isobutanol synthetic pathway. However, the toxicity of produced isobutanol inhibits cell growth, thereby restricting improvements in isobutanol titer, yield, and productivity. Therefore, the development of robust microbial strains tolerant to isobutanol is required. In this study, isobutanol-tolerant mutants were isolated from two E. coli parental strains, E. coli BL21(DE3) and MG1655(DE3), through adaptive laboratory evolution (ALE) under high isobutanol concentrations. Subsequently, 16 putative genes responsible for isobutanol tolerance were identified by transcriptomic analysis. When overexpressed in E. coli, four genes (fadB, dppC, acs, and csiD) conferred isobutanol tolerance. A fermentation study with a reverse engineered isobutanol-producing E. coli JK209 strain showed that fadB or dppC overexpression improved isobutanol titers by 1.5 times, compared to the control strain. Through coupling adaptive evolution with transcriptomic analysis, new genetic targets utilizable were identified as the basis for the development of an isobutanol-tolerant strain. Thus, these new findings will be helpful not only for a fundamental understanding of microbial isobutanol tolerance but also for facilitating industrially feasible isobutanol production.

The Mechanism of Biochemical NO-Sensing: Insights from Computational Chemistry

The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO-receptor in humans is soluble guanylyl cyclase (sGC) - a histidine-ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO-signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high-level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H-NOX protein from Shewanella oneidensis (So H-NOX) homologous to the NO-binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six-coordinate intermediate heme-NO complex, followed by relatively facile His decoordination yielding a five-coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.

What Drives Radical Halogenation versus Hydroxylation in Mononuclear Nonheme Iron Complexes? A Combined Experimental and Computational Study

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C–H bond of a substrate to convert it into C–X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe^III(OH)(X) intermediate in their catalytic cycle. The analogous Fe^III(OH) intermediate in nonheme iron hydroxylases transfers OH^• to give alcohol product, whereas the halogenases transfer X^• to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe^III(BNPA^Ph2O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R^•) is described. It is found that X^• transfer occurs with a secondary carbon radical, as opposed to OH^• transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C–X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

Structural Insight into the Catalytic Mechanism of the Endoperoxide Synthase FtmOx1

The 2-oxoglutarate (2OG)-dependent non-heme enzyme FtmOx1 catalyzes the endoperoxide biosynthesis of verruculogen. Although several mechanistic studies have been carried out, the catalytic mechanism of FtmOx1 is not well determined owing to the lack of a reliable complex structure of FtmOx1 with fumitremorgin B. Herein we provide the X-ray crystal structure of the ternary complex FtmOx1⋅2OG⋅fumitremorgin B at a resolution of 1.22 Å. Our structures show that the binding of fumitremorgin B induces significant compression of the active pocket and that Y68 is in close proximity to C26 of the substrate. Further MD simulation and QM/MM calculations support a CarC-like mechanism, in which Y68 acts as the H atom donor for quenching the C26-centered substrate radical. Our results are consistent with all available experimental data and highlight the importance of accurate complex structures in the mechanistic study of enzymatic catalysis.

Deciphering the Reaction Pathway of Mononuclear Iron Enzyme-Catalyzed N≡C Triple Bond Formation in Isocyanide Lipopeptide and Polyketide Biosynthesis

Despite the diversity of reactions catalyzed by 2-oxoglutarate-dependent nonheme iron (Fe/2OG) enzymes identified in recent years, only a limited number of these enzymes have been investigated in mechanistic detail. In particular, several Fe/2OG-dependent enzymes capable of catalyzing isocyanide formation have been reported. While the glycine moiety has been identified as a biosynthon for the isocyanide group, how the actual conversion is effected remains obscure. To elucidate the catalytic mechanism, we characterized two previously unidentified (AecA and AmcA) along with two known (ScoE and SfaA) Fe/2OG-dependent enzymes that catalyze N≡C triple bond installation using synthesized substrate analogues and potential intermediates. Our results indicate that isocyanide formation likely entails a two-step sequence involving an imine intermediate that undergoes decarboxylation-assisted desaturation to yield the isocyanide product. Results obtained from the in vitro experiments are further supported by mutagenesis, the product-bound enzyme structure, and in silico analysis.

Electrostatic Perturbations in the Substrate-Binding Pocket of Taurine/α-Ketoglutarate Dioxygenase Determine its Selectivity

Taurine/α-ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C1 -position. Recent computational studies showed that in the gas-phase the C2 -H bond of taurine is substantially weaker than the C1 -H bond, yet no evidence exists of 2-hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second-coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C1 -H and C2 -H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro-R C1 -H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

Local Charge Distributions, Electric Dipole Moments, and Local Electric Fields Influence Reactivity Patterns and Guide Regioselectivities in α-Ketoglutarate-Dependent Non-heme Iron Dioxygenases

Non-heme iron dioxygenases catalyze vital processes for human health related to the biosynthesis of essential products and the biodegradation of toxic metabolites. Often the natural product biosyntheses by these non-heme iron dioxygenases is highly regio- and chemoselective, which are commonly assigned to tight substrate-binding and positioning. However, recent high-level computational modeling has shown that substrate-binding and positioning is only part of the story and long-range electrostatic interactions can play a major additional role.In this Account, we review and summarize computational viewpoints on the high regio- and chemoselectivity of α-ketoglutarate-dependent non-heme iron dioxygenases and how external perturbations affect the catalysis. In particular, studies from our groups have shown that often a regioselectivity in enzymes can be accomplished by stabilization of the rate-determining transition state for the reaction through external charges, electric dipole moments, or local electric field effects. Furthermore, bond dissociation energies in molecules are shown to be influenced by an electric field effect, and through targeting a specific bond in an electric field, this can lead to an unusually specific reaction. For instance, in the carbon-induced starvation protein, we studied two substrate-bound conformations and showed that regardless of what C-H bond of the substrate is closest to the iron(IV)-oxo oxidant, the lowest hydrogen atom abstraction barrier is always for the pro-S C²-H abstraction due to an induced dipole moment of the protein that weakens this bond. In another example of the hygromycin biosynthesis enzyme, an oxidative ring-closure reaction in the substrate forms an ortho-δ-ester ring. Calculations on this enzyme show that the selectivity is guided by a protonated lysine residue in the active site that, through its positive charge, triggers a low energy hydrogen atom abstraction barrier. A final set of examples in this Account discuss the viomycin biosynthesis enzyme and the 2-(trimethylammonio)ethylphosphonate dioxygenase (TmpA) enzyme. Both of these enzymes are shown to possess a significant local dipole moment and local electric field effect due to charged residues surrounding the substrate and oxidant binding pockets. The protein dipole moment and local electric field strength changes the C-H bond strengths of the substrate as compared to the gas-phase triggers the regioselectivity of substrate activation. In particular, we show that in the gas phase and in a protein environment C-H bond strengths are different due to local electric dipole moments and electric field strengths. These examples show that enzymes have an intricately designed structure that enables a chemical reaction under ambient conditions through the positioning of positively and negatively charged residues that influence and enhance reaction mechanisms. These computational insights create huge possibilities in bioengineering to apply local electric field and dipole moments in proteins to achieve an unusual selectivity and specificity and trigger a fit-for-purpose biocatalyst for unique biotransformations.

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.