Transition States in a Protein Environment − ONIOM QM:MM Modeling of Isopenicillin N Synthesis

Structural and mechanistic insights into Quinolone Synthase to address its functional promiscuity

Quinolone synthase from Aegle marmelos (AmQNS) is a type III polyketide synthase that yields therapeutically effective quinolone and acridone compounds. Addressing the structural and molecular underpinnings of AmQNS and its substrate interaction in terms of its high selectivity and specificity can aid in the development of numerous novel compounds. This paper presents a high-resolution AmQNS crystal structure and explains its mechanistic role in synthetic selectivity. Additionally, we provide a model framework to comprehend structural constraints on ketide insertion and postulate that AmQNS's steric and electrostatic selectivity plays a role in its ability to bind to various core substrates, resulting in its synthetic diversity. AmQNS prefers quinolone synthesis and can accommodate large substrates because of its wide active site entrance. However, our research suggests that acridone is exclusively synthesized in the presence of high malonyl-CoA concentrations. Potential implications of functionally relevant residue mutations were also investigated, which will assist in harnessing the benefits of mutations for targeted polyketide production. The pharmaceutical industry stands to gain from these findings as they expand the pool of potential drug candidates, and these methodologies can also be applied to additional promising enzymes.

Synthesis of efficient bi-anchoring bifuran/biphenyl derivatives for dye-sensitized solar cell applications†

The synthesis, description, and demonstration of dye-sensitive solar cell sensitizers containing bifuran/biphenyl derivatives with cyanoacetic acid, barbiturate, thiobarbituric acid, and 4-carboxylcyanoacetamides have been reported. A photovoltaic performance measurement was conducted using the Ru(ii) dye N3 as a reference to examine the effects of different electron acceptor units and replacement of the π-spacer bifuran by biphenyl units on the photophysical, electrochemical, and photovoltaic properties of eight new distinct organic dyes HB-1–8. The new organic dyes HB-1–8 were prepared and compared with the N3 metal dye. Density functional theory calculations were carried out to explore the ground state geometrical structures and electronic structures of the eight dyes. Under standard global AM 1.5 solar condition, the solar cells based on HB-1–8 show the overall power conversion efficiencies of 2.93–5.51%. The presented research shows that the organic dye photovoltaic performances can vary greatly depending on the type of electron donor and acceptor used. Dye HB-3 exhibited the highest efficiency among the eight investigated dyes, reaching 5.51% with a V_OC value higher than N3.

We have synthesized new metal-free dyes (HB-1–8) containing bifuran/biphenyl derivatives for dye-sensitized solar cells. Dye HB-3 exhibited the highest efficiency among the eight investigated dyes, reaching 5.51% with a V_OC value higher than N3.

Properties of the Reactants and Their Interactions within and with the Enzyme Binding Cavity Determine Reaction Selectivities. The Case of Fe(II)/2-Oxoglutarate Dependent Enzymes

Fe(II)/2-oxoglutarate dependent dioxygenases (ODDs) share a double stranded beta helix (DSBH) fold and utilise a common reactive intermediate, ferryl species, to catalyse oxidative transformations of substrates. Despite the structural similarities, ODDs accept a variety of substrates and facilitate a wide range of reactions, that is hydroxylations, desaturations, (oxa)cyclisations and ring rearrangements. In this review we present and discuss the factors contributing to the observed (regio)selectivities of ODDs. They span from inherent properties of the reactants, that is, substrate molecule and iron cofactor, to the interactions between the substrate and the enzyme's binding cavity; the latter can counterbalance the effect of the former. Based on results of both experimental and computational studies dedicated to ODDs, we also line out the properties of the reactants which promote reaction outcomes other than the "default" hydroxylation. It turns out that the reaction selectivity depends on a delicate balance of interactions between the components of the investigated system.

X-ray free-electron laser studies reveal correlated motion during isopenicillin N synthase catalysis

Isopenicillin N synthase (IPNS) catalyzes the unique reaction of l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) with dioxygen giving isopenicillin N (IPN), the precursor of all natural penicillins and cephalosporins. X-ray free-electron laser studies including time-resolved crystallography and emission spectroscopy reveal how reaction of IPNS:Fe(II):ACV with dioxygen to yield an Fe(III) superoxide causes differences in active site volume and unexpected conformational changes that propagate to structurally remote regions. Combined with solution studies, the results reveal the importance of protein dynamics in regulating intermediate conformations during conversion of ACV to IPN. The results have implications for catalysis by multiple IPNS-related oxygenases, including those involved in the human hypoxic response, and highlight the power of serial femtosecond crystallography to provide insight into long-range enzyme dynamics during reactions presently impossible for nonprotein catalysts.

