Controlling Monoterpene Isomerization by Guiding Challenging Carbocation Rearrangement Reactions in Engineered Squalene-Hopene Cyclases

The interconversion of monoterpenes is facilitated by a complex network of carbocation rearrangement pathways. Controlling these isomerization pathways is challenging when using common Brønsted and Lewis acid catalysts, which often produce product mixtures that are difficult to separate. In contrast, natural monoterpene cyclases exhibit high control over the carbocation rearrangement reactions but are reliant on phosphorylated substrates. In this study, we present engineered squalene-hopene cyclases from Alicyclobacillus acidocaldarius (AacSHC) that catalyze the challenging isomerization of monoterpenes with unprecedented precision. Starting from a promiscuous isomerization of (+)-β-pinene, we first demonstrate noticeable shifts in the product distribution solely by introducing single point mutations. Furthermore, we showcase the tuneable cation steering by enhancing (+)-borneol selectivity from 1 % to >90 % (>99 % de) aided by iterative saturation mutagenesis. Our combined experimental and computational data suggest that the reorganization of key aromatic residues leads to the restructuring of the water network that facilitates the selective termination of the secondary isobornyl cation. This work expands our mechanistic understanding of carbocation rearrangements and sets the stage for target-oriented skeletal reorganization of broadly abundant terpenes.

Fatty acid synthase (FASN) signalome: A molecular guide for precision oncology

The initial excitement generated more than two decades ago by the discovery of drugs targeting fatty acid synthase (FASN)-catalyzed de novo lipogenesis for cancer therapy was short-lived. However, the advent of the first clinical-grade FASN inhibitor (TVB-2640; denifanstat), which is currently being studied in various phase II trials, and the exciting advances in understanding the FASN signalome are fueling a renewed interest in FASN-targeted strategies for the treatment and prevention of cancer. Here, we provide a detailed overview of how FASN can drive phenotypic plasticity and cell fate decisions, mitochondrial regulation of cell death, immune escape and organ-specific metastatic potential. We then present a variety of FASN-targeted therapeutic approaches that address the major challenges facing FASN therapy. These include limitations of current FASN inhibitors and the lack of precision tools to maximize the therapeutic potential of FASN inhibitors in the clinic. Rethinking the role of FASN as a signal transducer in cancer pathogenesis may provide molecularly driven strategies to optimize FASN as a long-awaited target for cancer therapeutics.

The shortest path method (SPM) webserver for computational enzyme design

SPMweb is the online webserver of the Shortest Path Map (SPM) tool for identifying the key conformationally-relevant positions of a given enzyme structure and dynamics. The server is built on top of the DynaComm.py code and enables the calculation and visualization of the SPM pathways. SPMweb is easy-to-use as it only requires three input files: the three-dimensional structure of the protein of interest, and the two matrices (distance and correlation) previously computed from a Molecular Dynamics simulation. We provide in this publication information on how to generate the files for SPM construction even for non-expert users and discuss the most relevant parameters that can be modified. The tool is extremely fast (it takes less than one minute per job), thus allowing the rapid identification of distal positions connected to the active site pocket of the enzyme. SPM applications expand from computational enzyme design, especially if combined with other tools to identify the preferred substitution at the identified position, but also to rationalizing allosteric regulation, and even cryptic pocket identification for drug discovery. The simple user interface and setup make the SPM tool accessible to the whole scientific community. SPMweb is freely available for academia at http://spmosuna.com/.

Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Retracing the Rapid Evolution of an Herbicide-Degrading Enzyme by Protein Engineering

The mechanisms underlying the rapid evolution of novel enzymatic activities from promiscuous side activities are poorly understood. Recently emerged enzymes catalyzing the catabolic degradation of xenobiotic substances that have been spread out into the environment during the last decades provide an exquisite opportunity to study these mechanisms. A prominent example is the herbicide atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine), which is degraded through a number of enzymatic reactions constituting the Atz pathway. Here, we analyzed the evolution of the hydroxyatrazine ethylaminohydrolase AtzB, a Zn(II)-dependent metalloenzyme that adopts the amidohydrolase fold and catalyzes the second step of the Atz pathway. We searched for promiscuous side activities of AtzB, which might point to the identity of its progenitor. These investigations revealed that AtzB has low promiscuous guanine deaminase activity. Furthermore, we found that the two closest AtzB homologues, which have not been functionally annotated up to now, are guanine deaminases with modest promiscuous hydroxyatrazine hydrolase activity. Based on sequence comparisons with the closest AtzB homologues, the guanine deaminase activity of AtzB could be increased by three orders of magnitude through the introduction of only four active site mutations. Interestingly, introducing the inverse four mutations into the AtzB homologues significantly enhanced their hydroxyatrazine hydrolase activity, and in one case is even equivalent to that of wild-type AtzB. Molecular dynamics simulations elucidated the structural and molecular basis for the mutation-induced activity changes. The example of AtzB highlights how novel enzymes with high catalytic proficiency can evolve from low promiscuous side activities by only few mutational events within a short period of time.

