Biocatalytic stereocontrolled head-to-tail cyclizations of unbiased terpenes as a tool in chemoenzymatic synthesis

Terpene synthesis stands at the forefront of modern synthetic chemistry and represents the state-of-the-art in the chemist's toolbox. Notwithstanding, these endeavors are inherently tied to the current availability of natural cyclic building blocks. Addressing this limitation, the stereocontrolled cyclization of abundant unbiased linear terpenes emerges as a valuable tool, which is still difficult to achieve with chemical catalysts. In this study, we showcase the remarkable capabilities of squalene-hopene cyclases (SHCs) in the chemoenzymatic synthesis of head-to-tail-fused terpenes. By combining engineered SHCs and a practical reaction setup, we generate ten chiral scaffolds with >99% ee and de, at up to decagram scale. Our mechanistic insights suggest how cyclodextrin encapsulation of terpenes may influence the performance of the membrane-bound enzyme. Moreover, we transform the chiral templates to valuable (mero)-terpenes using interdisciplinary synthetic methods, including a catalytic ring-contraction of enol-ethers facilitated by cooperative iodine/lipase catalysis.

Ancestral terpene cyclases: From fundamental science to applications in biosynthesis

Terpenes constitute one of the largest family of natural products with potent applications as renewable platform chemicals and medicines. The low activity, selectivity and stability displayed by terpene biosynthetic machineries can constitute an obstacle towards achieving expedient biosynthesis of terpenoids in processes that adhere to the 12 principles of green chemistry. Accordingly, engineering of terpene synthase enzymes is a prerequisite for industrial biotechnology applications, but obstructed by their complex catalysis that depend on reactive carbocationic intermediates that are prone to undergo bifurcation mechanisms. Rational redesign of terpene synthases can be tedious and requires high-resolution structural information, which is not always available. Furthermore, it has proven difficult to link sequence space of terpene synthase enzymes to specific product profiles. Herein, the author shows how ancestral sequence reconstruction (ASR) can favorably be used as a protein engineering tool in the redesign of terpene synthases without the need of a structure, and without excessive screening. A detailed workflow of ASR is presented along with associated limitations, with a focus on applying this methodology on terpene synthases. From selected examples of both class I and II enzymes, the author advocates that ancestral terpene cyclases constitute valuable assets to shed light on terpene-synthase catalysis and in enabling accelerated biosynthesis.

Tunnel engineering of gas-converting enzymes for inhibitor retardation and substrate acceleration

Carbon monoxide dehydrogenase (CODH), formate dehydrogenase (FDH), hydrogenase (H2ase), and nitrogenase (N2ase) are crucial enzymatic catalysts that facilitate the conversion of industrially significant gases such as CO, CO2, H2, and N2. The tunnels in the gas-converting enzymes serve as conduits for these low molecular weight gases to access deeply buried catalytic sites. The identification of the substrate tunnels is imperative for comprehending the substrate selectivity mechanism underlying these gas-converting enzymes. This knowledge also holds substantial value for industrial applications, particularly in addressing the challenges associated with separation and utilization of byproduct gases. In this comprehensive review, we delve into the emerging field of tunnel engineering, presenting a range of approaches and analyses. Additionally, we propose methodologies for the systematic design of enzymes, with the ultimate goal of advancing protein engineering strategies.

Evolution-Guided Biosynthesis of Terpenoid Inhibitors

Terpenoids, the largest and most structurally diverse group of natural products, include a striking variety of biologically active compounds, from flavors to medicines. Despite their well-documented biochemical versatility, the evolutionary processes that generate new functional terpenoids are poorly understood and difficult to recapitulate in engineered systems. This study uses a synthetic biochemical objective─a transcriptional system that links the inhibition of protein tyrosine phosphatase 1B (PTP1B), a human drug target, to the expression of a gene for antibiotic resistance in Escherichia coli (E. coli)─to evolve a terpene synthase to produce enzyme inhibitors. Site saturation mutagenesis of poorly conserved residues on γ-humulene synthase (GHS), a promicuous enzyme, yielded mutants that improved fitness (i.e., the antibiotic resistance of E. coli) by reducing GHS toxicity and/or by increasing inhibitor production. Intriguingly, a combination of two mutations enhanced the titer of a minority product─a terpene alcohol that inhibits PTP1B─by over 50-fold, and a comparison of similar mutants enabled the identification of a site where mutations permit efficient hydroxylation. Findings suggest that the plasticity of terpene synthases enables an efficient sampling of structurally distinct starting points for building new functional molecules and provide an experimental framework for exploiting this plasticity in activity-guided screens.

Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase

Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.

Key difference between transition state stabilization and ground state destabilization: increasing atomic charge densities before or during enzyme–substrate binding

