Search

One of the most important reactions in fullerene chemistry is the Diels-Alder (DA) reaction. In two previous experimental studies, the DA cycloaddition reactions of cyclopentadiene (Cp) and 1,2,3,4,5-pentamethylcyclopentadiene (Cp*) with La@ C2 v-C82 were investigated. The attack of Cp was proposed to occur on bond 19, whereas that of Cp* was confirmed by X-ray analysis to be over bond o. Moreover, the stabilities of the Cp and Cp* adducts were found to be significantly different, that is, the decomposition of La@ C2 v-C82Cp was one order of magnitude faster than that of La@ C2 v-C82Cp*. Herein, we computationally analyze these DA cycloadditions with two main goals: First, to compute the thermodynamics and kinetics of the cycloadditions of Cp and Cp* to different bonds of La@ C2 v-C82 to assess and compare the regioselectivity of these two reactions. Second, to understand the origin of the different thermal stabilities of the La@C82Cp and La@C82Cp* adducts. Our results show that the regioselectivity of the two DA cycloadditions is the same, with preferred attack on bond o. This result corrects the previous assumption of the regioselectivity of the Cp attack that was made based only on the shape of the La@C82 singly occupied molecular orbital. In addition, we show that the higher stability of the La@C82Cp* adduct is not due to the electronic effects of the methyl groups on the Cp ring, as previously suggested, but to higher long-range dispersion interactions in the Cp* case, which enhance the stabilization of the reactant complex, transition state, and products with respect to the separated reactants. This stabilization for the La@C82Cp* case decreases the Gibbs reaction energy, thus allowing competition between the direct and retro reactions and making dissociation more difficult. [ABSTRACT FROM AUTHOR]

The metal-directed supramolecular synthetic approach has paved the way for the development of functional nanosized molecules. In this work, we report the preparation of the new nanocapsule 3⋅(CF3SO3)8 with a A4B2 tetragonal prismatic geometry, where A corresponds to the dipalladium hexaazamacrocyclic complex Pd-1, and B corresponds to the tetraanionic form of palladium 5,10,15,20-tetrakis(4-carboxyphenyl)porphyrin ( 2). The large void space of the inner cavity and the supramolecular affinity for guest molecules towards porphyrin-based hosts converts this nanoscale molecular 3D structure into a good candidate for host-guest chemistry. The interaction between this nanocage and different guest molecules has been studied by means of NMR, UV/Vis, ESI-MS, and DOSY experiments, from which highly selective molecular recognition has been found for anionic, planar-shaped π guests with association constants ( Ka) higher than 109 M−1, in front of non-interacting aromatic neutral or cationic substrates. DFT theoretical calculations provided insights to further understand this strong interaction. Nanocage 3⋅(CF3SO3)8 can not only strongly host one single molecule of M(dithiolene)2 complexes (M=Au, Pt, Pd, and Ni), but also can finely tune their optical and redox properties. The very simple synthesis of both the supramolecular cage and the building blocks represents a step forward for the development of polyfunctional supramolecular nanovessels, which offer multiple applications as sensors or nanoreactors. [ABSTRACT FROM AUTHOR]

The chemical functionalization of endohedral (metallo)fullerenes has become a main focus of research in the last few years. It has been found that the reactivity of endohedral (metallo)fullerenes may be quite different from that of the empty fullerenes. Encapsulated species have an enormous influence on the thermodynamics, kinetics, and regiochemistry of the exohedral addition reactions undergone by these species. A detailed understanding of the changes in chemical reactivity due to incarceration of atoms or clusters of atoms is essential to assist the synthesis of new functionalized endohedral fullerenes with specific properties. Herein, we report the study of the Diels-Alder cycloaddition between 1,3-butadiene and all nonequivalent bonds of the Ti2C2@ D3 h-C78 metallic carbide endohedral metallofullerene (EMF) at the BP86/TZP//BP86/DZP level of theory. The results obtained are compared with those found by some of us at the same level of theory for the D3 h-C78 free cage and the M3N@ D3 h-C78 (M=Sc and Y) metallic nitride EMFs. It is found that the free cage is more reactive than the Ti2C2@ D3 h-C78 EMF and this, in turn, has a higher reactivity than M3N@ D3 h-C78. The results indicate that, for Ti2C2@ D3 h-C78, the corannulene-type [5, 6] bonds c and f, and the type B [6, 6] bond 3 are those thermodynamically and kinetically preferred. In contrast, the D3 h-C78 free cage has a preference for addition to the [6, 6] 1 and 6 bonds and the [5, 6] b bond, whereas M3N@ D3 h-C78 favors additions to the [6, 6] 6 (M=Sc) and [5, 6] d (M=Y) bonds. The reasons for the regioselectivity found in Ti2C2@ D3 h-C78 are discussed. [ABSTRACT FROM AUTHOR]

In this work we show that the regioselectivity of the Diels–Alder, 1,3-dipolar, and carbene cycloadditions to C60 changes from the usual [6,6] addition in neutral species to the [5,6] attack when the number of electrons added to the fullerenic cage increases. Changes in the aromaticity of the five- and six-membered rings of C60 during the reduction process provide a rationale to understand this regioselectivity change. [ABSTRACT FROM AUTHOR]

Methods for regioselective N‐methylation and ‐alkylation of unsaturated heterocycles with "off the shelf" reagents are highly sought‐after. This reaction could drastically simplify synthesis of privileged bioactive molecules. Here we report engineered and natural methyltransferases for challenging N‐(m)ethylation of heterocycles, including benzimidazoles, benzotriazoles, imidazoles and indazoles. The reactions are performed through a cyclic enzyme cascade that consists of two methyltransferases using only iodoalkanes or methyl tosylate as simple reagents. This method enables the selective synthesis of important molecules that are otherwise difficult to access, proceeds with high regioselectivity (r.r. up to >99 %), yield (up to 99 %), on a preparative scale, and with nearly equimolar concentrations of simple starting materials. [ABSTRACT FROM AUTHOR]