Designing Efficient Enzymes: Eight Predicted Mutations Convert a Hydroxynitrile Lyase into an Efficient Esterase

Hydroxynitrile lyase from rubber tree (HbHNL) shares 45% identical amino acid residues with the homologous esterase from tobacco, SABP2, but the two enzymes catalyze different reactions. The x-ray structures reveal a serine-histidine-aspartate catalytic triad in both enzymes along with several differing amino acid residues within the active site. Previous exchange of three amino acid residues in the active site of HbHNL with the corresponding amino acid residue in SABP2 (T11G-E79H-K236M) created variant HNL3, which showed low esterase activity toward p-nitrophenyl acetate. Further structure comparison reveals additional differences surrounding the active site. HbHNL contains an improperly positioned oxyanion hole residue and differing solvation of the catalytic aspartate. We hypothesized that correcting these structural differences would impart good esterase activity on the corresponding HbHNL variant. To predict the amino acid substitutions needed to correct the structure, we calculated shortest path maps for both HbHNL and SABP2, which reveal correlated movements of amino acids in the two enzymes. Replacing four amino acid residues (C81L-N104T-V106F-G176S) whose movements are connected to the movements of the catalytic residues yielded variant HNL7TV (stabilizing substitution H103V was also added), which showed an esterase catalytic efficiency comparable to that of SABP2. The x-ray structure of an intermediate variant, HNL6V, showed an altered solvation of the catalytic aspartate and a partially corrected oxyanion hole. This dramatic increase in catalytic efficiency demonstrates the ability of shortest path maps to predict which residues outside the active site contribute to catalytic activity.

C-H Bonds as Functional Groups: Simultaneous Generation of Multiple Stereocenters by Enantioselective Hydroxylation at Unactivated Tertiary C-H Bonds

Enantioselective C-H oxidation is a standing chemical challenge foreseen as a powerful tool to transform readily available organic molecules into precious oxygenated building blocks. Here, we describe a catalytic enantioselective hydroxylation of tertiary C-H bonds in cyclohexane scaffolds with H2O2, an evolved manganese catalyst that provides structural complementary to the substrate similarly to the lock-and-key recognition operating in enzymatic active sites. Theoretical calculations unveil that enantioselectivity is governed by the precise fitting of the substrate scaffold into the catalytic site, through a network of complementary weak non-covalent interactions. Stereoretentive C(sp3)-H hydroxylation results in a single-step generation of multiple stereogenic centers (up to 4) that can be orthogonally manipulated by conventional methods providing rapid access, from a single precursor to a variety of chiral scaffolds.

AlphaFold2 and Deep Learning for Elucidating Enzyme Conformational Flexibility and Its Application for Design

The recent success of AlphaFold2 (AF2) and other deep learning (DL) tools in accurately predicting the folded three-dimensional (3D) structure of proteins and enzymes has revolutionized the structural biology and protein design fields. The 3D structure indeed reveals key information on the arrangement of the catalytic machinery of enzymes and which structural elements gate the active site pocket. However, comprehending enzymatic activity requires a detailed knowledge of the chemical steps involved along the catalytic cycle and the exploration of the multiple thermally accessible conformations that enzymes adopt when in solution. In this Perspective, some of the recent studies showing the potential of AF2 in elucidating the conformational landscape of enzymes are provided. Selected examples of the key developments of AF2-based and DL methods for protein design are discussed, as well as a few enzyme design cases. These studies show the potential of AF2 and DL for allowing the routine computational design of efficient enzymes.

Harnessing the Structure and Dynamics of the Squalene-Hopene Cyclase for (-)-Ambroxide Production