Functionalized Au15 nanoclusters as luminescent probes for protein carbonylation detection

Atomically precise, ligand-protected gold nanoclusters (AuNCs) attract considerable attention as contrast agents in the biosensing field. However, the control of their optical properties and functionalization of surface ligands remain challenging. Here we report a strategy to tailor AuNCs for the precise detection of protein carbonylation-a causal biomarker of ageing. We produce Au₁₅SG₁₃ (SG for glutathione) with atomic precision and functionalize it with a thiolated aminooxy moiety to impart protein carbonyl-binding properties. Mass spectrometry and molecular modelling reveal the key structural features of Au₁₅SG₁₂-Aminooxy and its reactivity towards carbonyls. Finally, we demonstrate that Au₁₅SG₁₂-Aminooxy detects protein carbonylation in gel-based 1D electrophoresis by one- and two-photon excited fluorescence. Importantly, to our knowledge, this is the first application of an AuNC that detects a post-translational modification as a nonlinear optical probe. The significance of post-translational modifications in life sciences may open avenues for the use of Au₁₅SG₁₃ and other nanoclusters as contrast agents with tailored surface functionalization and optical properties.

Isopenicillin N Synthase: Crystallographic Studies

Isopenicillin N synthase (IPNS) is a non-heme iron oxidase (NHIO) that catalyses the cyclisation of tripeptide δ-(l-α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) to bicyclic isopenicillin N (IPN). Over the last 25 years, crystallography has shed considerable light on the mechanism of IPNS catalysis. The first crystal structure, for apo-IPNS with Mn bound in place of Fe at the active site, reported in 1995, was also the first structure for a member of the wider NHIO family. This was followed by the anaerobic enzyme-substrate complex IPNS-Fe-ACV (1997), this complex plus nitric oxide as a surrogate for co-substrate dioxygen (1997), and an enzyme product complex (1999). Since then, crystallography has been used to probe many aspects of the IPNS reaction mechanism, by crystallising the protein with a diversity of substrate analogues and triggering the oxidative reaction by using elevated oxygen pressures to force the gaseous co-substrate throughout protein crystals and maximise synchronicity of turnover in crystallo. In this way, X-ray structures have been elucidated for a range of complexes closely related to and/or directly derived from key intermediates in the catalytic cycle, thereby answering numerous mechanistic questions that had arisen from solution-phase experiments, and posing many new ones. The results of these crystallographic studies have, in turn, informed computational experiments that have brought further insight. These combined crystallographic and computational investigations augment and extend the results of earlier spectroscopic analyses and solution phase studies of IPNS turnover, to enrich our understanding of this important protein and the wider NHIO enzyme family.

Temperature-Dependent Reactivity of a Non-heme FeIII(OH)(SR) Complex: Relevance to Isopenicillin N Synthase

Non-heme iron complexes with cis-Fe^III(OH)(SAr/OAr) coordination were isolated and examined for their reactivity with a tertiary carbon radical. The sulfur-ligated complex shows a temperature dependence on ^•OH versus ArS^• transfer, whereas the oxygen-ligated complex does not. These results provide the first working model for C-S bond formation in isopenicillin N synthase and indicate that kinetic control may be a key factor in the selectivity of non-heme iron "rebound" processes.

Reveal the interaction mechanism of five old drugs targeting VEGFR2 through computational simulations