The stereochemical outcome of cis--trans isomerization of optically pure [60], [70], and endohedral H2O@C60 fulleropyrrolidines reveals that the electronic nature of substituents, fullerene size, and surprisingly the incarcerated water molecule plays a crucial role in this rearrangement process. Theoretical DFT calculations are in very good agreement with the experimental findings. On the basis of the experimental results and computational calculations, a plausible reaction mechanism involving the hydrogen-bonding assistance of the inner water molecule in the carbanion stabilization of endofullerene is proposed. [ABSTRACT FROM AUTHOR]

The mechanism by which molecular oxygen is activated by the organic cofactor pyridoxal phosphate (PLP) for oxidation reactions remains poorly understood. Recent work has identified arginine oxidases that catalyze desaturation or hydroxylation reactions. Here, we investigate a desaturase from the Pseudoalteromonas luteoviolacea indolmycin pathway. Our work, combining X-ray crystallographic, biochemical, spectroscopic, and computational studies, supports a shared mechanism with arginine hydroxylases, involving two rounds of single-electron transfer to oxygen and superoxide rebound at the 4′ carbon of the PLP cofactor. The precise positioning of a water molecule in the active site is proposed to control the final reaction outcome. This proposed mechanism provides a unified framework to understand how oxygen can be activated by PLP-dependent enzymes for oxidation of arginine and elucidates a shared mechanistic pathway and intertwined evolutionary history for arginine desaturases and hydroxylases. [ABSTRACT FROM AUTHOR]

Cytochromes P450 (P450s) are nature's catalysts of choice for performing demanding and physiologically vital oxidation reactions. Biochemical characterization of these enzymes over the past decades has provided detailed mechanistic insight and highlighted the diversity of substrates P450s accommodate and the spectrum of oxidative transformations they catalyze. Previously, we discovered that the bacterial P450 MycCI from the mycinamicin biosynthetic pathway in Micromonospora griseorubida possesses an unusually broad substrate scope, whereas the homologous P450 from tylosin-producing Streptomyces fradiae (TylHI) exhibits a high degree of specificity for its native substrate. Here, using biochemical, structural, and computational approaches, we aimed to understand the molecular basis for the disparate reactivity profiles of these two P450s. Turnover and equilibrium binding experiments with substrate analogs revealed that TylHI strictly prefers 16-membered ring macrolides bearing the deoxyamino sugar mycaminose. To help rationalize these results, we solved the X-ray crystal structure of TylHI in complex with its native substrate at 1.99-Å resolution and assayed several site-directed mutants. We also conducted molecular dynamics simulations of TylHI and MycCI and biochemically characterized a third P450 homolog from the chalcomycin biosynthetic pathway in Streptomyces bikiniensis. These studies provided a basis for constructing P450 chimeras to gain further insight into the features dictating the differences in reaction profile among these structurally and functionally related enzymes, ultimately unveiling the central roles of key loop regions in influencing substrate binding and turnover. Our work highlights the complex nature of P450/substrate interactions and raises interesting questions regarding the evolution of functional diversity among biosynthetic enzymes. [ABSTRACT FROM AUTHOR]

Microbial conversion of aromatic compounds is an emerging and promising strategy for valorization of the plant biopolymer lignin. A critical and often rate-limiting reaction in aromatic catabolism is O-aryl-demethylation of the abundant aromatic methoxy groups in lignin to form diols, which enables subsequent oxidative aromatic ring-opening. Recently, a cytochrome P450 system, GcoAB, was discovered to demethylate guaiacol (2-methoxyphenol), which can be produced from coniferyl alcohol-derived lignin, to form catechol. However, native GcoAB has minimal ability to demethylate syringol (2,6-dimethoxyphenol), the analogous compound that can be produced from sinapyl alcohol-derived lignin. Despite the abundance of sinapyl alcohol-based lignin in plants, no pathway for syringol catabolism has been reported to date. Here we used structureguided protein engineering to enable microbial syringol utilization with GcoAB. Specifically, a phenylalanine residue (GcoA-F169) interferes with the binding of syringol in the active site, and on mutation to smaller amino acids, efficient syringol O-demethylation is achieved. Crystallography indicates that syringol adopts a productive binding pose in the variant, which molecular dynamics simulations trace to the elimination of steric clash between the highly flexible side chain of GcoA-F169 and the additional methoxy group of syringol. Finally, we demonstrate in vivo syringol turnover in Pseudomonas putida KT2440 with the GcoA-F169A variant. Taken together, our findings highlight the significant potential and plasticity of cytochrome P450 aromatic O-demethylases in the biological conversion of lignin-derived aromatic compounds. [ABSTRACT FROM AUTHOR]

The first chemical modification on the brand new endohedral HF@C60 is reported. In particular, the isomerization from optically pure (2S,5S)-cis-pyrrolidino[3,4:1,2][60]fullerene 2b to (2S,5R)-trans-pyrrolidino[3,4:1,2][60]fullerene 2b has been studied and compared with empty C60 (2a) and endohedral H2O@C60 (3). The comparative study shows a kinetic order for the isomerization process of H2O@C60 ＞ HF@C60 ＞ C60, thus confirming the effect of the incarcerated species on the zwitterionic intermediate stability. [ABSTRACT FROM AUTHOR]

Enzymes are renowned for their catalytic efficiency and selectivity, but many classical transformations in organic synthesis have no biocatalytic counterpart. Aldolases are prodigious C-C bond-forming enzymes, but their reactivity has only been extended past activated carbonyl electrophiles in special cases. To probe the mechanistic origins of this limitation, we use a pair of aldolases whose activity is dependent on pyridoxal phosphate. Our results reveal how aldolases are limited by kinetically favourable proton transfer with solvent, which undermines aldol addition into ketones. We show how a transaldolase can circumvent this limitation, enabling efficient addition into unactivated ketones. The resulting products are highly sought non-canonical amino acids with side chains that contain chiral tertiary alcohols. Mechanistic analysis reveals that transaldolase activity is an intrinsic feature of pyridoxal phosphate chemistry and identifies principles for extending aldolase catalysis beyond its previous limits to enable convergent, enantioselective C-C bond formation from simple starting materials., Competing Interests: Competing interests: S.K.B., A.M., T.J.D. and A.R.B. are the inventors on a provisional patent application (number 18/351,559) submitted by the University of Wisconsin related to this work. The other authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