Terpene cyclases offer enormous synthetic potential, given their unique ability to forge complex hydrocarbon scaffolds from achiral precursors within a single cationic rearrangement cascade. Harnessing their synthetic power, however, proved to be challenging due to their generally low catalytic performance. In this study, we unveil the catalytic potential of the squalene-hopene cyclase (SHC) by harnessing its structure and dynamics. First, we synergistically tailored the enzyme's active site and entrance tunnel to generate a 397-fold improved (-)-ambroxide synthase. Our computational investigations explain how the introduced mutations work in concert to improve substrate acquisition, flow, and chaperoning. Kinetics, however, disclosed a terpene-induced inactivation of the membrane-bound SHC as the major turnover limitation in vivo. Merging these insights with the improved and stereoselective catalysis of the enzyme, we applied a feeding strategy to exceed 105 total turnovers. We believe that our results may bridge the gap for broader application of SHCs in synthetic chemistry.

Brave new surfactant world revisited by thermoalkalophilic lipases: computational insights into the role of SDS as a substrate analog

Non-ionic surfactants were shown to stabilize the active conformation of thermoalkalophilic lipases by mimicking the lipid substrate while the catalytic interactions formed by anionic surfactants have not been well characterized. In this study, we combined μs-scale molecular dynamics (MD) simulations and lipase activity assays to analyze the effect of ionic surfactant, sodium dodecyl sulfate (SDS), on the structure and activity of thermoalkalophilic lipases. Both the open and closed lipase conformations that differ in geometry were recruited to the MD analysis to provide a broader understanding of the molecular effect of SDS on the lipase structure. Simulations at 298 K showed the potential of SDS for maintaining the active lipase through binding to the sn-1 acyl-chain binding pocket in the open conformation or transforming the closed conformation to an open-like state. Consistent with MD findings, experimental analysis showed increased lipase activity upon SDS incubation at ambient temperature. Notably, the lipase cores stayed intact throughout 2 μs regardless of an increase in the simulation temperature or SDS concentration. However, the surface structures were unfolded in the presence of SDS and at elevated temperature for both conformations. Simulations of the dimeric lipase were also carried out and showed reduced flexibility of the surface structures which were unfolded in the monomer, indicating the insulating role of dimer interactions against SDS. Taken together, this study provides insights into the possible substrate mimicry by the ionic surfactant SDS for the thermoalkalophilic lipases without temperature elevation, underscoring SDS's potential for interfacial activation at ambient temperatures.

Synthesis, In Vitro Profiling, and In Vivo Evaluation of Benzohomoadamantane-Based Ureas for Visceral Pain: A New Indication for Soluble Epoxide Hydrolase Inhibitors

The soluble epoxide hydrolase (sEH) has been suggested as a pharmacological target for the treatment of several diseases, including pain-related disorders. Herein, we report further medicinal chemistry around new benzohomoadamantane-based sEH inhibitors (sEHI) in order to improve the drug metabolism and pharmacokinetics properties of a previous hit. After an extensive in vitro screening cascade, molecular modeling, and in vivo pharmacokinetics studies, two candidates were evaluated in vivo in a murine model of capsaicin-induced allodynia. The two compounds showed an anti-allodynic effect in a dose-dependent manner. Moreover, the most potent compound presented robust analgesic efficacy in the cyclophosphamide-induced murine model of cystitis, a well-established model of visceral pain. Overall, these results suggest painful bladder syndrome as a new possible indication for sEHI, opening a new range of applications for them in the visceral pain field.

Estimating conformational heterogeneity of tryptophan synthase with a template‐based Alphafold2 approach

The three‐dimensional structure of the enzymes provides very relevant information on the arrangement of the catalytic machinery and structural elements gating the active site pocket. The recent success of the neural network Alphafold2 in predicting the folded structure of proteins from the primary sequence with high levels of accuracy has revolutionized the protein design field. However, the application of Alphafold2 for understanding and engineering function directly from the obtained single static picture is not straightforward. Indeed, understanding enzymatic function requires the exploration of the ensemble of thermally accessible conformations that enzymes adopt in solution. In the present study, we evaluate the potential of Alphafold2 in assessing the effect of the mutations on the conformational landscape of the beta subunit of tryptophan synthase (TrpB). Specifically, we develop a template‐based Alphafold2 approach for estimating the conformational heterogeneity of several TrpB enzymes, which is needed for enhanced stand‐alone activity. Our results show the potential of Alphafold2, especially if combined with molecular dynamics simulations, for elucidating the changes induced by mutation in the conformational landscapes at a rather reduced computational cost, thus revealing its plausible application in computational enzyme design.