VEGFR2, vascular endothelial growth factor receptor 2, plays an important role in anti-angiogenesis and is an effective target for inhibiting tumor cell proliferation and metastasis. Many small molecule inhibitors have so far exhibited fine therapeutic effects but do not rule out some adverse reactions. From the perspective of the new use of old drugs, we use a combination of two different docking methods, molecular dynamics simulations and quantum-chemical calculations to acquire potential anti-angiogenesis inhibitors from the library of FDA-approved drugs. We attain five FDA-approved old drugs from Drugbank as potential inhibitors against VEGFR2. Therein, the anti-tumor effects of three compounds, including vilazodone (psychiatric drug), pranlukast and zafirlukast (asthma drugs), have been reported by previous experiments but no anti-tumor data is available for the other two compounds, including antrafenine (analgesic and anti-inflammatory drug) and iloperidone (psychiatric drug). These five compounds exhibit more stable interaction than sorafenib as a market-oriented drug targeting VEGFR2. In parallel, there is a most stable interaction for zafirlukast while a weakest interaction for iloperidone with VEGFR2. We show that these five compounds bind with the hydrophobic cavity of VEGFR2, then forming hydrogen bond interactions with three key residues, Glu-885, Cys-919 and Asp-1046. Lys-868 and Phe-1047 play an important role in stabilizing the interaction conformation. The binding poses of pranlukast and vilazodone are similar to that of sorafenib, whereas antrafenine and zafirlukast act differently from sorafenib, focusing on the direction difference of the respective ring structure. This work may help to develop new and effective anti-angiogenic inhibitors.

The Driving Force for the Acylation of β-Lactam Antibiotics by L,D-Transpeptidase 2: Quantum Mechanics/Molecular Mechanics (QM/MM) Study

β-lactam antibiotics, which are used to treat infectious diseases, are currently the most widely used class of antibiotics. This study focused on the chemical reactivity of five- and six-membered ring systems attached to the β-lactam ring. The ring strain energy (RSE), force constant (FC) of amide (C-N), acylation transition states and second-order perturbation stabilization energies of 13 basic structural units of β-lactam derivatives were computed using the M06-2X and G3/B3LYP multistep method. In the ring strain calculations, an isodesmic reaction scheme was used to obtain the total energies. RSE is relatively greater in the five-(1a-2c) compared to the six-membered ring systems except for 4b, which gives a RSE that is comparable to five-membered ring lactams. These variations were also observed in the calculated inter-atomic amide bond distances (C-N), which is why the six-membered ring lactams C-N bond are more rigid than those with five-membered ring lactams. The calculated ΔG^# values from the acylation reaction of the lactams (involving the S-H group of the cysteine active residue from L,D transpeptidase 2) revealed a faster rate of C-N cleavage in the five-membered ring lactams especially in the 1-2 derivatives (17.58 kcal mol^-1 ). This observation is also reflected in the calculated amide bond force constant (1.26 mDyn/A) indicating a weaker bond strength, suggesting that electronic factors (electron delocalization) play more of a role on reactivity of the β-lactam ring, than ring strain.

Structural basis for the recognition of sulfur in phosphorothioated DNA

There have been very few reports on protein domains that specifically recognize sulfur. Here we present the crystal structure of the sulfur-binding domain (SBD) from the DNA phosphorothioation (PT)-dependent restriction endonuclease ScoMcrA. SBD contains a hydrophobic surface cavity that is formed by the aromatic ring of Y164, the pyrolidine ring of P165, and the non-polar side chains of four other residues that serve as lid, base, and wall of the cavity. The SBD and PT-DNA undergo conformational changes upon binding. The S¹⁸⁷RGRR¹⁹¹ loop inserts into the DNA major groove to make contacts with the bases of the G_PSGCC core sequence. Mutating key residues of SBD impairs PT-DNA association. More than 1000 sequenced microbial species from fourteen phyla contain SBD homologs. We show that three of these homologs bind PT-DNA in vitro and restrict PT-DNA gene transfer in vivo. These results show that SBD-like PT-DNA readers exist widely in prokaryotes.

Spectroscopic and Electronic Structure Study of ETHE1: Elucidating the Factors Influencing Sulfur Oxidation and Oxygenation in Mononuclear Nonheme Iron Enzymes

