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/.

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.

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.

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, has proved to be challenging owing to their generally low catalytic performance. In this study, we unveiled the catalytic potential of the squalene-hopene cyclase (SHC) by harnessing its structure and dynamics. First, we synergistically tailored the active site and entrance tunnel of the enzyme 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, showed terpene-induced inactivation of the membrane-bound SHC to be the major turnover limitation in vivo. Merging this insight 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.

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.

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.

Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics.

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.

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.

Enzyme Conformation Influences the Performance of Lipase-powered Nanomotors.

Enzyme-powered micro/nanomotors have myriads of potential applications in various areas. To efficiently reach those applications, it is necessary and critical to understand the fundamental aspects affecting the motion dynamics. Herein, we explored the impact of enzyme orientation on the performance of lipase-powered nanomotors by tuning the lipase immobilization strategies. The influence of the lipase orientation and lid conformation on substrate binding and catalysis was analyzed using molecular dynamics simulations. Besides, the motion performance indicates that the hydrophobic binding (via OTES) represents the best orienting strategy, providing 48.4 % and 95.4 % increase in diffusion coefficient compared to hydrophilic binding (via APTES) and Brownian motion (no fuel), respectively (with C[triacetin] of 100 mm). This work provides vital evidence for the importance of immobilization strategy and corresponding enzyme orientation for the catalytic activity and in turn, the motion performance of nanomotors, and is thus helpful to future applications.