Enhancing luciferase activity and stability through generative modeling of natural enzyme sequences.

The availability of natural protein sequences synergized with generative AI provides new paradigms to engineer enzymes. Although active enzyme variants with numerous mutations have been designed using generative models, their performance often falls short of their wild type counterparts. Additionally, in practical applications, choosing fewer mutations that can rival the efficacy of extensive sequence alterations is usually more advantageous. Pinpointing beneficial single mutations continues to be a formidable task. In this study, using the generative maximum entropy model to analyze Renilla luciferase (RLuc) homologs, and in conjunction with biochemistry experiments, we demonstrated that natural evolutionary information could be used to predictively improve enzyme activity and stability by engineering the active center and protein scaffold, respectively. The success rate to improve either luciferase activity or stability of designed single mutants is ~50%. This finding highlights nature's ingenious approach to evolving proficient enzymes, wherein diverse evolutionary pressures are preferentially applied to distinct regions of the enzyme, ultimately culminating in an overall high performance. We also reveal an evolutionary preference in RLuc toward emitting blue light that holds advantages in terms of water penetration compared to other light spectra. Taken together, our approach facilitates navigation through enzyme sequence space and offers effective strategies for computer-aided rational enzyme engineering.

Exploring the Light-Emitting Agents in Renilla Luciferases by an Effective QM/MM Approach.

Bioluminescence is a fascinating natural phenomenon, wherein organisms produce light through specific biochemical reactions. Among these organisms, Renilla luciferase (RLuc) derived from the sea pansy Renilla reniformis is notable for its blue light emission and has potential applications in bioluminescent tagging. Our study focuses on RLuc8, a variant of RLuc with eight amino acid substitutions. Recent studies have shown that the luminescent emitter coelenteramide can adopt multiple protonation states, which may be influenced by nearby residues at the enzyme's active site, demonstrating a complex interplay between protein structure and bioluminescence. Herein, using the quantum mechanical consistent force field method and the semimacroscopic protein dipole-Langevin dipole method with linear response approximation, we show that the phenolate state of coelenteramide in RLuc8 is the primary light-emitting species in agreement with experimental results. Our calculations also suggest that the proton transfer (PT) from neutral coelenteramide to Asp162 plays a crucial role in the bioluminescence process. Additionally, we reproduced the observed emission maximum for the amide anion in RLuc8-D120A and the pyrazine anion in the presence of a Na+ counterion in RLuc8-D162A, suggesting that these are the primary emitters. Furthermore, our calculations on the neutral emitter in the engineered AncFT-D160A enzyme, structurally akin to RLuc8-D162A but with a considerably blue-shifted emission peak, aligned with the observed data, possibly explaining the variance in emission peaks. Overall, this study demonstrates an effective approach to investigate chromophores' bimolecular states while incorporating the PT process in emission spectra calculations, contributing valuable insights for future studies of PT in photoproteins.

Reactivity and Enantioselectivity in NHC Organocatalysis Provide Evidence for the Complex Role of Modifications at the Secondary Sphere.

Secondary-sphere interactions are often harnessed to control reactivity and selectivity in organometallic and enzymatic catalysis. Yet, such strategies have only recently been explicitly applied in the context of organocatalytic systems. Although increased stability, reproducibility, and selectivity were obtained in previous work using this approach, the precise mechanistic pathway promoted by secondary-sphere modification in organocatalysis remained unclear. Herein, we report a comprehensive mechanistic study on the origin of the unique reactivity patterns and stereocontrol observed with boronic acids (BAs) as secondary-sphere modifiers of N-heterocyclic carbene (NHC) organocatalysts. Kinetic experiments revealed partial order in catalyst upon the addition of BA and unusual preactivation behavior, indicating the presence of stable off-cycle catalyst aggregation and BA-base adducts. These hypotheses were supported both by computations and by a series of NMR and nonlinear effect experiments. Furthermore, computations indicated a rate-limiting, water-assisted hydrogen atom transfer mechanism. This finding led to a considerable enhancement in the experimental reaction rate while maintaining excellent enantioselectivity by adding catalytic amounts of water. Finally, computations and racemization experiments uncovered an uncommon Curtin-Hammett-controlled enantioselectivity in the presence of secondary-sphere modifiers.

Thioxobimanes.

Dioxobimanes, colloquially known as bimanes, are a well-established family of N-heterobicyclic compounds that share a characteristic core structure, 1,5-diazabicyclo[3.3.0]octadienedione, bearing two endocyclic carbonyl groups. By sequentially thionating these carbonyls in the syn and anti isomers of the known (Me,Me)dioxobimane, we were able to synthesize a series of thioxobimanes, representing the first heavy-chalcogenide bimane variants. These new compounds were extensively characterized spectroscopically and crystallographically, and their aromaticity was probed computationally. Their potential role as ligands for transition metals was demonstrated by synthesizing a representative gold(I)-thioxobimane complex.

Thermally-Activated Tunneling in the Two-Water Bridge Catalyzed Tautomerization of Phosphinylidene Compounds.