In recent years there has been a surge in the development of methods for the synthesis of organofluorine compounds. However, enzymatic methods for C-F bond formation have been limited to nucleophilic fluoride substitution. Here, we report the incorporation of iron-catalysed radical fluorine transfer, a reaction mechanism that is not used in naturally occurring enzymes, into enzymatic catalysis for the development of biocatalytic enantioselective C( sp 3 )-F bond formation. Using this strategy, we repurposed ( S )-2-hydroxypropylphosphonate epoxidase from Streptomyces viridochromogenes ( Sv HppE) to catalyse an N -fluoroamide directed C( sp 3 )-H fluorination. Directed evolution has enabled Sv HppE to be optimized, forming diverse chiral benzylic fluoride products with turnover numbers of up to 180 and with excellent enantiocontrol (up to 94% e.e.). Mechanistic investigations showed that the N-F bond activation is the rate-determining step, and the strong preference for fluorination in the presence of excess NaN 3 can be attributed to the spatial proximity of the carbon-centered radical to the iron-bound fluoride., Competing Interests: Competing interests A provisional patent application (No. PCT/US2023/022431) covering enantioselective biocatalytic C–F bond formation has been filed through the Johns Hopkins University with Q.Z., X.H., Z.C., and J.R. as inventors. Authors J.S., N.J., E.Y., X.C., Y.Y., and M.G.B. declare no competing interests.

S-Adenosyl-l-methionine (SAM) is an important cosubstrate in various biochemical processes, including selective methyl transfer reactions. Simple methods for the (re)generation of SAM analogs could expand the chemistry accessible with SAM-dependent transferases and go beyond methylation reactions. Here we present an efficient enzyme engineering strategy to synthesize different SAM analogs from "off-the-shelf" iodoalkanes through enzymatic alkylation of S-adenosyl-l-homocysteine (SAH). This was achieved by mutating multiple hydrophobic and structurally dynamic amino acids simultaneously. Combinatorial mutagenesis was guided by the natural amino acid diversity and generated a highly functional mutant library. This approach increased the speed as well as the scale of enzyme engineering by providing a panel of optimized enzymes with orders of magnitude higher activities for multiple substrates in just one round of enzyme engineering. The optimized enzymes exhibit catalytic efficiencies up to 31 M -1  s -1 , convert various iodoalkanes, including substrates bearing cyclopropyl or aromatic moieties, and catalyze S-alkylation of SAH with very high stereoselectivities (>99 % de). We further report a high throughput chromatographic screening system for reliable and rapid SAM analog analysis. We believe that the methods and enzymes described herein will further advance the field of selective biocatalytic alkylation chemistry by enabling SAM analog regeneration with "off-the-shelf" reagents., (© 2024 The Authors. ChemBioChem published by Wiley-VCH GmbH.)

P450 enzymes naturally perform selective hydroxylations and epoxidations of unfunctionalized hydrocarbon substrates, among other reactions. The adaptation of P450 enzymes to a particular oxidative reaction involving alkenes is of great interest for the design of new synthetically useful biocatalysts. However, the mechanism that these enzymes utilize to precisely modulate the chemoselectivity and distinguishing between competing alkene double bond epoxidations and allylic C-H hydroxylations is sometimes not clear, which hampers the rational design of specific biocatalysts. In a previous work, a P450 from Labrenzia aggregata (P450 LA1 ) was engineered in the laboratory using directed evolution to catalyze the direct oxidation of trans-β-methylstyrene to phenylacetone. The final variant, KS, was able to overcome the intrinsic preference for alkene epoxidation to directly generate a ketone product via the formation of a highly reactive carbocation intermediate. Here, additional library screening along this evolutionary lineage permitted to serendipitously detect a mutation that overcomes epoxidation and carbonyl formation by exhibiting a large selectivity of 94 % towards allylic C-H hydroxylation. A multiscalar computational methodology was applied to reveal the molecular basis towards this hydroxylation preference. Enzyme modelling suggests that introduction of a bulky substitution dramatically changes the accessible conformations of the substrate in the active site, thus modifying the enzymatic selectivity towards terminal hydroxylation and avoiding the competing epoxidation pathway, which is sterically hindered., (© 2024 The Authors. ChemBioChem published by Wiley-VCH GmbH.)

Catalysis with engineered enzymes has provided more efficient routes for the production of active pharmaceutical agents. However, the potential of biocatalysis to assist in early-stage drug discovery campaigns remains largely untapped. In this study, we have developed a biocatalytic strategy for the construction of sp 3 -rich polycyclic compounds via the intramolecular cyclopropanation of benzothiophenes and related heterocycles. Two carbene transferases with complementary regioisomer selectivity were evolved to catalyse the stereoselective cyclization of benzothiophene substrates bearing diazo ester groups at the C2 or C3 position of the heterocycle. The detailed mechanisms of these reactions were elucidated by a combination of crystallographic and computational analyses. Leveraging these insights, the substrate scope of one of the biocatalysts could be expanded to include previously unreactive substrates, highlighting the value of integrating evolutionary and rational strategies to develop enzymes for new-to-nature transformations. The molecular scaffolds accessed here feature a combination of three-dimensional and stereochemical complexity with 'rule-of-three' properties, which should make them highly valuable for fragment-based drug discovery campaigns., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Alkylidene cyclopropanes (ACPs) are valuable synthetic intermediates because of their constrained structure and opportunities for further diversification. Although routes to ACPs are known, preparations of ACPs with control of both the configuration of the cyclopropyl ( R vs S ) group and the geometry of the alkene ( E vs Z ) are unknown. We describe enzymatic cyclopropanation of allenes with ethyl diazoacetate (EDA) catalyzed by an iridium-containing cytochrome (Ir(Me)-CYP119) that controls both stereochemical elements. Two mutants of Ir(Me)-CYP119 identified by 6-codon (6c, VILAFG) saturation mutagenesis catalyze the formation of ( E )-ACPs with -93% to >99% ee and >99:1 E / Z ratio with just three rounds of 96 mutants. By four additional rounds of mutagenesis, an enzyme variant was identified that forms ( Z )-ACPs with up to 94% ee and a 28:72 E/Z ratio. Computational studies show that the orientation of the carbene unit dictated by the mutated positions accounts for the stereoselectivity.