ETHE1 is a member of a growing subclass of nonheme Fe enzymes that catalyzes transformations of sulfur-containing substrates without a cofactor. ETHE1 dioxygenates glutathione persulfide (GSSH) to glutathione (GSH) and sulfite in a reaction which is similar to that of cysteine dioxygenase (CDO), but with monodentate (vs bidentate) substrate coordination and a 2-His/1-Asp (vs 3-His) ligand set. In this study, we demonstrate that GSS^- binds directly to the iron active site, causing coordination unsaturation to prime the site for O₂ activation. Nitrosyl complexes without and with GSSH were generated and spectroscopically characterized as unreactive analogues for the invoked ferric superoxide intermediate. New spectral features from persulfide binding to the Fe^III include the appearance of a low-energy Fe^III ligand field transition, an energy shift of a NO^- to Fe^III CT transition, and two new GSS^- to Fe^III CT transitions. Time-dependent density functional theory calculations were used to simulate the experimental spectra to determine the persulfide orientation. Correlation of these spectral features with those of monodentate cysteine binding in isopenicillin N synthase (IPNS) shows that the persulfide is a poorer donor but still results in an equivalent frontier molecular orbital for reactivity. The ETHE1 persulfide dioxygenation reaction coordinate was calculated, and while the initial steps are similar to the reaction coordinate of CDO, an additional hydrolysis step is required in ETHE1 to break the S-S bond. Unlike ETHE1 and CDO, which both oxygenate sulfur, IPNS oxidizes sulfur through an initial H atom abstraction. Thus, factors that determine oxygenase vs oxidase reactivity were evaluated. In general, sulfur oxygenation is thermodynamically favored and has a lower barrier for reactivity. However, in IPNS, second-sphere residues in the active site pocket constrain the substrate, raising the barrier for sulfur oxygenation relative to oxidation via H atom abstraction.

The catalytic role of water in the binding site of l,d-transpeptidase 2 within acylation mechanism: A QM/MM (ONIOM) modelling

Mycobacterium tuberculosis is the causative agent of Tuberculosis. Formation of 3 → 3 crosslinks in the peptidoglycan layer of M. tuberculosis is catalyzed by l,d-transpeptidases. These enzymes can confer resistance against classical β-lactams that inhibit enzymes that generate 4 → 3 peptidoglycan crosslinks. The focus of this study is to investigate the catalytic role of water molecules in the acylation mechanism of the β-lactam ring within two models; 4- and 6-membered ring systems using two-layered our Own N-layer integrated Molecular Mechanics ONIOM (B3LYP/6-311++G(2d,2p): AMBER) model. The obtained thermochemical parameters revealed that the 6-membered ring model best describes the inhibition mechanism of acylation which indicates the role of water in the preference of 6-membered ring reaction pathway. This finding is in accordance with experimental data for the rate-limiting step of cysteine protease with the same class of inhibitor and binding affinity for both inhibitors. As expected, the ΔG^# results also reveal that the 6-membered ring reaction pathway is the most favourable. The electrostatic potential (ESP) and the natural bond orbital analysis (NBO) showed stronger interactions in 6-membered ring transition state (TS-6) mechanism involving water in the active site of the enzyme. This study could be helpful in the development of novel antibiotics against l,d-transpeptidase.

Roles of 2-oxoglutarate oxygenases and isopenicillin N synthase in β-lactam biosynthesis

Covering: up to 2017 2-Oxoglutarate (2OG) dependent oxygenases and the homologous oxidase isopenicillin N synthase (IPNS) play crucial roles in the biosynthesis of β-lactam ring containing natural products. IPNS catalyses formation of the bicyclic penicillin nucleus from a tripeptide. 2OG oxygenases catalyse reactions that diversify the chemistry of β-lactams formed by both IPNS and non-oxidative enzymes. Reactions catalysed by the 2OG oxygenases of β-lactam biosynthesis not only involve their typical hydroxylation reactions, but also desaturation, epimerisation, rearrangement, and ring-forming reactions. Some of the enzymes involved in β-lactam biosynthesis exhibit remarkable substrate and product selectivities. We review the roles of 2OG oxygenases and IPNS in β-lactam biosynthesis, highlighting opportunities for application of knowledge of their roles, structures, and mechanisms.

Parametrization of Combined Quantum Mechanical and Molecular Mechanical Methods: Bond-Tuned Link Atoms

Combined quantum mechanical and molecular mechanical (QM/MM) methods are the most powerful available methods for high-level treatments of subsystems of very large systems. The treatment of the QM−MM boundary strongly affects the accuracy of QM/MM calculations. For QM/MM calculations having covalent bonds cut by the QM−MM boundary, it has been proposed previously to use a scheme with system-specific tuned fluorine link atoms. Here, we propose a broadly parametrized scheme where the parameters of the tuned F link atoms depend only on the type of bond being cut. In the proposed new scheme, the F link atom is tuned for systems with a certain type of cut bond at the QM−MM boundary instead of for a specific target system, and the resulting link atoms are call bond-tuned link atoms. In principle, the bond-tuned link atoms can be as convenient as the popular H link atoms, and they are especially well adapted for high-throughput and accurate QM/MM calculations. Here, we present the parameters for several kinds of cut bonds along with a set of validation calculations that confirm that the proposed bond-tuned link-atom scheme can be as accurate as the system-specific tuned F link-atom scheme.