The origin of the enormous catalytic power of enzymes has been extensively studied through experimental and computational approaches. Although precise mechanisms are still subject to much debate, enzymes are thought to catalyze reactions by stabilizing transition states (TSs) or destabilizing ground states (GSs). By exploring the catalysis of various types of enzyme–substrate noncovalent interactions, we found that catalysis by TS stabilization and the catalysis by GS destabilization share common features by reducing the free energy barriers (ΔG^‡s) of reactions, but are different in attaining the requirement for ΔG^‡ reduction. Irrespective of whether enzymes catalyze reactions by TS stabilization or GS destabilization, they reduce ΔG^‡s by enhancing the charge densities of catalytic atoms that experience a reduction in charge density between GSs and TSs. Notably, in TS stabilization, the charge density of catalytic atoms is enhanced prior to enzyme–substrate binding; whereas in GS destabilization, the charge density of catalytic atoms is enhanced during the enzyme–substrate binding. Results show that TS stabilization and GS destabilization are not contradictory to each other and are consistent in reducing the ΔG^‡s of reactions. The full mechanism of enzyme catalysis includes the mechanism of reducing ΔG^‡ and the mechanism of enhancing atomic charge densities. Our findings may help resolve the debate between TS stabilization and GS destabilization and assist our understanding of catalysis and the design of artificial enzymes.

Transition state stabilization and ground state destabilization utilize the same molecular mechanism when lowering the free energy barriers (ΔG^‡s) of reactions, but differ in achieving the requirement for ΔG^‡ reduction.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

Engineering enzyme access tunnels

Enzymes are efficient and specific catalysts for many essential reactions in biotechnological and pharmaceutical industries. Many times, the natural enzymes do not display the catalytic efficiency, stability or specificity required for these industrial processes. The current enzyme engineering methods offer solutions to this problem, but they mainly target the buried active site where the chemical reaction takes place. Despite being many times ignored, the tunnels and channels connecting the environment with the active site are equally important for the catalytic properties of enzymes. Changes in the enzymatic tunnels and channels affect enzyme activity, specificity, promiscuity, enantioselectivity and stability. This review provides an overview of the emerging field of enzyme access tunnel engineering with case studies describing design of all the aforementioned properties. The software tools for the analysis of geometry and function of the enzymatic tunnels and channels and for the rational design of tunnel modifications will also be discussed. The combination of new software tools and enzyme engineering strategies will provide enzymes with access tunnels and channels specifically tailored for individual industrial processes.

Exploring Solanum tuberosum Epoxide Hydrolase Internal Architecture by Water Molecules Tracking

Several different approaches are used to describe the role of protein compartments and residues in catalysis and to identify key residues suitable for the modification of the activity or selectivity of the desired enzyme. In our research, we applied a combination of molecular dynamics simulations and a water tracking approach to describe the water accessible volume of Solanum tuberosum epoxide hydrolase. Using water as a molecular probe, we were able to identify small cavities linked with the active site: (i) one made up of conserved amino acids and indispensable for the proper positioning of catalytic water and (ii) two others in which modification can potentially contribute to enzyme selectivity and activity. Additionally, we identified regions suitable for de novo tunnel design that could also modify the catalytic properties of the enzyme. The identified hot-spots extend the list of the previously targeted residues used for modification of the regioselectivity of the enzyme. Finally, we have provided an example of a simple and elegant process for the detailed description of the network of cavities and tunnels, which can be used in the planning of enzyme modifications and can be easily adapted to the study of any other protein.

Protonation-Initiated Cyclization by a Class II Terpene Cyclase Assisted by Tunneling

Terpenes represent one of the most diversified classes of natural products with potent biological activities. The key to the myriad of polycyclic terpene skeletons with crucial functions in organisms from all kingdoms of life are terpene cyclase enzymes. These biocatalysts enable stereospecific cyclization of relatively simple, linear, prefolded polyisoprenes by highly complex, partially concerted, electrophilic cyclization cascades that remain incompletely understood. Herein, additional mechanistic light is shed on terpene biosynthesis by kinetic studies in mixed H2 O/D2 O buffers of a class II bacterial ent-copalyl diphosphate synthase. Mass spectrometry determination of the extent of deuterium incorporation in the bicyclic product, reminiscent of initial carbocation formation by protonation, resulted in a large kinetic isotope effect of up to seven. Kinetic analysis at different temperatures confirmed that the isotope effect was independent of temperature, which is consistent with hydrogen tunneling.

MD Simulations Reveal Complex Water Paths in Squalene-Hopene Cyclase: Tunnel-Obstructing Mutations Increase the Flow of Water in the Active Site