We report a computationally driven approach to access enantiodivergent enzymatic carbene N-H insertions catalyzed by P411 enzymes. Computational modeling was employed to rationally guide engineering efforts to control the accessible conformations of a key lactone-carbene (LAC) intermediate in the enzyme active site by installing a new H-bond anchoring point. This H-bonding interaction controls the relative orientation of the reactive carbene intermediate, orienting it for an enantioselective N-nucleophilic attack by the amine substrate. By combining MD simulations and site-saturation mutagenesis and screening targeted to only two key residues, we were able to reverse the stereoselectivity of previously engineered S-selective P411 enzymes. The resulting variant, L5_FL-B3, accepts a broad scope of amine substrates for N-H insertion with excellent yields (up to >99 %), high efficiency (up to 12 300 TTN), and good enantiocontrol (up to 7 : 93 er)., (© 2023 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Trifluoromethyl-substituted cyclopropanes (CF 3 -CPAs) constitute an important class of compounds for drug discovery. While several methods have been developed for synthesis of trans-CF 3 -CPAs, stereoselective production of corresponding cis-diastereomers remains a formidable challenge. We report a biocatalyst for diastereo- and enantio-selective synthesis of cis-CF 3 -CPAs with activity on a variety of alkenes. We found that an engineered protoglobin from Aeropyrnum pernix (ApePgb) can catalyze this unusual reaction at preparative scale with low-to-excellent yield (6-55 %) and enantioselectivity (17-99 % ee), depending on the substrate. Computational studies revealed that the steric environment in the active site of the protoglobin forced iron-carbenoid and substrates to adopt a pro-cis near-attack conformation. This work demonstrates the capability of enzyme catalysts to tackle challenging chemistry problems and provides a powerful means to expand the structural diversity of CF 3 -CPAs for drug discovery., (© 2022 Wiley-VCH GmbH.)

The aerobic oxidation of alkenes to carbonyls is an important and challenging transformation in synthesis. Recently, a new P450-based enzyme (aMOx) has been evolved in the laboratory to directly oxidize styrenes to their corresponding aldehydes with high activity and selectivity. The enzyme utilizes a heme-based, high-valent iron-oxo species as a catalytic oxidant that normally epoxidizes alkenes, similar to other catalysts. How the evolved aMOx enzyme suppresses the commonly preferred epoxidation and catalyzes direct carbonyl formation is currently not well understood. Here, we combine computational modelling together with mechanistic experiments to study the reaction mechanism and unravel the molecular basis behind the selectivity achieved by aMOx. Our results describe that although both pathways are energetically accessible diverging from a common covalent radical intermediate, intrinsic dynamic effects determine the strong preference for epoxidation. We discovered that aMOx overrides these intrinsic preferences by controlling the accessible conformations of the covalent radical intermediate. This disfavors epoxidation and facilitates the formation of a carbocation intermediate that generates the aldehyde product through a fast 1,2-hydride migration. Electrostatic preorganization of the enzyme active site also contributes to the stabilization of the carbocation intermediate. Computations predicted that the hydride migration is stereoselective due to the enzymatic conformational control over the intermediate species. These predictions were corroborated by experiments using deuterated styrene substrates, which proved that the hydride migration is cis - and enantioselective. Our results demonstrate that directed evolution tailored a highly specific active site that imposes strong steric control over key fleeting biocatalytic intermediates, which is essential for accessing the carbonyl forming pathway and preventing competing epoxidation.

New-to-nature radical biocatalysis has recently emerged as a powerful strategy to tame fleeting open-shell intermediates for stereoselective transformations. In 2021, we introduced a novel metalloredox biocatalysis strategy that leverages the innate redox properties of the heme cofactor of P450 enzymes, furnishing new-to-nature atom-transfer radical cyclases (ATRCases) with excellent activity and stereoselectivity. Herein, we report a combined computational and experimental study to shed light on the mechanism and origins of enantioselectivity for this system. Molecular dynamics and quantum mechanics/molecular mechanics (QM/MM) calculations revealed an unexpected role of the key beneficial mutation I263Q. The glutamine residue serves as an essential hydrogen bond donor that engages with the carbonyl moiety of the substrate to promote bromine atom abstraction and enhance the enantioselectivity of radical cyclization. Therefore, the evolved ATRCase is a bifunctional biocatalyst, wherein the heme cofactor enables atom-transfer radical biocatalysis, while the hydrogen bond donor residue further enhances the activity and enantioselectivity. Unlike many enzymatic stereocontrol rationales based on a rigid substrate binding model, our computations demonstrate a high degree of rotational flexibility of the allyl moiety in an enzyme-substrate complex and succeeding intermediates. Therefore, the enantioselectivity is controlled by the radical cyclization transition states rather than the substrate orientation in ground-state complexes in the preceding steps. During radical cyclization, anchoring effects of the Q263 residue and steric interactions with the heme cofactor concurrently control the π-facial selectivity, allowing for highly enantioselective C-C bond formation. Our computational findings are corroborated by experiments with ATRCase mutants generated from site-directed mutagenesis.

Protein-ligand binding processes often involve changes in protonation states that can be key to recognize and orient the ligand in the binding site. The pathways through which (bio)molecules interplay to attain productively bound complexes are intricate and involve a series of interconnected intermediate and transition states. Molecular dynamics (MD) simulations and enhanced sampling techniques are commonly used to characterize the spontaneous binding of a ligand to its receptor. However, the effect of protonation state changes of in-pathway residues in spontaneous binding MD simulations remained mostly unexplored. Here, we used molecular dynamics simulations to reconstruct the trypsin-benzamidine binding pathway considering different protonation states of His57. This residue is part of the trypsin catalytic triad and is located more than 10 Å away from Asp189, which is responsible for benzamidine binding in the trypsin S1 pocket. Our MD simulations showed that the binding pathways that benzamidine follow to target the S1 binding site are critically dependent on the His57 protonation state. Binding of benzamidine frequently occurs when His57 is protonated in the delta nitrogen while the binding process is significantly less frequent when His57 is positively charged. Constant-pH MD simulations retrieved the equilibrium populations of His57 protonation states at trypsin active pH offering a clearer picture of benzamidine recognition and binding. These results indicate that properly accounting for protonation states of distal residues can be important in spontaneous binding MD simulations., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Girame, Garcia-Borràs and Feixas.)