Unique chemistry of non-heme iron enzymes in fungal biosynthetic pathways

Reactions by non-heme iron enzymes in structurally intriguing fungal natural products pathways are summarized and discussed.

Terminally Truncated Isopenicillin N Synthase Generates a Dithioester Product: Evidence for a Thioaldehyde Intermediate during Catalysis and a New Mode of Reaction for Non-Heme Iron Oxidases

Isopenicillin N synthase (IPNS) catalyses the four-electron oxidation of a tripeptide, l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV), to give isopenicillin N (IPN), the first-formed β-lactam in penicillin and cephalosporin biosynthesis. IPNS catalysis is dependent upon an iron(II) cofactor and oxygen as a co-substrate. In the absence of substrate, the carbonyl oxygen of the side-chain amide of the penultimate residue, Gln330, co-ordinates to the active-site metal iron. Substrate binding ablates the interaction between Gln330 and the metal, triggering rearrangement of seven C-terminal residues, which move to take up a conformation that extends the final α-helix and encloses ACV in the active site. Mutagenesis studies are reported, which probe the role of the C-terminal and other aspects of the substrate binding pocket in IPNS. The hydrophobic nature of amino acid side-chains around the ACV binding pocket is important in catalysis. Deletion of seven C-terminal residues exposes the active site and leads to formation of a new type of thiol oxidation product. The isolated product is shown by LC-MS and NMR analyses to be the ene-thiol tautomer of a dithioester, made up from two molecules of ACV linked between the thiol sulfur of one tripeptide and the oxidised cysteinyl β-carbon of the other. A mechanism for its formation is proposed, supported by an X-ray crystal structure, which shows the substrate ACV bound at the active site, its cysteinyl β-carbon exposed to attack by a second molecule of substrate, adjacent. Formation of this product constitutes a new mode of reaction for IPNS and non-heme iron oxidases in general.

Multistructural microiteration technique for geometry optimization and reaction path calculation in large systems

We propose a multistructural microiteration (MSM) method for geometry optimization and reaction path calculation in large systems. MSM is a simple extension of the geometrical microiteration technique. In conventional microiteration, the structure of the non-reaction-center (surrounding) part is optimized by fixing atoms in the reaction-center part before displacements of the reaction-center atoms. In this method, the surrounding part is described as the weighted sum of multiple surrounding structures that are independently optimized. Then, geometric displacements of the reaction-center atoms are performed in the mean field generated by the weighted sum of the surrounding parts. MSM was combined with the QM/MM-ONIOM method and applied to chemical reactions in aqueous solution or enzyme. In all three cases, MSM gave lower reaction energy profiles than the QM/MM-ONIOM-microiteration method over the entire reaction paths with comparable computational costs. © 2017 Wiley Periodicals, Inc.

Dioxygen activation by nonheme iron enzymes with the 2-His-1-carboxylate facial triad that generate high-valent oxoiron oxidants

The 2-His-1-carboxylate facial triad is a widely used scaffold to bind the iron center in mononuclear nonheme iron enzymes for activating dioxygen in a variety of oxidative transformations of metabolic significance. Since the 1990s, over a hundred different iron enzymes have been identified to use this platform. This structural motif consists of two histidines and the side chain carboxylate of an aspartate or a glutamate arranged in a facial array that binds iron(II) at the active site. This triad occupies one face of an iron-centered octahedron and makes the opposite face available for the coordination of O₂ and, in many cases, substrate, allowing the tailoring of the iron-dioxygen chemistry to carry out a plethora of diverse reactions. Activated dioxygen-derived species involved in the enzyme mechanisms include iron(III)-superoxo, iron(III)-peroxo, and high-valent iron(IV)-oxo intermediates. In this article, we highlight the major crystallographic, spectroscopic, and mechanistic advances of the past 20 years that have significantly enhanced our understanding of the mechanisms of O₂ activation and the key roles played by iron-based oxidants.

Beyond ferryl-mediated hydroxylation: 40 years of the rebound mechanism and C-H activation