Time Evolution of the Millisecond Allosteric Activation of Imidazole Glycerol Phosphate Synthase

Deciphering the molecular mechanisms of enzymatic allosteric regulation requires the structural characterization of functional states and also their time evolution toward the formation of the allosterically activated ternary complex. The transient nature and usually slow millisecond time scale interconversion between these functional states hamper their experimental and computational characterization. Here, we combine extensive molecular dynamics simulations, enhanced sampling techniques, and dynamical networks to describe the allosteric activation of imidazole glycerol phosphate synthase (IGPS) from the substrate-free form to the active ternary complex. IGPS is a heterodimeric bienzyme complex whose HisH subunit is responsible for hydrolyzing glutamine and delivering ammonia for the cyclase activity in HisF. Despite significant advances in understanding the underlying allosteric mechanism, essential molecular details of the long-range millisecond allosteric activation of IGPS remain hidden. Without using a priori information of the active state, our simulations uncover how IGPS, with the allosteric effector bound in HisF, spontaneously captures glutamine in a catalytically inactive HisH conformation, subsequently attains a closed HisF:HisH interface, and finally forms the oxyanion hole in HisH for efficient glutamine hydrolysis. We show that the combined effector and substrate binding dramatically decreases the conformational barrier associated with oxyanion hole formation, in line with the experimentally observed 4500-fold activity increase in glutamine hydrolysis. The allosteric activation is controlled by correlated time-evolving dynamic networks connecting the effector and substrate binding sites. This computational strategy tailored to describe millisecond events can be used to rationalize the effect of mutations on the allosteric regulation and guide IGPS engineering efforts.

Discovery and In Vivo Proof of Concept of a Highly Potent Dual Inhibitor of Soluble Epoxide Hydrolase and Acetylcholinesterase for the Treatment of Alzheimer's Disease

With innumerable clinical failures of target-specific drug candidates for multifactorial diseases, such as Alzheimer's disease (AD), which remains inefficiently treated, the advent of multitarget drug discovery has brought a new breath of hope. Here, we disclose a class of 6-chlorotacrine (huprine)-TPPU hybrids as dual inhibitors of the enzymes soluble epoxide hydrolase (sEH) and acetylcholinesterase (AChE), a multitarget profile to provide cumulative effects against neuroinflammation and memory impairment. Computational studies confirmed the gorge-wide occupancy of both enzymes, from the main site to a secondary site, including a so far non-described AChE cryptic pocket. The lead compound displayed in vitro dual nanomolar potencies, adequate brain permeability, aqueous solubility, human microsomal stability, lack of neurotoxicity, and it rescued memory, synaptic plasticity, and neuroinflammation in an AD mouse model, after low dose chronic oral administration.

In Silico Identification and Experimental Validation of Distal Activity-Enhancing Mutations in Tryptophan Synthase

Allostery is a central mechanism for the regulation of multi-enzyme complexes. The mechanistic basis that drives allosteric regulation is poorly understood but harbors key information for enzyme engineering. In the present study, we focus on the tryptophan synthase complex that is composed of TrpA and TrpB subunits, which allosterically activate each other. Specifically, we develop a rational approach for identifying key amino acid residues of TrpB distal from the active site. Those residues are predicted to be crucial for shifting the inefficient conformational ensemble of the isolated TrpB to a productive ensemble through intra-subunit allosteric effects. The experimental validation of the conformationally driven TrpB design demonstrates its superior stand-alone activity in the absence of TrpA, comparable to those enhancements obtained after multiple rounds of experimental laboratory evolution. Our work evidences that the current challenge of distal active site prediction for enhanced function in computational enzyme design has become within reach.

SBMOpenMM: A Builder of Structure-Based Models for OpenMM

Molecular dynamics (MD) simulations have become a standard tool to correlate the structure and function of biomolecules and significant advances have been made in the study of proteins and their complexes. A major drawback of conventional MD simulations is the difficulty and cost of obtaining converged results, especially when exploring potential energy surfaces containing considerable energy barriers. This limits the wide use of MD calculations to determine the thermodynamic properties of biomolecular processes. Indeed, this is true when considering the conformational entropy of such processes, which is ultimately critical in assessing the simulations' convergence. Alternatively, a wide range of structure-based models (SBMs) has been used in the literature to unravel the basic mechanisms of biomolecular dynamics. These models introduce simplifications that focus on the relevant aspects of the physical process under study. Because of this, SBMs incorporate the need to modify the force field definition and parameters to target specific biophysical simulations. Here we introduce SBMOpenMM, a Python library to build force fields for SBMs, that uses the OpenMM framework to create and run SBM simulations. The code is flexible and user-friendly and profits from the high customizability and performance provided by the OpenMM platform.