Propargyl amines are versatile synthetic intermediates with numerous applications in the pharmaceutical industry. An attractive strategy for efficient preparation of these compounds is nitrene propargylic C(sp 3 )-H insertion. However, achieving this reaction with good chemo-, regio-, and enantioselective control has proven to be challenging. Here, we report an enzymatic platform for the enantioselective propargylic amination of alkynes using a hydroxylamine derivative as the nitrene precursor. Cytochrome P450 variant PA-G8 catalyzing this transformation was identified after eight rounds of directed evolution. A variety of 1-aryl-2-alkyl alkynes are accepted by PA-G8 , including those bearing heteroaromatic rings. This biocatalytic process is efficient and selective (up to 2610 total turnover number (TTN) and 96% ee ) and can be performed on preparative scale.

Chiral amines can be made by insertion of a carbene into an N-H bond using two-catalyst systems that combine a transition metal-based carbene-transfer catalyst and a chiral proton-transfer catalyst to enforce stereocontrol. Haem proteins can effect carbene N-H insertion, but asymmetric protonation in an active site replete with proton sources is challenging. Here we describe engineered cytochrome P450 enzymes that catalyse carbene N-H insertion to prepare biologically relevant α-amino lactones with high activity and enantioselectivity (up to 32,100 total turnovers, >99% yield and 98% e.e.). These enzymes serve as dual-function catalysts, inducing carbene transfer and promoting the subsequent proton transfer with excellent stereoselectivity in a single active site. Computational studies uncover the detailed mechanism of this new-to-nature enzymatic reaction and explain how active-site residues accelerate this transformation and provide stereocontrol., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Iterative P450 enzymes are powerful biocatalysts for selective late-stage C-H oxidation of complex natural product scaffolds. These enzymes represent useful tools for selectivity and cascade reactions, facilitating direct access to core structure diversification. Recently, we reported the structure of the multifunctional bacterial P450 TamI and elucidated the molecular basis of its substrate binding and strict reaction sequence at distinct carbon atoms of the substrate. Here, we report the design and characterization of a toolbox of TamI biocatalysts, generated by mutations at Leu101, Leu244, and/or Leu295, that alter the native selectivity, step sequence, and number of reactions catalyzed, including the engineering of a variant capable of catalyzing a four-step oxidative cascade without the assistance of the flavoprotein and oxidative partner TamL. The tuned enzymes override inherent substrate reactivity, enabling catalyst-controlled C-H functionalization and alkene epoxidation of the tetramic acid-containing natural product tirandamycin. Five bioactive tirandamycin derivatives (6-10) were generated through TamI-mediated enzymatic synthesis. Quantum mechanics calculations and MD simulations provide important insights into the basis of altered selectivity and underlying biocatalytic mechanisms for enhanced continuous oxidation of the iterative P450 TamI., Competing Interests: The authors declare no competing financial interest.

A cytochrome c heme protein was recently engineered to catalyze the formation of carbon-silicon bonds via carbene insertion into Si-H bonds, a reaction that was not previously known to be catalyzed by a protein. High chemoselectivity toward C-Si bond formation over competing C-N bond formation was achieved, although this trait was not screened for during directed evolution. Using computational and experimental tools, we now establish that activity and chemoselectivity are modulated by conformational dynamics of a protein loop that covers the substrate access to the iron-carbene active species. Mutagenesis of residues computationally predicted to control the loop conformation altered the protein's chemoselectivity from preferred silylation to preferred amination of a substrate containing both N-H and Si-H functionalities. We demonstrate that information on protein structure and conformational dynamics, combined with knowledge of mechanism, leads to understanding of how non-natural and selective chemical transformations can be introduced into the biological world.

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Biocatalysis offers an expanding and powerful strategy to construct and diversify complex molecules by C─H bond functionalization. Due to their high selectivity, enzymes have become an essential tool for C─H bond functionalization and offer complementary reactivity to small-molecule catalysts. Hemoproteins, particularly cytochromes P450, have proven effective for selective oxidation of unactivated C─H bonds. Previously, we reported the in vitro characterization of an oxidative tailoring cascade in which TamI, a multifunctional P450 functions co-dependently with the TamL flavoprotein to catalyze regio- and stereoselective hydroxylations and epoxidation to yield tirandamycin A and tirandamycin B. TamI follows a defined order including 1) C10 hydroxylation, 2) C11/C12 epoxidation, and 3) C18 hydroxylation. Here we present a structural, biochemical, and computational investigation of TamI to understand the molecular basis of its substrate binding, diverse reactivity, and specific reaction sequence. The crystal structure of TamI in complex with tirandamycin C together with molecular dynamics simulations and targeted mutagenesis suggest that hydrophobic interactions with the polyene chain of its natural substrate are critical for molecular recognition. QM calculations and molecular dynamics simulations of TamI with variant substrates provided detailed information on the molecular basis of sequential reactivity, and pattern of regio- and stereo-selectivity in catalyzing the three-step oxidative cascade., Competing Interests: Conflict of Interest: The authors declare no competing interests.