Since our initial report in 1976, the oxygen rebound mechanism has become the consensus mechanistic feature for an expanding variety of enzymatic C-H functionalization reactions and small molecule biomimetic catalysts. For both the biotransformations and models, an initial hydrogen atom abstraction from the substrate (R-H) by high-valent iron-oxo species (Feⁿ=O) generates a substrate radical and a reduced iron hydroxide, [Fe^n-1-OH ·R]. This caged radical pair then evolves on a complicated energy landscape through a number of reaction pathways, such as oxygen rebound to form R-OH, rebound to a non-oxygen atom affording R-X, electron transfer of the incipient radical to yield a carbocation, R⁺, desaturation to form olefins, and radical cage escape. These various flavors of the rebound process, often in competition with each other, give rise to the wide range of C-H functionalization reactions performed by iron-containing oxygenases. In this review, we first recount the history of radical rebound mechanisms, their general features, and key intermediates involved. We will discuss in detail the factors that affect the behavior of the initial caged radical pair and the lifetimes of the incipient substrate radicals. Several representative examples of enzymatic C-H transformations are selected to illustrate how the behaviors of the radical pair [Fe^n-1-OH ·R] determine the eventual reaction outcome. Finally, we discuss the powerful potential of "radical rebound" processes as a general paradigm for developing novel C-H functionalization reactions with synthetic, biomimetic catalysts. We envision that new chemistry will continue to arise by bridging enzymatic "radical rebound" with synthetic organic chemistry.

Go it alone: four-electron oxidations by mononuclear non-heme iron enzymes

This review discusses the current mechanistic understanding of a group of mononuclear non-heme iron-dependent enzymes that catalyze four-electron oxidation of their organic substrates without the use of any cofactors or cosubstrates. One set of enzymes acts on α-ketoacid-containing substrates, coupling decarboxylation to oxygen activation. This group includes 4-hydroxyphenylpyruvate dioxygenase, 4-hydroxymandelate synthase, and CloR involved in clorobiocin biosynthesis. A second set of enzymes acts on substrates containing a thiol group that coordinates to the iron. This group is comprised of isopenicillin N synthase, thiol dioxygenases, and enzymes involved in the biosynthesis of ergothioneine and ovothiol. The final group of enzymes includes HEPD and MPnS that both carry out the oxidative cleavage of the carbon-carbon bond of 2-hydroxyethylphosphonate but generate different products. Commonalities amongst many of these enzymes are discussed and include the initial substrate oxidation by a ferric-superoxo-intermediate and a second oxidation by a ferryl species.

Artificial Force Induced Reaction Method for Systematic Determination of Complex Reaction Mechanisms

Nowadays, computational studies are very important for the elucidation of reaction mechanisms and selectivity of complex reactions. However, traditional computational methods usually require an estimated reaction path, mainly driven by limited experimental implications, intuition, and assumptions of stationary points. However, the artificial force induced reaction (AFIR) method in the global reaction route mapping (GRRM) strategy can be used for unbiased and automatic reaction path searches for complex reactions. In this account, we highlight applications of the AFIR method to a variety of reactions (organic, organometallic, enzymatic, and photochemical) of complex molecular systems. In addition, the AFIR method has been successfully used to rationalise the origin of stereo- and regioselectivity. The AFIR method can be applied from small to large molecular systems, and will be a very useful tool for the study of complex molecular problems in many areas of chemistry, biology, and material sciences.

Spectroscopic Evidence for the Two C-H-Cleaving Intermediates of Aspergillus nidulans Isopenicillin N Synthase

The enzyme isopenicillin N synthase (IPNS) installs the β-lactam and thiazolidine rings of the penicillin core into the linear tripeptide l-δ-aminoadipoyl-l-Cys-d-Val (ACV) on the pathways to a number of important antibacterial drugs. A classic set of enzymological and crystallographic studies by Baldwin and co-workers established that this overall four-electron oxidation occurs by a sequence of two oxidative cyclizations, with the β-lactam ring being installed first and the thiazolidine ring second. Each phase requires cleavage of an aliphatic C-H bond of the substrate: the pro-S-CCys,β-H bond for closure of the β-lactam ring, and the CVal,β-H bond for installation of the thiazolidine ring. IPNS uses a mononuclear non-heme-iron(II) cofactor and dioxygen as cosubstrate to cleave these C-H bonds and direct the ring closures. Despite the intense scrutiny to which the enzyme has been subjected, the identities of the oxidized iron intermediates that cleave the C-H bonds have been addressed only computationally; no experimental insight into their geometric or electronic structures has been reported. In this work, we have employed a combination of transient-state-kinetic and spectroscopic methods, together with the specifically deuterium-labeled substrates, A[d2-C]V and AC[d8-V], to identify both C-H-cleaving intermediates. The results show that they are high-spin Fe(III)-superoxo and high-spin Fe(IV)-oxo complexes, respectively, in agreement with published mechanistic proposals derived computationally from Baldwin's founding work.