Phosphinylidenes are an important class of organophosphorus compounds that can exhibit tautomerization between tricoordinated P(III) hydroxide (R1 R2 POH) and a pentacoordinated P(V) oxide (R1 R2 P(O)H) form. Herein we show, using the canonical variational transition state theory combined with multidimensional small-curvature tunneling approximation, the dominance of proton tunneling in the two-water-bridged tautomerizations of phosphinous acid and model phosphinylidenes comprising phosphosphinates, H-phosphonates, H-phosphinates and secondary phosphine oxides. Based on the studied system, the contribution of thermally-activated tunneling is predicted to speed up the semiclassical reaction rate by ca. threefold to as large as two orders of magnitude at 298.15 K in the gas phase. The large KIE and the concavity in the Arrhenius plots are further fingerprints of tunneling. The simulations also predicted that the rapid tunneling rate and short half-life span for the forward reaction, as opposed to the reverse reaction in fluorinated secondary phosphine oxides, would result in P(V) being elusive and only P(III) being isolable, which agrees with previous experiments where only P(III) was detected by IR and NMR spectroscopy. We also explored the role of solvent and predicted tunneling to be substantial.

Boron Tunneling in the "Weak" Bond-Stretch Isomerization of N-B Lewis Adducts.

Some nitrile-boron halide adducts exhibit a double-well potential energy surface with two distinct minima: a "long bond" geometry (LB, a van der Waals interaction mostly based on electrostatics, but including a residual charge transfer component) and a "short bond" structure (SB, a covalent dative bond). This behavior can be considered as a "weak" form of bond stretch isomerism. Our computations reveal that complexes RCN-BX3 (R=CH3 , FCH2 , BrCH2 , and X=Cl, Br) exhibit a fast interconversion from LB to SB geometries even close to the absolute zero thanks to a boron atom tunneling mechanism. The computed half-lives of the meta-stable LB compounds vary between minutes to nanoseconds at cryogenic conditions. Accordingly, we predict that the long bond structures are practically impossible to isolate or characterize, which agrees with previous matrix-isolation experiments.

History and Future of Dative Bonds.

Dative Bond (IUPAC): "The coordination bond formed upon interaction between molecular species, one of which serves as a donor and the other as an acceptor of the electron pair to be shared in the complex formed… The distinctive feature of dative bonds is that their minimum-energy rupture in the gas phase or in inert solvent follows the heterolytic bond cleavage path." This definition encompasses an immense number of molecules such as Lewis adducts, transition-metal complexes, supposedly hypervalent or hypovalent systems, and many molecules with multifaceted Lewis structures. Still, there is a large reticence to include dative bonds in the regular depiction of molecules, and even a larger unawareness of the dative bond arrow in many chemical circles. Herein we will discuss in simple chemical terms the past, present and future of such bonds. In addition, we will try to provide cleaner options to represent intricate molecules without sacrificing physical accuracy.

Carbon Tunneling in the Automerization of Cyclo[18]carbon.

Cyclo[18]carbon (C18 ), a recently synthesized carbon allotrope, was found to have a polyynic ground-state structure with D9h symmetry and formally alternating single and triple bonds. Yet, under less influencing experimental conditions this molecule might undergo an automerization reaction between its two degenerate geometries through a cumulenic (non-alternating, adjacent double bonds) D18h transition state. Herein, we discuss the role of quantum mechanical tunneling (QMT) in this degenerate reaction. Our computations predict that at the experimental temperature (5 K) the reaction in the gas phase is completely driven by an extremely rapid heavy atom tunneling (k=2.1×108  s-1 ). Even when approaching room temperature, the QMT rate is still an order of magnitude faster than the semi-classical one. We propose an experimental test to support our prediction, by measuring a characteristic tunneling energy splitting within the radio wave region. Additionally, we examine the role of QMT in other hypothetical C4n+2 carbon clusters.

Ping-Pong Tunneling Reactions: Can Fluoride Jump at Absolute Zero?

In a recent study, Scheiner designed a double-germanium-based fluoride receptor that binds the halogen by means of strong tetrel bonds (Chem. Eur. J. 2016, 22, 18850). In this system the F- binds to the germanium atoms in an asymmetric fashion, thereby producing a double-well potential in which the fluoride can jump from one germanium to the other as in a ping-pong game. Herein we prove through the use of computational tools that at cryogenic temperatures this rearrangement occurs by heavy-atom quantum mechanical tunneling. The inductive strength of the substituents and the polarity of the solvent can modify the barrier and the tunneling rate. But the strongest effect is observed upon modification of the geometry of the molecule by specific substitutions that affect the barrier width, the most critical factor in a tunneling mechanism. We postulate two experimental tests, one by microwave spectroscopy and one by cryogenic NMR spectroscopy, that can prove the predicted fluoride tunneling.

Isotope-Controlled Selectivity by Quantum Tunneling: Hydrogen Migration versus Ring Expansion in Cyclopropylmethylcarbenes.