From the Design to the In Vivo Evaluation of Benzohomoadamantane-Derived Soluble Epoxide Hydrolase Inhibitors for the Treatment of Acute Pancreatitis

The pharmacological inhibition of soluble epoxide hydrolase (sEH) is efficient for the treatment of inflammatory and pain-related diseases. Numerous potent sEH inhibitors (sEHIs) present adamantyl or phenyl moieties, such as the clinical candidates AR9281 or EC5026. Herein, in a new series of sEHIs, these hydrophobic moieties have been merged in a benzohomoadamantane scaffold. Most of the new sEHIs have excellent inhibitory activities against sEH. Molecular dynamics simulations suggested that the addition of an aromatic ring into the adamantane scaffold produced conformational rearrangements in the enzyme to stabilize the aromatic ring of the benzohomoadamantane core. A screening cascade permitted us to select a candidate for an in vivo efficacy study in a murine model of cerulein-induced acute pancreatitis. The administration of 22 improved the health status of the animals and reduced pancreatic damage, demonstrating that the benzohomoadamantane unit is a promising scaffold for the design of novel sEHIs.

Insights into the molecular determinants of thermal stability in halohydrin dehalogenase HheD2

Halohydrin dehalogenases (HHDHs) are promising enzymes for application in biocatalysis due to their promiscuous epoxide ring-opening activity with various anionic nucleophiles. So far, seven different HHDH subtypes A to G have been reported with subtype D containing the by far largest number of enzymes. Moreover, several characterized members of subtype D have been reported to display outstanding characteristics such as high catalytic activity, broad substrate spectra or remarkable thermal stability. Yet, no structure of a D-type HHDH has been reported to date that could be used to investigate and understand those features on a molecular level. We therefore solved the crystal structure of HheD2 from gamma proteobacterium HTCC2207 at 1.6 Å resolution and used it as a starting point for targeted mutagenesis in combination with molecular dynamics (MD) simulation, in order to study the low thermal stability of HheD2 in comparison with other members of subtype D. This revealed a hydrogen bond between conserved residues Q160 and D198 to be connected with a high catalytic activity of this enzyme. Moreover, a flexible surface region containing two α-helices was identified to impact thermal stability of HheD2. Exchange of this surface region by residues of HheD3 yielded a variant with 10 °C higher melting temperature and reaction temperature optimum. Overall, our results provide important insights into the structure-function relationship of HheD2 and presumably for other D-type HHDHs. DATABASES: Structural data are available in PDB database under the accession number 7B73.

2-Oxaadamant-1-yl Ureas as Soluble Epoxide Hydrolase Inhibitors: In Vivo Evaluation in a Murine Model of Acute Pancreatitis

In vivo pharmacological inhibition of soluble epoxide hydrolase (sEH) reduces inflammatory diseases, including acute pancreatitis (AP). Adamantyl ureas are very potent sEH inhibitors, but the lipophilicity and metabolism of the adamantane group compromise their overall usefulness. Herein, we report that the replacement of a methylene unit of the adamantane group by an oxygen atom increases the solubility, permeability, and stability of three series of urea-based sEH inhibitors. Most of these oxa-analogues are nanomolar inhibitors of both the human and murine sEH. Molecular dynamics simulations rationalize the molecular basis for their activity and suggest that the presence of the oxygen atom on the adamantane scaffold results in active site rearrangements to establish a weak hydrogen bond. The 2-oxaadamantane 22, which has a good solubility, microsomal stability, and selectivity for sEH, was selected for further in vitro and in vivo studies in models of cerulein-induced AP. Both in prophylactic and treatment studies, 22 diminished the overexpression of inflammatory and endoplasmic reticulum stress markers induced by cerulein and reduced the pancreatic damage.

Regio- and Stereoselective Steroid Hydroxylation at C7 by Cytochrome P450 Monooxygenase Mutants