Radicalization and Radical Catalysis of Biomass Sugars: Insights from First-principles Studies

Ab initio and density functional calculations are conducted to investigate the radicalization processes and radical catalysis of biomass sugars. Structural alterations due to radicalization generally focus on the radicalized sites, and radicalization affects H-bonds in D-fructofuranose more than in D-glucopyranose, potentially with outcome of new H-bonds. Performances of different functionals and basis sets are evaluated for all radicalization processes, and enthalpy changes and Gibbs free energies for these processes are presented with high accuracy, which can be referenced for subsequent experimental and theoretical studies. It shows that radicalization can be utilized for direct transformation of biomass sugars, and for each sugar, C rather than O sites are always preferred for radicalization, thus suggesting the possibility to activate C-H bonds of biomass sugars. Radical catalysis is further combined with Brønsted acids, and it clearly states that functionalization fundamentally regulates the catalytic effects of biomass sugars. In presence of explicit water molecules, functionalization significantly affects the activation barriers and reaction energies of protonation rather than dehydration steps. Tertiary butyl and phenyl groups with large steric hindrances or hydroxyl and amino groups resulting in high stabilities for protonation products drive the protonation steps to occur facilely at ambient conditions.

DFT-ONIOM study of the dopamine–β-CD complex: NBO and AIM analysis

In this work, we conducted a systematic search of the minimum energy of a dopamine–β-cyclodextrin complex via different ONIOM approaches using both mixed (DFT-HF) and nonmixed (DFT-DFT) combinations. Different density functionals were employed: B3LYP, MPW1PW91, M05-2X, M06-2X, and ωB97X-D. Two different basis sets were used on the dopamine–β-cyclodextrin complex; a lower basis set (3-21G*) is used on β-cyclodextrin and a higher basis set (6-31G(d)) is used on dopamine. To fulfill this task, complexation and deformation energies were determined. The relative performance of these functionals was compared with that obtained with the DFT method at the M06-2X/6-31G(d) level that is chosen as a reference method. Although we did not clearly establish an assessment of the relative performance of these density functionals, the efficacy of the ONIOM combination (DFT-HF) is shown compared with a nonmixed ONIOM combination (DFT-DFT). The intermolecular hydrogen bond interactions of the complex obtained with the full M06-2X/6-31G(d) have been analyzed with the atoms-in-molecules and natural bond rrbital methodologies.

Molecular Dynamics Simulation of Nitrobenzene Dioxygenase Using AMBER Force Field

Molecular dynamics simulation of the oxygenase component of nitrobenzene dioxygenase (NBDO) system, a member of the naphthalene family of Rieske nonheme iron dioxygenases, has been carried out using the AMBER force field combined with a new set of parameters for the description of the mononuclear nonheme iron center and iron-sulfur Rieske cluster. Simulation results provide information on the structure and dynamics of nitrobenzene dioxygenase in an aqueous environment and shed light on specific interactions that occur in its catalytic center. The results suggest that the architecture of the active site is stabilized by key hydrogen bonds, and Asn258 positions the substrate for oxidation. Analysis of protein-water interactions reveal the presence of a network of solvent molecules at the entrance to the active site, which could be of potential catalytic importance.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

A DFT and ONIOM study of C–H hydroxylation catalyzed by nitrobenzene 1,2-dioxygenase

The DFT and ONIOM calculations show that C–H hydroxylation by nitrobenzene 1,2-dioxygenase proceeds through a HO–Fe^VO intermediate.

Standard role for a conserved aspartate or more direct involvement in deglycosylation? An ONIOM and MD investigation of adenine-DNA glycosylase