Squalene-hopene cyclase catalyzes the cyclization of squalene to hopanoids. A previous study has identified a network of tunnels in the protein, where water molecules have been indicated to move. Blocking these tunnels by site-directed mutagenesis was found to change the activation entropy of the catalytic reaction from positive to negative with a concomitant lowering of the activation enthalpy. As a consequence, some variants are faster and others are slower than the wild type (wt) in vitro under optimal reaction conditions for the wt. In this study, molecular dynamics (MD) simulations have been performed for the wt and the variants to investigate how the mutations affect the protein structure and the water flow in the enzyme, hypothetically influencing the activation parameters. Interestingly, the tunnel-obstructing variants are associated with an increased flow of water in the active site, particularly close to the catalytic residue Asp376. MD simulations with the substrate present in the active site indicate that the distance for the rate-determining proton transfer between Asp376 and the substrate is longer in the tunnel-obstructing protein variants than in the wt. On the basis of the previous experimental results and the current MD results, we propose that the tunnel-obstructing variants, at least partly, could operate by a different catalytic mechanism, where the proton transfer may have contributions from a Grotthuss-like mechanism.

Enzymatic Addition of Alcohols to Terpenes by Squalene Hopene Cyclase Variants

Squalene-hopene cyclases (SHCs) catalyze the polycyclization of squalene into a mixture of hopene and hopanol. Recently, amino-acid residues lining the catalytic cavity of the SHC from Alicyclobacillus acidocaldarius were replaced by small and large hydrophobic amino acids. The alteration of leucine 607 to phenylalanine resulted in increased enzymatic activity towards the formation of an intermolecular farnesyl-farnesyl ether product from farnesol. Furthermore, the addition of small-chain alcohols acting as nucleophiles led to the formation of non-natural ether-linked terpenoids and, thus, to significant alteration of the product pattern relative to that obtained with the wild type. It is proposed that the mutation of leucine at position 607 may facilitate premature quenching of the intermediate by small alcohol nucleophiles. This mutagenesis-based study opens the field for further intermolecular bond-forming reactions and the generation of non-natural products.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Separating Thermodynamics from Kinetics-A New Understanding of the Transketolase Reaction

Transketolase catalyzes asymmetric C−C bond formation of two highly polar compounds. Over the last 30 years, the reaction has unanimously been described in literature as irreversible because of the concomitant release of CO₂ if using lithium hydroxypyruvate (LiHPA) as a substrate. Following the reaction over a longer period of time however, we have now found it to be initially kinetically controlled. Contrary to previous suggestions, for the non‐natural conversion of synthetically more interesting apolar substrates, the complete change of active‐site polarity is therefore not necessary. From docking studies it was revealed that water and hydrogen‐bond networks are essential for substrate binding, thus allowing aliphatic aldehydes to be converted in the charged active site of transketolase.

Squalene-hopene cyclases-evolution, dynamics and catalytic scope

Herein we highlight recent mechanistic findings on the impact of solvent dynamics on catalysis displayed by squalene-hopene cyclases (SHCs). These fascinating biocatalysts that appeared early during the evolution of terpene biosynthetic machineries exploit a catalytic aspartic acid donating the anti-oriented proton to the terminal CC double bond of pre-folded isoprenoid substrates. We review how the unusual strength of this Brønsted acid can be used to harness a plethora of non-natural protonation-driven reactions in a plastic enzyme fold. Moreover, recent results underline how the reaction termination by deprotonation or water addition is governed by the spatial location of water in the active site. Site-directed mutagenesis of amino acids located in the hydrophobic binding pocket allows for the generation of novel catalytic function by active site reshaping with relatively small enzyme libraries. A deepened understanding of triterpene cyclase dynamics in concert with chemical expertise thus have a great potential to allow for the biocatalytic manufacturing of tailored building bricks that would expand the chemical repertoire currently found in nature.

Unraveling Entropic Rate Acceleration Induced by Solvent Dynamics in Membrane Enzymes

Enzyme catalysis evolved in an aqueous environment. The influence of solvent dynamics on catalysis is, however, currently poorly understood and usually neglected. The study of water dynamics in enzymes and the associated thermodynamical consequences is highly complex and has involved computer simulations, nuclear magnetic resonance (NMR) experiments, and calorimetry. Water tunnels that connect the active site with the surrounding solvent are key to solvent displacement and dynamics. The protocol herein allows for the engineering of these motifs for water transport, which affects specificity, activity and thermodynamics. By providing a biophysical framework founded on theory and experiments, the method presented herein can be used by researchers without previous expertise in computer modeling or biophysical chemistry. The method will advance our understanding of enzyme catalysis on the molecular level by measuring the enthalpic and entropic changes associated with catalysis by enzyme variants with obstructed water tunnels. The protocol can be used for the study of membrane-bound enzymes and other complex systems. This will enhance our understanding of the importance of solvent reorganization in catalysis as well as provide new catalytic strategies in protein design and engineering.

Exploring water as building bricks in enzyme engineering

A de novo designed water pattern is used to achieve a 34-fold accelerated promiscuous enzyme catalysis by efficient transition state stabilization.