Using the tunneling-controlled reactivity of cyclopropylmethylcarbene, we demonstrate the viability of isotope-controlled selectivity (ICS), a novel control element of chemical reactivity where a molecular system with two conceivable products of tunneling exclusively produces one or the other, depending only on isotopic composition. Our multidimensional small-curvature tunneling (SCT) computations indicate that, under cryogenic conditions, 1-methoxycyclopropylmethylcarbene shows rapid H-migration to 1-methoxy-1-vinylcyclopropane, whereas deuterium-substituted 1-methoxycyclopropyl-d3-methylcarbene undergoes ring expansion to 1-d3-methylcyclobutene. This predicted change in reactivity constitutes the first example of a kinetic isotope effect that discriminates between the formation of two products.

Assessing the Catalytic Role of Native Glucagon Amyloid Fibrils.

Glucagon stands out as a pivotal peptide hormone, instrumental in controlling blood glucose levels and lipid metabolism. While the formation of glucagon amyloid fibrils has been documented, their biological functions remain enigmatic. Recently, we demonstrated experimentally that glucagon amyloid fibrils can act as catalysts in several biological reactions including esterolysis, lipid hydrolysis, and dephosphorylation. Herein, we present a multiscale quantum mechanics/molecular mechanics (QM/MM) simulation of the acylation step in the esterolysis of para-nitrophenyl acetate (p-NPA), catalyzed by native glucagon amyloid fibrils, serving as a model system to elucidate their catalytic function. This step entails a concerted mechanism, involving proton transfer from serine to histidine, followed by the nucleophilic attack of the serine oxy anion on the carbonyl carbon of p-NPA. We computed the binding energy and free-energy profiles of this reaction using the protein-dipole Langevin-dipole (PDLD) within the linear response approximation (LRA) framework (PDLD/S-LRA-2000) and the empirical valence bond (EVB) methods. This included simulations of the reaction in an aqueous environment and in the fibril, enabling us to estimate the catalytic effect of the fibril. Our EVB calculations obtained a barrier of 23.4 kcal mol-1 for the enzyme-catalyzed reaction compared to the experimental value of 21.9 kcal mol-1 (and a calculated catalytic effect of 3.2 kcal mol-1 compared to the observed effect of 4.7 kcal mol-1). This close agreement together with the barrier reduction when transitioning from the reference solution reaction to the amyloid fibril provides supporting evidence to the catalytic role of glucagon amyloid fibrils. Moreover, employing the PDLD/S-LRA-2000 approach further reinforced exclusively the enzyme's catalytic role. The results presented in this study contribute significantly to our understanding of the catalytic role of glucagon amyloid fibrils, marking, to the best of our knowledge, the first-principles mechanistic investigation of fibrils using QM/MM methods. Therefore, our findings offer fruitful insights for future research into the mechanisms of related amyloid catalysis.

Enhancing Luciferase Activity and Stability through Generative Modeling of Natural Enzyme Sequences.

The availability of natural protein sequences synergized with generative artificial intelligence (AI) provides new paradigms to create enzymes. Although active enzyme variants with numerous mutations have been produced using generative models, their performance often falls short compared to their wild-type counterparts. Additionally, in practical applications, choosing fewer mutations that can rival the efficacy of extensive sequence alterations is usually more advantageous. Pinpointing beneficial single mutations continues to be a formidable task. In this study, using the generative maximum entropy model to analyze Renilla luciferase homologs, and in conjunction with biochemistry experiments, we demonstrated that natural evolutionary information could be used to predictively improve enzyme activity and stability by engineering the active center and protein scaffold, respectively. The success rate of designed single mutants is ~50% to improve either luciferase activity or stability. These finding highlights nature's ingenious approach to evolving proficient enzymes, wherein diverse evolutionary pressures are preferentially applied to distinct regions of the enzyme, ultimately culminating in an overall high performance. We also reveal an evolutionary preference in Renilla luciferase towards emitting blue light that holds advantages in terms of water penetration compared to other light spectrum. Taken together, our approach facilitates navigation through enzyme sequence space and offers effective strategies for computer-aided rational enzyme engineering.

Heavy-Atom Tunneling in the Covalent/Dative Bond Complexation of Cyclo[18]carbon-Piperidine.

Recent quantum chemical computations demonstrated the electron-acceptance behavior of this highly reactive cyclo[18]carbon (C18) ring with piperidine (pip). The C18-pip complexation exhibited a double-well potential along the N-C reaction coordinate, forming a van der Waals (vdW) adduct and a more stable, strong covalent/dative bond (DB) complex by overcoming a low activation barrier. By means of direct dynamical computations using canonical variational transition state theory (CVT), including the small-curvature tunneling (SCT), we show the conspicuous role of heavy atom quantum mechanical tunneling (QMT) in the transformation of vdW to DB complex in the solvent phase near absolute zero. Below 50 K, the reaction is entirely driven by QMT, while at 30 K, the QMT rate is too rapid (kT ∼ 0.02 s-1), corresponding to a half-life time of 38 s, indicating that the vdW adduct will have a fleeting existence. We also explored the QMT rates of other cyclo[n]carbon-pip systems. This study sheds light on the decisive role of QMT in the covalent/DB formation of the C18-pip complex at cryogenic temperatures.