Steroidal C7β alcohols and their respective esters have shown significant promise as neuroprotective and anti-inflammatory agents to treat chronic neuronal damage like stroke, brain trauma, and cerebral ischemia. Since C7 is spatially far away from any functional groups that could direct C-H activation, these transformations are not readily accessible using modern synthetic organic techniques. Reported here are P450-BM3 mutants that catalyze the oxidative hydroxylation of six different steroids with pronounced C7 regioselectivities and β stereoselectivities, as well as high activities. These challenging transformations were achieved by a focused mutagenesis strategy and application of a novel technology for protein library construction based on DNA assembly and USER (Uracil-Specific Excision Reagent) cloning. Upscaling reactions enabled the purification of the respective steroidal alcohols in moderate to excellent yields. The high-resolution X-ray structure and molecular dynamics simulations of the best mutant unveil the origin of regio- and stereoselectivity.

Regioselective Synthesis and Characterization of Tris- and Tetra-Prato Adducts of M3N@C80 (M = Y, Gd)

The tris- and tetra-adducts of M₃N@I_h-C₈₀ metallofullerenes were synthesized and characterized for the first time. The 1,3-dipolar cycloaddition (Prato reaction) of Y₃N@I_h-C₈₀ and Gd₃N@I_h-C₈₀ 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₃N@I_h-C₈₀ and Gd₃N@I_h-C₈₀. 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₃N@I_h-C₈₀ 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₃N@I_h-C₈₀ 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₃N@C₈₀ suggested that the rotation of the endohedral metal cluster slowed upon increase of the addition numbers to C₈₀ cage, which is favored for accommodating the Gd atoms of the relatively large Gd₃N cluster inner space at the sp³ addition sites. This is presumably related to the high regioselectivity in the Prato addition reaction driven by the strain release of the Gd₃N@C₈₀ fullerene structure.

In Vivo Selection for Formate Dehydrogenases with High Efficiency and Specificity toward NADP. ACS Catal 2020;10:7512-7525. [PMID: 32733773 PMCID: PMC7384739 DOI: 10.1021/acscatal.0c01487]

The efficient regeneration of cofactors is vital for the establishment of biocatalytic processes. Formate is an ideal electron donor for cofactor regeneration due to its general availability, low reduction potential, and benign byproduct (CO₂). However, formate dehydrogenases (FDHs) are usually specific to NAD⁺, such that NADPH regeneration with formate is challenging. Previous studies reported naturally occurring FDHs or engineered FDHs that accept NADP⁺, but these enzymes show low kinetic efficiencies and specificities. Here, we harness the power of natural selection to engineer FDH variants to simultaneously optimize three properties: kinetic efficiency with NADP⁺, specificity toward NADP⁺, and affinity toward formate. By simultaneously mutating multiple residues of FDH from Pseudomonas sp. 101, which exhibits practically no activity toward NADP⁺, we generate a library of >10⁶ variants. We introduce this library into an E. coli strain that cannot produce NADPH. By selecting for growth with formate as the sole NADPH source, we isolate several enzyme variants that support efficient NADPH regeneration. We find that the kinetically superior enzyme variant, harboring five mutations, has 5-fold higher efficiency and 14-fold higher specificity in comparison to the best enzyme previously engineered, while retaining high affinity toward formate. By using molecular dynamics simulations, we reveal the contribution of each mutation to the superior kinetics of this variant. We further determine how nonadditive epistatic effects improve multiple parameters simultaneously. Our work demonstrates the capacity of in vivo selection to identify highly proficient enzyme variants carrying multiple mutations which would be almost impossible to find using conventional screening methods.

Computational NMR Spectra of o-Benzyne and Stable Guests and Their Hemicarceplexes

The incarceration of o-benzyne and 27 other guest molecules within hemicarcerand 1, as reported experimentally by Warmuth, and Cram and co-workers, has been studied by density functional theory (DFT). The ¹ H NMR chemical shifts, rotational mobility, and conformational preference of the guests within the supramolecular cage were determined, which showed intriguing correlations of the chemical shifts with structural parameters of the host-guest system. Furthermore, based on the computed chemical shifts reassignments of some NMR signals are proposed. This affects, in particular, the putative characterization of the volatile benzyne molecule inside a hemicarcerand, for which our CCSD(T) and KT2 results indicate that the experimentally observed signals are most likely not resulting from an isolated o-benzyne within the supramolecular host. Instead, it is shown that the guest reacted with an aromatic ring of the host, and this adduct is responsible for the experimentally observed signals.