8-Oxoguanine (OG) is one of the most frequently occurring forms of DNA damage and is particularly deleterious since it forms a stable Hoogsteen base pair with adenine (A). The repair of an OG:A mispair is initiated by adenine-DNA glycosylase (MutY), which hydrolyzes the sugar-nucleobase bond of the adenine residue before the lesion is processed by other proteins. MutY has been proposed to use a two-part chemical step involving protonation of the adenine nucleobase, followed by SN1 hydrolysis of the glycosidic bond. However, differences between a recent (fluorine recognition complex, denoted as the FLRC) crystal structure and the structure on which most mechanistic conclusions have been based to date (namely, the lesion recognition complex or LRC) raise questions regarding the mechanism used by MutY and the discrete role of various active-site residues. The present work uses both molecular dynamics (MD) and quantum mechanical (ONIOM) models to compare the active-site conformational dynamics in the two crystal structures, which suggests that only the understudied FLRC leads to a catalytically competent reactant. Indeed, all previous computational studies on MutY have been initiated from the LRC structure. Subsequently, for the first time, various mechanisms are examined with detailed ONIOM(M06-2X:PM6) reaction potential energy surfaces (PES) based on the FLRC structure, which significantly extends the mechanistic picture. Specifically, our work reveals that the reaction proceeds through a different route than the commonly accepted mechanism and the catalytic function of various active-site residues (Geobacillus stearothermophilus numbering). Specifically, contrary to proposals based on the LRC, E43 is determined to solely be involved in the initial adenine protonation step and not the deglycosylation reaction as the general base. Additionally, a novel catalytic role is proposed for Y126, whereby this residue plays a significant role in stabilizing the highly charged active site, primarily through interactions with E43. More importantly, D144 is found to explicitly catalyze the nucleobase dissociation step through partial nucleophilic attack. Although this is a more direct role than previously proposed for any other DNA glycosylase, comparison to previous work on other glycosylases justifies the larger contribution in the case of MutY and allows us to propose a unified role for the conserved Asp/Glu in the DNA glycosylases, as well as other enzymes that catalyze nucleotide deglycosylation reactions.

Model iron-oxo species and the oxidation of imidazole: insights into the mechanism of OvoA and EgtB?

A density functional theory cluster and first-principles quantum and statistical mechanics approach have been used to investigate the ability of iron-oxygen intermediates to oxidize a histidine cosubstrate, which may then allow for the possible formation of 2- and 5-histidylcysteine sulfoxide, respectively. Namely, the ability of ferric superoxo (Fe(III)O(2)(•-)), Fe(IV)═O, and ferrous peroxysulfur (Fe(III)OOS) complexes to oxidize the imidazole of histidine via an electron transfer (ET) or a proton-coupled electron transfer (PCET) was considered. While the high-valent mononuclear Fe(IV)═O species is generally considered the ultimate biooxidant, the free energies for its reduction (via ET or PCET) suggest that it is unable to directly oxidize histidine's imidazole. Instead, only the ferrous peroxysulfur complexes are sufficiently powerful enough oxidants to generate a histidyl-derived radical via a PCET process. Furthermore, while this process preferably forms a HisN(δ)(-H)(•) radical, several such oxidants are also suggested to be capable of generating the higher-energy HisC(δ)(-H)(•) and HisC(ε)(-H)(•) radicals. Importantly, the present results suggest that formation of the sulfoxide-containing products (seen in both OvoA and EgtB) is a consequence of the reduction of a powerful Fe(III)OOS oxidant via a PCET.

Combined QM/MM (ONIOM) and QSAR approach to the study of complex formation of matrix metalloproteinase‑9 with a series of biphenylsulfonamides−LERE-QSAR analysis (V)

We previously proposed a novel QSAR (quantitative structure-activity relationship) procedure called LERE (linear expression by representative energy terms)-QSAR involving molecular calculations such as an ab initio fragment molecular orbital ones. In the present work, we applied LERE-QSAR to complex formation of matrix metalloproteinase-9 (MMP-9) with a series of substituted biphenylsulfonamides. The results shows that the overall free-energy change accompanying complex formation is due to predominantly the contribution from the electrostatic interaction with the zinc atom in the active site of MMP-9. Carbonic anhydrase (CA) belongs to the zinc-containing protease family. In contrast to the current case of MMP-9, the overall free-energy change during complex formation of CA with a series of benzenesulfonamides is due to the contributions from the solvation and dissociation free-energy changes, as previously reported. Comparison of the two sets of results indicates quantitative differences in the relative contributions of free-energy components to the overall free-energy change between the two data sets, corresponding with those in the respective classical QSAR equations. The LERE-QSAR procedure was demonstrated to quantitatively reveal differences in the binding mechanisms between the two cases involving similar but different zinc-containing proteins at the electronic and atomic levels.