Influence of Water and Enzyme on the Post-Transition State Bifurcation of NgnD-Catalyzed Ambimodal [6+4]/[4+2] Cycloaddition

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.

Dynamically or Kinetically Controlled? Computational Study of the Mechanisms of Electrophilic Aminoalkenylation of Heteroaromatics with Keteniminium Ions

Quantum chemical calculations and molecular dynamics simulations were applied to study the electrophilic aminoalkenylation of heteroaromatics with keniminium ions. Post-transition state bifurcation (PTSB) was found in the electrophilic addition step for the aminoalkenylation of pyrroles and indoles, and the selectivity for these reactions was dynamically controlled. However, the aminoalkenylation of furan was kinetically controlled because no apparent PTSB was found in the electrophilic addition step. The substituents on the keteniminium ions can also affect the dynamic results for the aminoalkenylations to pyrroles: the C2-aminoalkenylated product is much more favored over the C3-aminoalkenylated product for keteniminium ions with electron-donating substituents, while the product ratio (C2 product/C3 product) decreased when stronger electron-withdrawing substituents were applied.

Exploring the mechanism of the PTP1B inhibitors by molecular dynamics and experimental study

Protein tyrosine phosphatase 1B (PTP1B) has proven to be an attractive target for the treatment of cancer, diabetes and other diseases. Although many PTP1B inhibitors with various scaffolds have been developed, there is still a lack of PTP1B inhibitor with high specificity and acceptable pharmacological properties. Therefore, it is urgent to develop more methods to explore complex action mode of PTP1B and ligands for designing ideal PTP1B modulators. In this work, we developed a potential molecular dynamics (MD) analytic mode to analyze the mechanism of active compounds 6a and 6e against PTP1B from different perspectives, including the stable ability, interactions and binding site of ligand and protein, the binding energy, relative movement between residues and changes in protein internal interactions. The simulated results demonstrated that compound 6a bound more stably to the active pocket of PTP1B than 6e due to its smaller molecular volume (326 Å3), matched electronegativity, and enhanced the positive correlation motion of residues, especially for WPD loop and P loop. Lastly, compound 6a as a competitive inhibitor for PTP1B was verified by enzyme kinetic assay. This work successfully studied the mechanism of compound 6a against PTP1B from various aspects, enriched the analysis of interaction mode between PTP1B and inhibitors. In summary, we hope that this work could provide more theoretical information for designing and developing more novel and ideal PTP1B inhibitors in the future.

Symmetric Post-Transition State Bifurcation Reactions with Berry Pseudomagnetic Fields

We investigate how the Berry force (i.e., the pseudomagnetic force operating on nuclei as induced by electronic degeneracy and spin-orbit coupling (SOC)) might modify a post-transition state bifurcation (PTSB) reaction path and affect product selectivity for situations when multiple products share the same transition state. To estimate the magnitude of this effect, Langevin dynamics are performed on a model system with a valley-ridge inflection (VRI) point in the presence of a magnetic field (that mimics the Berry curvature). We also develop an analytic model for such selectivity that depends on key parameters such as the surface topology, the magnitude of the Berry force, and the nuclear friction. Within this dynamical model, static electronic structure calculations (at the level of generalized Hartree-Fock with spin-orbit coupling (GHF+SOC) theory) suggest that electronic spin induced Berry force effects may indeed lead to noticeable changes in methoxy radical isomerization.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

Mechanism of the Stereoselective Catalysis of Diels-Alderase PyrE3 Involved in Pyrroindomycin Biosynthesis

The biosynthesis of pyrroindomycins A and B features a complexity-building [4 + 2] cycloaddition cascade, which generates the spirotetramate core under the catalytic effects of monofunctional Diels-Alderases PyrE3 and PyrI4. We recently showed that the main functions of PyrI4 include acid catalysis and induced-fit/conformational selection. We now present quantum mechanical and molecular dynamics studies implicating a different mode of action by PyrE3, which prearranges an anionic polyene substrate into a high-energy reactive conformation at which an inverse-electron-demand Diels-Alder reaction can occur with a low barrier. Stereoselection is realized by strong binding interactions at the endo stereochemical relationship and a local steric constraint on the endo-1,3-diene unit. These findings, illustrating distinct mechanisms for PyrE3 and PyrI4, highlight how nature has evolved multiple ways to catalyze Diels-Alder reactions.

Kinetic Analysis of a Reaction Path Network Including Ambimodal Transition States: A Case Study of an Intramolecular Diels-Alder Reaction

This study proposes a methodology for the kinetic analysis of a reaction path network including ambimodal transition states (TSs), through which an ensemble of trajectories bifurcates to multiple minima in a phenomenon called dynamical bifurcation. The proposed methodology consists of three techniques: an automated reaction path search to construct a reaction path network including ambimodal TSs, an ab initio molecular dynamics simulation to evaluate the branching ratio, and the definition of rate constants incorporating this ratio. Applying the procedure to a Diels-Alder reaction, it was found that the inclusion of dynamical bifurcations is necessary to explain the experimental reaction yield of a byproduct. In addition, it was verified that the products take 10¹³ s to reach thermal equilibrium and that the experimental selectivity is determined by the dynamical bifurcations.