Protein-directed crystalline 2D fullerene assemblies

Water soluble 2D crystalline monolayers of fullerenes grow on planar assemblies of engineered consensus tetratricopeptide repeat proteins. Designed fullerene-coordinating tyrosine clamps on the protein introduce specific fullerene binding sites, which facilitate fullerene nucleation. Through reciprocal interactions between the components, the hybrid material assembles into two-dimensional 2 nm thick structures with crystalline order, that conduct photo-generated charges. Thus, the protein-fullerene hybrid material is a demonstration of the developments toward functional materials with protein-based precision control of functional elements.

Site-Selectivity of Prato Additions to C70: Experimental and Theoretical Studies of a New Thermodynamic Product at the dd-[5,6]-Junction

Three Prato monoadduct isomers were synthesized and structurally characterized by ¹H, ¹³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₇₀. 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.

Structures of Gd3N@C80 Prato Bis-Adducts: Crystal Structure, Thermal Isomerization, and Computational Study

The structures of two bis-ethylpyrrolidinoadducts of Gd₃N@I_h-C₈₀, 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₂-symmetric carbon cage with [6,6][6,6]-addition sites and with an endohedral Gd₃N cluster that is completely flattened. This is the first example of a crystal structure of Gd₃N@I_h-C₈₀ 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₃N@I_h-C₈₀ known to have an asymmetric [6,6][6,6]-structure. Based on experimental results showing that the minor-bis-adduct of Gd₃N@I_h-C₈₀ isomerized to the major-adduct, a possible second addition site was elucidated with support from density functional theory calculations.

Intrinsic enzymatic properties modulate the self-propulsion of micromotors

Bio-catalytic micro- and nanomotors self-propel by the enzymatic conversion of substrates into products. Despite the advances in the field, the fundamental aspects underlying enzyme-powered self-propulsion have rarely been studied. In this work, we select four enzymes (urease, acetylcholinesterase, glucose oxidase, and aldolase) to be attached on silica microcapsules and study how their turnover number and conformational dynamics affect the self-propulsion, combining both an experimental and molecular dynamics simulations approach. Urease and acetylcholinesterase, the enzymes with higher catalytic rates, are the only enzymes capable of producing active motion. Molecular dynamics simulations reveal that urease and acetylcholinesterase display the highest degree of flexibility near the active site, which could play a role on the catalytic process. We experimentally assess this hypothesis for urease micromotors through competitive inhibition (acetohydroxamic acid) and increasing enzyme rigidity (β-mercaptoethanol). We conclude that the conformational changes are a precondition of urease catalysis, which is essential to generate self-propulsion.

Self-propulsion of biocatalytic micro- and nanomotors is facilitated by enzymes converting substrates into products. Here, the authors show that intrinsic enzymatic properties such as conformational changes are crucial for the self-propulsion of silica microcapsules modified with urease.

Role of conformational dynamics in the evolution of novel enzyme function

The free energy landscape concept that describes enzymes as an ensemble of differently populated conformational sub-states in dynamic equilibrium is key for evaluating enzyme activity, enantioselectivity, and specificity. Mutations introduced in the enzyme sequence can alter the populations of the pre-existing conformational states, thus strongly modifying the enzyme ability to accommodate alternative substrates, revert its enantiopreferences, and even increase the activity for some residual promiscuous reactions. In this feature article, we present an overview of the current experimental and computational strategies to explore the conformational free energy landscape of enzymes. We provide a series of recent publications that highlight the key role of conformational dynamics for the enzyme evolution towards new functions and substrates, and provide some perspectives on how conformational dynamism should be considered in future computational enzyme design protocols.

Exploring the reversal of enantioselectivity on a zinc-dependent alcohol dehydrogenase

Alcohol Dehydrogenase (ADH) enzymes catalyse the reversible reduction of prochiral ketones to the corresponding alcohols. These enzymes present two differently shaped active site pockets, which dictate their substrate scope and selectivity. In this study, we computationally evaluate the effect of two commonly reported active site mutations (I86A, and W110T) on a secondary alcohol dehydrogenase from Thermoanaerobacter brockii (TbSADH) through Molecular Dynamics simulations. Our results indicate that the introduced mutations induce dramatic changes in the shape of the active site, but most importantly they impact the substrate-enzyme interactions. We demonstrate that the combination of Molecular Dynamics simulations with the tools POVME and NCIplot corresponds to a powerful strategy for rationalising and engineering the stereoselectivity of ADH variants.