Compared to the biological world's rich chemistry for functionalizing carbon, enzymatic transformations of the heavier homologue silicon are rare. We report that a wild-type cytochrome P450 monooxygenase (P450 BM3 from Bacillus megaterium, CYP102A1) has promiscuous activity for oxidation of hydrosilanes to give silanols. Directed evolution was applied to enhance this non-native activity and create a highly efficient catalyst for selective silane oxidation under mild conditions with oxygen as the terminal oxidant. The evolved enzyme leaves C-H bonds present in the silane substrates untouched, and this biotransformation does not lead to disiloxane formation, a common problem in silanol syntheses. Computational studies reveal that catalysis proceeds through hydrogen atom abstraction followed by radical rebound, as observed in the native C-H hydroxylation mechanism of the P450 enzyme. This enzymatic silane oxidation extends nature's impressive catalytic repertoire., (© 2020 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The tris- and tetra-adducts of M 3 N@ I h - C 80 metallofullerenes were synthesized and characterized for the first time. The 1,3-dipolar cycloaddition (Prato reaction) of Y 3 N@ I h -C 80 and Gd 3 N@ I h -C 80 with an excess of N -ethylglycine and formaldehyde provided tris- and tetra-fulleropyrrolidine adducts in a regioselective manner. Purification by HPLC and analyses of the isolated peaks by NMR, MS, and vis-NIR spectra revealed that the major products were four tris- and one tetra-isomers for both Y 3 N@ I h -C 80 and Gd 3 N@ I h -C 80 . Considering the large number of possible isomers (e.g., at least 1140 isomers for the tris-adduct), the limited number of isomers obtained indicated that the reactions proceeded with high regioselectivity. NMR analyses of the Y 3 N@ I h -C 80 adducts found that the tris-adducts were all-[6,6]- or [6,6][6,6][5,6]-isomers and that some showed mutual isomerization or remained intact at room temperature. The tetra-adduct obtained as a major product was all-[6,6] and stable. For the structural elucidation of Gd 3 N@ I h -C 80 tris- and tetra-adducts, density functional theory (DFT) calculations were performed to estimate the relative stabilities of tris- and tetra-adducts formed upon Prato functionalization of the most pyramidalized regions of the fullerene structure. The most stable structures corresponded to additions on the most pyramidalized (i.e., strained) bonds. Taking together the experimental vis-NIR spectra, NMR assignments, and the computed relative DFT stabilities of the potential tris- and tetra-adducts, the structures of the isolated adducts were elucidated. Electron resonance (ESR) measurements measurements of pristine, bis-, and tris-adducts of Gd 3 N@C 80 suggested that the rotation of the endohedral metal cluster slowed upon increase of the addition numbers to C 80 cage, which is favored for accommodating the Gd atoms of the relatively large Gd 3 N cluster inner space at the sp 3 addition sites. This is presumably related to the high regioselectivity in the Prato addition reaction driven by the strain release of the Gd 3 N@C 80 fullerene structure.

Evolution has converged on cation-π interactions for recognition of quaternary alkyl ammonium groups such as trimethyllysine (Kme3). While computational modelling indicates that Trp provides the strongest cation-π interaction of the native aromatic amino acids, there is limited corroborative data from measurements within proteins. Herein we investigate a Tyr to Trp mutation in the binding pocket of the HP1 chromodomain, a reader protein that recognizes Kme3. Binding studies demonstrate that the Trp-mediated cation-π interaction is about -5 kcal mol -1 stronger, and the Y24W crystal structure shows that the mutation is not perturbing. Quantum mechanical calculations indicate that greater enthalpic binding is predominantly due to increased cation-π interactions. NMR studies indicate that differences in the unbound state of the Y24W mutation lead to enthalpy-entropy compensation. These results provide direct experimental quantification of Trp versus Tyr in a cation-π interaction and afford insight into the conservation of aromatic cage residues in Kme3 reader domains., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

The structures of two bis-ethylpyrrolidinoadducts of Gd 3 N@ I h - C 80 , obtained by regioselective 1,3-dipolar cycloadditions, were elucidated by single crystal X-ray, visible-near infrared (vis-NIR) spectra, studies on their thermal isomerization, and theoretical calculations. The structure of the minor-bis-adduct reveals a C 2 -symmetric carbon cage with [6,6][6,6]-addition sites and with an endohedral Gd 3 N cluster that is completely flattened. This is the first example of a crystal structure of Gd 3 N@ I h - C 80 derivatives. The structure of the major-bis-adduct was inferred by the vis-NIR spectrum being corresponded to the structure of a previously reported major-bis-adduct of Y 3 N@ I h - C 80 known to have an asymmetric [6,6][6,6]-structure. Based on experimental results showing that the minor-bis-adduct of Gd 3 N@ I h - C 80 isomerized to the major-adduct, a possible second addition site was elucidated with support from density functional theory calculations.

Three Prato monoadduct isomers were synthesized and structurally characterized by 1 H, 13 C NMR spectra and single-crystal X-ray diffraction, and one adduct on the dd-[5,6]-bond was found as the first example of a Prato [5,6]-adduct of C 70 . To investigate the mechanism in the generation of this dd-[5,6]-adduct, computational studies were employed to show that it was thermodynamically obtained by sigmatropic rearrangement from the presumed initial kinetic product de-[6,6]-adduct.

There are few biocatalytic transformations that produce fluorine-containing molecules prevalent in modern pharmaceuticals. To expand the scope of biocatalysis for organofluorine synthesis, we have developed an enzymatic platform for highly enantioselective carbene B-H bond insertion to yield versatile α- trifluoromethylated (α-CF 3 ) organoborons, an important class of organofluorine molecules that contain stereogenic centers bearing both CF 3 and boron groups. In contrast to current "carbene transferase" enzymes that use a limited set of simple diazo compounds as carbene precursors, this system based on Rhodothermus marinus cytochrome c ( Rma cyt c ) can accept a broad range of trifluorodiazo alkanes and deliver versatile chiral α-CF 3 organoborons with total turnovers up to 2870 and enantiomeric ratios up to 98.5:1.5. Computational modeling reveals that this broad diazo scope is enabled by an active-site environment that directs the alkyl substituent on the heme CF 3 -carbene intermediate toward the solvent-exposed face, thereby allowing the protein to accommodate diazo compounds with diverse structural features., Competing Interests: The authors declare no competing financial interest.

A novel naphthalene-based 5·(BArF)8 capsule allows for the size-selective inclusion of C60 from fullerene mixtures. Its size selectivity towards C60 has been rationalized by its dynamic adaptability in solution that has been investigated by molecular dynamics. Additionally, 5·(BArF)8 encapsulates C60-derivatives such as C60-PCBM and N-methylpyrrolidine-C60. The latter can be separated from C60 since 5·(BArF)8 displays distinct affinity for them.

The conformational landscape of Bacillus megaterium epoxide hydrolase (BmEH) and how it is altered by mutations that confer the enzyme the ability to accept bulky epoxide substrates has been investigated. Extensive molecular dynamics (MD) simulations coupled to active site volume calculations have unveiled relevant features of the enzyme conformational dynamics and function. Our long-timescale MD simulations identify key conformational states not previously observed by means of X-ray crystallography and short MD simulations that present the loop containing one of the catalytic residues, Asp239, in a wide-open conformation, which is likely involved in the binding of the epoxide substrate. Introduction of mutations M145S and F128A dramatically alters the conformational landscape of the enzyme. These singly mutated variants can accept bulky epoxide substrates due to the disorder induced by mutation in the α-helix containing the catalytic Tyr144 and some parts of the lid domain. These changes impact the enzyme active site, which is substantially wider and more complementary to the bulky pharmacologically relevant epoxide substrates., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Controlling the regioselectivity of Baeyer-Villiger (BV) reactions remains an ongoing issue in organic chemistry, be it by synthetic catalysts or enzymes of the type Baeyer-Villiger monooxygenases (BVMOs). Herein, we address the challenging problem of switching normal to abnormal BVMO regioselectivity by directed evolution using three linear ketones as substrates, which are not structurally biased toward abnormal reactivity. Upon applying iterative saturation mutagenesis at sites lining the binding pocket of the thermostable BVMO from Thermocrispum municipale DSM 44069 (TmCHMO) and using 4-phenyl-2-butanone as substrate, the regioselectivity was reversed from 99:1 (wild-type enzyme in favor of the normal product undergoing 2-phenylethyl migration) to 2:98 in favor of methyl migration when applying the best mutant. This also stands in stark contrast to the respective reaction using the synthetic reagent m-CPBA, which provides solely the normal product. Reversal of regioselectivity was also achieved in the BV reaction of two other linear ketones. Kinetic parameters and melting temperatures revealed that most of the evolved mutants retained catalytic activity, as well as thermostability. In order to shed light on the origin of switched regioselectivity in reactions of 4-phenyl-2-butanone and phenylacetone, extensive QM/MM and MD simulations were performed. It was found that the mutations introduced by directed evolution induce crucial changes in the conformation of the respective Criegee intermediates and transition states in the binding pocket of the enzyme. In mutants that destabilize the normally preferred migration transition state, a reversal of regioselectivity is observed. This conformational control of regioselectivity overrides electronic control, which normally causes preferential migration of the group that is best able to stabilize positive charge. The results can be expected to aid future protein engineering of BVMOs.

Electrides are ionic substances containing isolated electrons. These confined electrons are topologically characterised by a quasi-atom, that is, a non-nuclear attractor (NNA) of the electron density. The electronic structure of the octahedral 4 A 1g Li 6 + and 5 A 1g Be 6 species shows that these species have a large number of NNAs. These NNAs have highly delocalised electron densities and, as a result, the chemical bonding pattern of these systems is reminiscent of that in solid metals, in which metal cations are surrounded by a "sea" of delocalised valence electrons. We propose the term metal cluster electrides to refer to this new class of compounds. In this study, we establish a computational protocol to identify, characterize, and design metal cluster electrides and we elucidate the intricate bonding patterns of this particular type of species., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Recently, heme proteins have been discovered and engineered by directed evolution to catalyze chemical transformations that are biochemically unprecedented. Many of these nonnatural enzyme-catalyzed reactions are assumed to proceed through a catalytic iron porphyrin carbene (IPC) intermediate, although this intermediate has never been observed in a protein. Using crystallographic, spectroscopic, and computational methods, we have captured and studied a catalytic IPC intermediate in the active site of an enzyme derived from thermostable Rhodothermus marinus ( Rma ) cytochrome c High-resolution crystal structures and computational methods reveal how directed evolution created an active site for carbene transfer in an electron transfer protein and how the laboratory-evolved enzyme achieves perfect carbene transfer stereoselectivity by holding the catalytic IPC in a single orientation. We also discovered that the IPC in Rma cytochrome c has a singlet ground electronic state and that the protein environment uses geometrical constraints and noncovalent interactions to influence different IPC electronic states. This information helps us to understand the impressive reactivity and selectivity of carbene transfer enzymes and offers insights that will guide and inspire future engineering efforts., Competing Interests: The authors declare no conflict of interest.

Microbial aromatic catabolism offers a promising approach to convert lignin, a vast source of renewable carbon, into useful products. Aryl-O-demethylation is an essential biochemical reaction to ultimately catabolize coniferyl and sinapyl lignin-derived aromatic compounds, and is often a key bottleneck for both native and engineered bioconversion pathways. Here, we report the comprehensive characterization of a promiscuous P450 aryl-O-demethylase, consisting of a cytochrome P450 protein from the family CYP255A (GcoA) and a three-domain reductase (GcoB) that together represent a new two-component P450 class. Though originally described as converting guaiacol to catechol, we show that this system efficiently demethylates both guaiacol and an unexpectedly wide variety of lignin-relevant monomers. Structural, biochemical, and computational studies of this novel two-component system elucidate the mechanism of its broad substrate specificity, presenting it as a new tool for a critical step in biological lignin conversion.

Duclauxins are dimeric and heptacyclic fungal polyketides with notable bioactivities. We characterized the cascade of redox transformations in the biosynthetic pathway of duclauxin from Talaromyces stipitatus. The redox reaction sequence is initiated by a cupin family dioxygenase DuxM that performs an oxidative cleavage of the peri-fused tricyclic phenalenone and affords a transient hemiketal-oxaphenalenone intermediate. Additional redox enzymes then morph the oxaphenoalenone into either an anhydride or a dihydrocoumarin-containing monomeric building block that is found in dimeric duxlauxins. Oxidative coupling between the monomers to form the initial C-C bond was shown to be catalyzed by a P450 monooxygenase, although the enzyme responsible for the second C-C bond formation was not found in the pathway. Collectively, the number and variety of redox enzymes used in the duclauxin pathway showcase Nature's strategy to generate structural complexity during natural product biosynthesis.

Controlling the regioselectivity in the exohedral functionalization of fullerenes and endohedral metallofullerenes is essential to produce specific desired fullerene derivatives. In this work, using density functional theory (DFT) calculations, we show that the regioselectivity of the Diels-Alder (DA) cycloaddition of cyclopentadiene to 2S+1C60 changes from the usual [6,6] addition in the singlet ground state to the [5,6] attack in high spin states of C60. Changes in the aromaticity of the five- and six-membered rings when going from singlet to high spin C60 provide a rationale to understand this regioselectivity change. Experimentally, however, we find that the DA cycloaddition of isoindene to triplet C60 yields the usual [6,6] adduct. Further DFT calculations and computational analysis give an explanation to this unanticipated experimental result by showing the presence of an intersystem crossing close to the formed triplet biradical intermediate.

Hapalindole alkaloids are a structurally diverse class of cyanobacterial natural products defined by their varied polycyclic ring systems and diverse biological activities. These complex metabolites are generated from a common biosynthetic intermediate by the Stig cyclases in three mechanistic steps: a rare Cope rearrangement, 6-exo-trig cyclization, and electrophilic aromatic substitution. Here we report the structure of HpiC1, a Stig cyclase that catalyzes the formation of 12-epi-hapalindole U in vitro. The 1.5-Å structure revealed a dimeric assembly with two calcium ions per monomer and with the active sites located at the distal ends of the protein dimer. Mutational analysis and computational methods uncovered key residues for an acid-catalyzed [3,3]-sigmatropic rearrangement, as well as specific determinants that control the position of terminal electrophilic aromatic substitution, leading to a switch from hapalindole to fischerindole alkaloids.

Malbrancheamide is a dichlorinated fungal indole alkaloid isolated from both Malbranchea aurantiaca and Malbranchea graminicola that belongs to a family of natural products containing a characteristic bicyclo[2.2.2]diazaoctane core. The introduction of chlorine atoms on the indole ring of malbrancheamide differentiates it from other members of this family and contributes significantly to its biological activity. In this study, we characterized the two flavin-dependent halogenases involved in the late-stage halogenation of malbrancheamide in two different fungal strains. MalA and MalA' catalyze the iterative dichlorination and monobromination of the free substrate premalbrancheamide as the final steps in the malbrancheamide biosynthetic pathway. Two unnatural bromo-chloro-malbrancheamide analogues were generated through MalA-mediated chemoenzymatic synthesis. Structural analysis and computational studies of MalA' in complex with three substrates revealed that the enzyme represents a new class of zinc-binding flavin-dependent halogenases and provides new insights into a potentially unique reaction mechanism.

Epoxides are highly useful synthons and biosynthons for the construction of complex natural products during total synthesis and biosynthesis, respectively. Among enzyme-catalyzed epoxide transformations, a reaction that is notably missing, in regard to the synthetic toolbox, is cationic rearrangement that takes place under strong acid. This is a challenging transformation for enzyme catalysis, as stabilization of the carbocation intermediate upon epoxide cleavage is required. Here, we discovered two Brønsted acid enzymes that can catalyze two unprecedented epoxide transformations in biology. PenF from the penigequinolone pathway catalyzes a cationic epoxide rearrangement under physiological conditions to generate a quaternary carbon center, while AsqO from the aspoquinolone pathway catalyzes a 3-exo-tet cyclization to forge a cyclopropane-tetrahydrofuran ring system. The discovery of these new epoxide-modifying enzymes further highlights the versatility of epoxides in complexity generation during natural product biosynthesis.

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

Diels-Alder cycloaddition is one of the most powerful tools for the functionalization of single-walled carbon nanotubes (SWCNTs). Density functional theory at the B3-LYP level of theory has been used to investigate the reactivity of different-diameter SWCNTs (4-9,5) in Diels-Alder reactions with 1,3-butadiene; the reactivity was found to decrease with increasing SWCNT diameter. Distortion/interaction analysis along the whole reaction pathway was found to be a better way to explore the reactivity of this type of reaction. The difference in interaction energy along the reaction pathway is larger than that of the corresponding distortion energy. However, the distortion energy plots for these reactions show the same trend. Therefore, the formation of the transition state can be determined from the interaction energy. A lower interaction energy leads to an earlier transition state, which indicates a lower activation energy. The computational results also indicate that the original distortion of the SWCNTs leads to an increase in the reactivity of the SWCNTs., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)