Inducing high activity of a thermophilic enzyme at ambient temperatures by directed evolution

The long-standing problem of achieving high activity of a thermophilic enzyme at low temperatures and short reaction times with little tradeoff in thermostability has been solved by directed evolution, an alcohol dehydrogenase found in hot springs serving as the catalyst in enantioselective ketone reductions.

Effect of incarcerated HF on the exohedral chemical reactivity of HF@C60

The first chemical modification on the brand new endohedral HF@C₆₀ 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 C₆₀ (2a) and endohedral H₂O@C₆₀ (3). The comparative study shows a kinetic order for the isomerization process of H₂O@C₆₀ > HF@C₆₀ > C₆₀, thus confirming the effect of the incarcerated species on the zwitterionic intermediate stability.

Exploring the origins of selectivity in soluble epoxide hydrolase from Bacillus megaterium

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. In this work, the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative are computationally investigated using density functional theory (DFT).

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. These enzymes present different enantio- and regioselectivity depending upon either the substrate structure or the substitution pattern of the epoxide ring. In this study, we computationally investigate the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative using density functional theory (DFT) and an active site cluster model consisting of 195 and 197 atoms, respectively. Full reaction mechanisms for epoxide ring opening were evaluated considering the attack at both oxirane carbons and considering two possible orientations of the substrate at the BmEH active site. Our results indicate that for both SO and p-NSO substrates the BmEH enantio- and regioselectivity is opposite to the inherent (R)-BmEH selectivity, the attack at the benzylic position (C1) of the (S)-enantiomer being the most favoured chemical outcome.

Role of Conformational Dynamics in the Evolution of Retro-Aldolase Activity

Enzymes exist as ensembles of conformations that are important for function. Tuning these populations of conformational states through mutation enables evolution toward additional activities. Here we computationally evaluate the population shifts induced by distal and active site mutations in a family of computationally designed and experimentally optimized retro-aldolases. The conformational landscape of these enzymes was significantly altered during evolution, as pre-existing catalytically active conformational substates became major states in the most evolved variants. We further demonstrate that key residues responsible for these substate conversions can be predicted computationally. Significantly, the identified residues coincide with those positions mutated in the laboratory evolution experiments. This study establishes that distal mutations that affect enzyme catalytic activity can be predicted computationally and thus provides the enzyme (re)design field with a rational strategy to determine promising sites for enhancing activity through mutation.

Rationalizing the relative abundances of trimetallic nitride template-based endohedral metallofullerenes from aromaticity measures

The synthesis of endohedral metallofullerenes (EMFs) from a carbon soot sample of an arc discharge leads to a variety of EMFs that are obtained in different relative abundances. In the present work, we show that these abundances can be predicted from aromaticity calculations. In particular, we use the normalized Additive Local Aromaticity (ALAN) index. Our results show that the most abundant Sc3N-based and Y3N-based EMFs in fullerene soot are the most aromatic. This study reinforces the idea that aromaticity plays a key role in determining the stability of EMFs.

Computational tools for the evaluation of laboratory-engineered biocatalysts

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.

Reactivity of Single-Walled Carbon Nanotubes in the Diels-Alder Cycloaddition Reaction: Distortion-Interaction Analysis along the Reaction Pathway

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.

Reaction Mechanism and Regioselectivity of the Bingel-Hirsch Addition of Dimethyl Bromomalonate to La@C2v -C82

We quantum chemically explore the thermodynamics and kinetics of all 65 possible mechanistic pathways of the Bingel-Hirsch addition of dimethyl bromomalonate to the endohedral metallofullerene La@C2v -C82 that result from the combination of 24 nonequivalent carbon atoms and 35 different bonds present in La@C2v -C82 by using dispersion-corrected DFT calculations. Experimentally, this reaction leads to four singly bonded derivatives and one fulleroid adduct. Of these five products, only the singly bonded derivative on C23 could be experimentally identified unambiguously. Our calculations show that La@C2v -C82 is not particularly regioselective under Bingel-Hirsch conditions. From the obtained results, however, it is possible to make a tentative assignment of the products observed experimentally. We propose that the observed fulleroid adduct results from the attack at bond 19 and that the singly bonded derivatives correspond to the C2, C19, C21, and C23 initial attacks. However, other possibilities cannot be ruled out completely.