Competition between N and O: use of diazine N-oxides as a test case for the Marcus theory rationale for ambident reactivity

The preferred site of alkylation of diazine N-oxides by representative hard and soft alkylating agents was established conclusively using the ¹H–¹⁵N HMBC NMR technique in combination with other NMR spectroscopic methods. Alkylation of pyrazine N-oxides (1 and 2) occurs preferentially on nitrogen regardless of the alkylating agent employed, while O-methylation of pyrimidine N-oxide (3) is favoured in its reaction with MeOTf. As these outcomes cannot be explained in the context of the hard/soft acid/base (HSAB) principle, we have instead turned to Marcus theory to rationalise these results. Marcus intrinsic barriers (ΔG^‡₀) and Δ_rG° values were calculated at the DLPNO-CCSD(T)/def2-TZVPPD/SMD//M06-2X-D3/6-311+G(d,p)/SMD level of theory for methylation reactions of 1 and 3 by MeI and MeOTf, and used to derive Gibbs energies of activation (ΔG^‡) for the processes of N- and O-methylation, respectively. These values, as well as those derived directly from the DFT calculations, closely reproduce the observed experimental N- vs. O-alkylation selectivities for methylation reactions of 1 and 3, indicating that Marcus theory can be used in a semi-quantitative manner to understand how the activation barriers for these reactions are constructed. It was found that N-alkylation of 1 is favoured due to the dominant contribution of Δ_rG° to the activation barrier in this case, while O-alkylation of 3 is favoured due to the dominant contribution of the intrinsic barrier (ΔG^‡₀) for this process. These results are of profound significance in understanding the outcomes of reactions of ambident reactants in general.

Marcus theory enables rationalisation and quantification of selectivities in reactions of ambident nucleophiles for which the HSAB principle cannot operate.     

Synthesis of a Ni Complex Chelated by a [2.2]Paracyclophane-Functionalized Diimine Ligand and Its Catalytic Activity for Olefin Oligomerization

A diimine ligand having two [2.2]paracyclophanyl substituents at the N atoms (L1) was prepared from the reaction of amino[2.2]paracyclophane with acenaphtenequinone. The ligand reacts with NiBr₂(dme) (dme: 1,2-dimethoxyethane) to form the dibromonickel complex with (R,R) and (S,S) configuration, NiBr₂(L1). The structure of the complex was confirmed by X-ray crystallography. NiBr₂(L1) catalyzes oligomerization of ethylene in the presence of methylaluminoxane (MAO) co-catalyst at 10–50 °C to form a mixture of 1- and 2-butenes after 3 h. The reactions for 6 h and 8 h at 25 °C causes further increase of 2-butene formed via isomerization of 1-butene and formation of hexenes. Reaction of 1-hexene catalyzed by NiBr₂(L1)–MAO produces 2-hexene via isomerization and C12 and C18 hydrocarbons via oligomerization. Consumption of 1-hexene of the reaction obeys first-order kinetics. The kinetic parameters were obtained to be ΔG^‡ = 93.6 kJ mol⁻¹, ΔH^‡ = 63.0 kJ mol⁻¹, and ΔS^‡ = −112 J mol⁻¹deg⁻¹. NiBr₂(L1) catalyzes co-dimerization of ethylene and 1-hexene to form C8 hydrocarbons with higher rate and selectivity than the tetramerization of ethylene.     

Computational study of the mechanism and selectivity of ruthenium-catalyzed hydroamidations of terminal alkynes

Density functional theory calculations were performed to elucidate the mechanism of the ruthenium-catalyzed hydroamidation of terminal alkynes, a powerful and sustainable method for the stereoselective synthesis of enamides. The results provide an explanation for the puzzling experimental finding that with tri-n-butylphosphine (P(Bu)₃) as the ligand, the E-configured enamides are obtained, whereas the stereoselectivity is inverted in favor of the Z-configured enamides with (dicyclohexylphosphino)methane (dcypm) ligands. Using the addition of pyrrolidinone to 1-hexyne as a model reaction, various pathways were investigated, among which a catalytic cycle turned out to be most advantageous for both ligand systems that consists of: (a) oxidative addition, (b) alkyne coordination, (c) alkyne insertion (d) vinyl-vinylidene rearrangement, (e) nucleophilic transfer and finally (f) reductive elimination. The stereoselectivity of the reaction is decided in the nucleophilic transfer step. For the P(ⁿBu)₃ ligand, the butyl moiety is oriented anti to the incoming 2-pyrolidinyl unit during the nucleophilic transfer step, whereas for the dcypm ligand, steric repulsion between the butyl and cyclohexyl groups turns it into a syn orientation. Overall, the formation of E-configured product is favorable by 4.8 kcal mol^–1 (Δ^‡ GSDL) for the catalytic cycle computed with P(Bu)₃ as ancillary ligand, whereas for the catalytic cycle computed with dcypm ligands, the Z-product is favored by 7.0 kcal mol^–1 (Δ^‡ GSDL). These calculations are in excellent agreement with experimental findings.     

Factors determining the enzyme catalytic power caused by noncovalent interactions: Charge alterations in enzyme active sites

Understanding the origin of the enormous catalytic power of enzymes is very important. Electrostatic interactions and desolvation are the phenomena that are most proposed to explain the catalysis of enzymes; however, they also decelerate enzymatic reactions. How enzymes catalyze reactions through noncovalent interactions is still not well-understood. In this study, we explored how enzyme-substrate noncovalent interactions affect the free energy barriers (ΔG³s) of reactions by using a theoretical derivation approach. We found that enzymes reduce ΔG³s of reactions by decreasing positive charges and/or increasing negative charges in the electron-donating centers and by decreasing negative charges and/or increasing positive charges in the electron-accepting centers of reactions. Enzyme-substrate noncovalent interactions are essential approaches through which the charge alterations lead to ΔG³ reductions. Validations with reported experimental data demonstrated that this charge alteration mechanism can explain the catalyses caused by diverse types of noncovalent interactions. Electrostatic interactions and desolvation are the most observed noncovalent interactions essential for ΔG³ reductions. This mechanism does not contradict any specific enzymatic catalysis and overcomes the shortages of the electrostatic interaction and desolvation mechanisms. This study can provide useful guidance in exploring enzymatic catalysis and designing catalyst.     

Strategies for switching the mechanism of proton-coupled electron transfer reactions illustrated by mechanistic zone diagrams

The mechanism by which proton-coupled electron transfer (PCET) occurs is of fundamental importance and has great consequences for applications, e.g. in catalysis. However, determination and tuning of the PCET mechanism is often non-trivial. Here, we apply mechanistic zone diagrams to illustrate the competition between concerted and stepwise PCET-mechanisms in the oxidation of 4-methoxyphenol by Ru(bpy)₃³⁺-derivatives in the presence of substituted pyridine bases. These diagrams show the dominating mechanism as a function of driving force for electron and proton transfer (ΔG⁰_ET and ΔG⁰_PT) respectively [Tyburski et al., J. Am. Chem. Soc., 2021, 143, 560]. Within this framework, we demonstrate strategies for mechanistic tuning, namely balancing of ΔG⁰_ET and ΔG⁰_PT, steric hindrance of the proton-transfer coordinate, and isotope substitution. Sterically hindered pyridine bases gave larger reorganization energy for concerted PCET, resulting in a shift towards a step-wise electron first-mechanism in the zone diagrams. For cases when sufficiently strong oxidants are used, substitution of protons for deuterons leads to a switch from concerted electron–proton transfer (CEPT) to an electron transfer limited (ETPT_lim) mechanism. We thereby, for the first time, provide direct experimental evidence, that the vibronic coupling strength affects the switching point between CEPT and ETPT_lim, i.e. at what driving force one or the other mechanism starts dominating. Implications for solar fuel catalysis are discussed.

The mechanism by which proton-coupled electron transfer (PCET) occurs is of fundamental importance and has great consequences for applications, e.g. in catalysis.     

Enantioselective synthesis of chiral porphyrin macrocyclic hosts and kinetic enantiorecognition of viologen guests

The construction of macromolecular hosts that are able to thread chiral guests in a stereoselective fashion is a big challenge. We herein describe the asymmetric synthesis of two enantiomeric C₂-symmetric porphyrin macrocyclic hosts that thread and bind different viologen guests. Time-resolved fluorescence studies show that these hosts display a factor 3 kinetic preference (ΔΔG^‡_on = 3 kJ mol⁻¹) for threading onto the different enantiomers of a viologen guest appended with bulky chiral 1-phenylethoxy termini. A smaller kinetic selectivity (ΔΔG^‡_on = 1 kJ mol⁻¹) is observed for viologens equipped with small chiral sec-butoxy termini. Kinetic selectivity is absent when the C₂-symmetric hosts are threaded onto chiral viologens appended with chiral tails in which the chiral moieties are located in the centers of the chains, rather than at the chain termini. The reason is that the termini of the latter guests, which engage in the initial stages of the threading process (entron effect), cannot discriminate because they are achiral, in contrast to the chiral termini of the former guests. Finally, our experiments show that the threading and de-threading rates are balanced in such a way that the observed binding constants are highly similar for all the investigated host–guest complexes, i.e. there is no thermodynamic selectivity.

Chiral guests display kinetic stereoselective threading through chiral porphyrin cages if their chirality is located at the chain ends and not in the centers, supporting the previously reported entron effect of threading.     

Catalyst control of selectivity in the C–O bond alumination of biomass derived furans

Non-catalysed and catalysed reactions of aluminium reagents with furans, dihydrofurans and dihydropyrans were investigated and lead to ring-expanded products due to the insertion of the aluminium reagent into a C–O bond of the heterocycle. Specifically, the reaction of [{(ArNCMe)₂CH}Al] (Ar = 2,6-di-iso-propylphenyl, 1) with furans proceeded between 25 and 80 °C leading to dearomatised products due to the net transformation of a sp² C–O bond into a sp² C–Al bond. The kinetics of the reaction of 1 with furan were found to be 1st order with respect to 1 with activation parameters ΔH^‡ = +19.7 (±2.7) kcal mol⁻¹, ΔS^‡ = −18.8 (±7.8) cal K⁻¹ mol⁻¹ and ΔG^‡_{298 K} = +25.3 (±0.5) kcal mol⁻¹ and a KIE of 1.0 ± 0.1. DFT calculations support a stepwise mechanism involving an initial (4 + 1) cycloaddition of 1 with furan to form a bicyclic intermediate that rearranges by an α-migration. The selectivity of ring-expansion is influenced by factors that weaken the sp² C–O bond through population of the σ*-orbital. Inclusion of [Pd(PCy₃)₂] as a catalyst in these reactions results in expansion of the substrate scope to include 2,3-dihydrofurans and 3,4-dihydropyrans and improves selectivity. Under catalysed conditions, the C–O bond that breaks is that adjacent to the sp²C–H bond. The aluminium(iii) dihydride reagent [{(MesNCMe)₂CH}AlH₂] (Mes = 2,4,6-trimethylphenyl, 2) can also be used under catalytic conditions to effect a dehydrogenative ring-expansion of furans. Further mechanistic analysis shows that C–O bond functionalisation occurs via an initial C–H bond alumination. Kinetic products can be isolated that are derived from installation of the aluminium reagent at the 2-position of the heterocycle. C–H alumination occurs with a KIE of 4.8 ± 0.3 consistent with a turnover limiting step involving oxidative addition of the C–H bond to the palladium catalyst. Isomerisation of the kinetic C–H aluminated product to the thermodynamic C–O ring expansion product is an intramolecular process that is again catalysed by [Pd(PCy₃)₂]. DFT calculations suggest that the key C–O bond breaking step involves attack of an aluminium based metalloligand on the 2-palladated heterocycle. The new methodology has been applied to important platform chemicals from biomass.

Non-catalysed and catalysed reactions of aluminium reagents with furans, dihydrofurans and dihydropyrans were investigated and lead to ring-expanded products due to the insertion of the aluminium reagent into a C–O bond of the heterocycle.     

Butatrienes as extended alkenes: barriers to internal rotation and substitution effects on the stabilities of the ground states and transition states

The barriers to internal rotation of methylated, ethynylated, and vinylated butatrienes and alkenes were calculated at the CASPT2/6-31G(d)//B3LYP/6-31G(d) level. Calculated butatriene rotational barriers are lower than those of analogous alkenes, but there is a larger variance in rotational barrier for alkenes than for butatrienes. The barriers to rotation were analyzed by isodesmic equations designed to estimate the substituent effects in the ground (GS) and transition (TS) states individually. The GSs of both series are stabilized to roughly the same extent. In contrast, the TSs of butatrienes are more stabilized overall than those of alkenes. Much of the stabilization in the TS of butatrienes comes from the internal triple bond and not from the substituent. Estimation of the substituent stabilization alone reveals the TSs of ethylenes to be more stabilized by substitution than butatrienes.     

Glycoside hydrolase stabilization of transition state charge: new directions for inhibitor design

Carbasugars are structural mimics of naturally occurring carbohydrates that can interact with and inhibit enzymes involved in carbohydrate processing. In particular, carbasugars have attracted attention as inhibitors of glycoside hydrolases (GHs) and as therapeutic leads in several disease areas. However, it is unclear how the carbasugars are recognized and processed by GHs. Here, we report the synthesis of three carbasugar isotopologues and provide a detailed transition state (TS) analysis for the formation of the initial GH-carbasugar covalent intermediate, as well as for hydrolysis of this intermediate, using a combination of experimentally measured kinetic isotope effects and hybrid QM/MM calculations. We find that the α-galactosidase from Thermotoga maritima effectively stabilizes TS charge development on a remote C5-allylic center acting in concert with the reacting carbasugar, and catalysis proceeds via an exploded, or loose, S_N2 transition state with no discrete enzyme-bound cationic intermediate. We conclude that, in complement to what we know about the TS structures of enzyme-natural substrate complexes, knowledge of the TS structures of enzymes reacting with non-natural carbasugar substrates shows that GHs can stabilize a wider range of positively charged TS structures than previously thought. Furthermore, this enhanced understanding will enable the design of new carbasugar GH transition state analogues to be used as, for example, chemical biology tools and pharmaceutical lead compounds.

Positive charge stabilized on remote C5-allylic center with catalysis occurring via a loose S_N2 transition state.     

Catalytic resonance theory: parallel reaction pathway control

Catalytic enhancement of chemical reactions via heterogeneous materials occurs through stabilization of transition states at designed active sites, but dramatically greater rate acceleration on that same active site can be achieved when the surface intermediates oscillate in binding energy. The applied oscillation amplitude and frequency can accelerate reactions orders of magnitude above the catalytic rates of static systems, provided the active site dynamics are tuned to the natural frequencies of the surface chemistry. In this work, differences in the characteristics of parallel reactions are exploited via selective application of active site dynamics (0 < ΔU < 1.0 eV amplitude, 10⁻⁶ < f < 10⁴ Hz frequency) to control the extent of competing reactions occurring on the shared catalytic surface. Simulation of multiple parallel reaction systems with broad range of variation in chemical parameters revealed that parallel chemistries are highly tunable in selectivity between either pure product, even when specific products are not selectively produced under static conditions. Two mechanisms leading to dynamic selectivity control were identified: (i) surface thermodynamic control of one product species under strong binding conditions, or (ii) catalytic resonance of the kinetics of one reaction over the other. These dynamic parallel pathway control strategies applied to a host of simulated chemical conditions indicate significant potential for improving the catalytic performance of many important industrial chemical reactions beyond their existing static performance.

Branched catalytic reaction networks with oscillating chemical pathways perfectly select for reaction products at varying frequency.     

State-selective frustration as a key driver of allosteric pluripotency

Allosteric pluripotency arises when an allosteric effector switches from agonist to antagonist depending on the experimental conditions. For example, the Rp-cAMPS ligand of Protein Kinase A (PKA) switches from agonist to antagonist as the MgATP concentration increases and/or the kinase substrate affinity or concentration decreases. Understanding allosteric pluripotency is essential to design effective allosteric therapeutics with minimal side effects. Allosteric pluripotency of PKA arises from divergent allosteric responses of two homologous tandem cAMP-binding domains, resulting in a free energy landscape for the Rp-cAMPS-bound PKA regulatory subunit R1a in which the ground state is kinase inhibition-incompetent and the kinase inhibition-competent state is excited. The magnitude of the free energy difference between the ground non-inhibitory and excited inhibitory states (ΔG_R,Gap) relative to the effective free energy of R1a binding to the catalytic subunit of PKA (ΔG_R:C) dictates whether the antagonism-to-agonism switch occurs. However, the key drivers of ΔG_R,Gap are not fully understood. Here, by analyzing an R1a mutant that selectively silences allosteric pluripotency, we show that a major determinant of ΔG_R,Gap unexpectedly arises from state-selective frustration in the ground inhibition-incompetent state of Rp-cAMPS-bound R1a. Such frustration is caused by steric clashes between the phosphate-binding cassette and the helices preceding the lid, which interact with the phosphate and base of Rp-cAMPS, respectively. These clashes are absent in the excited inhibitory state, thus reducing the ΔG_R,Gap to values comparable to ΔG_R:C, as needed for allosteric pluripotency to occur. The resulting model of allosteric pluripotency is anticipated to assist the design of effective allosteric modulators.

The Rp-cAMPS ligand of protein kinase A switches from agonist to antagonist depending on metabolite and proteomic contexts. We show that the state-selective frustration is a key driver of this allosteric pluripotency phenomenon.     

The Effect of Intramolecular Hydrogen Bond Type on the Gas-Phase Deprotonation of ortho-Substituted Benzenesulfonic Acids. A Density Functional Theory Study

Structural factors have been identified that determine the gas-phase acidity of ortho-substituted benzenesulfonic acid, 2-XC₆H₄–SO₃H, (X = –SO₃H, –COOH, –NO₂, –SO₂F, –C≡N, –NH₂, –CH₃, –OCH₃, –N(CH₃)₂, –OH). The DFT/B3LYP/cc-pVTZ method was used to perform conformational analysis and study the structural features of the molecular and deprotonated forms of these compounds. It has been shown that many of the conformers may contain anintramolecular hydrogen bond (IHB) between the sulfonic group and the substituent, and the sulfonic group can be an IHB donor or an acceptor. The Gibbs energies of gas-phase deprotonation Δ_rG⁰₂₉₈ (kJ mol^–1) were calculated for all compounds. It has been set that in ortho-substituted benzenesulfonic acids, the formation of various types of IHB is possible, having a significant effect on the Δ_rG⁰₂₉₈ values of gas-phase deprotonation. If the –SO₃H group is the IHB donor, then an ion without an IHB is formed upon deprotonation, and the deprotonation energy increases. If this group is an IHB acceptor, then a significant decrease in Δ_rG⁰₂₉₈ of gas-phase deprotonation is observed due to an increase in IHB strength and the A⁻ anion additional stabilization. A proton donor ability comparative characteristic of the –SO₃H group in the studied ortho-substituted benzenesulfonic acids is given, and the Δ_rG⁰₂₉₈ energies are compared with the corresponding values of ortho-substituted benzoic acids.     

Characteristics of transition state and reactivity of hydrocarbon molecules in reaction of hydrogen atom abstraction by peroxide radicals

In the MINDO/3 approximation, transition states (TSs) have been found for the reaction of hydrogen atom abstraction by hydroperoxide radicals from molecules containing a C-H bond in the -position relative to the double bond. It has been shown that the TSs have a structure with partial charge transfer, and the reactivity of the molecules depends on their electron-donor properties and the TS structure.L. M. Litvinenko Institute of Physical Organic Chemistry and Coal-Tar Chemistry, Academy of Sciences of the Ukrainian SSR, Donetsk. Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 27, No. 4, pp. 479-482, July–August, 1991. Original article submitted February 25, 1991.     

Computational Study on the Mechanism of the Photouncaging Reaction of Vemurafenib: Toward an Enhanced Photoprotection Approach for Photosensitive Drugs

The photochemical behavior of the photosensitive first-line anticancer drug vemurafenib (VFB) is of great interest due to the impact of such behavior on its pharmacological activity. In this work, we computationally elucidated the mechanism of the photoinduced release of VFB from the 4,5-dimethoxy-2-nitrobenzene (DMNB) photoprotecting group by employing various density functional theory (DFT)/time-dependent DFT (TD-DFT) approaches. The computational investigations included a comparative assessment of the influence of the position of the photoprotecting group as a substituent on the thermodynamics and kinetics of the photouncaging reactions of two VFB-DMNB prodrugs, namely pyrrole (N_P) and sulfonamide (N_S). With the aid of the DFT calculations concerning the activation energy barrier (∆G^‡), the obtained results suggest that the step of the photoinduced intramolecular proton transfer of the DMNB moiety is not detrimental concerning the overall reaction profile of the photouncaging reaction of both prodrugs. However, the obtained results suggested that the position of the substitution position of the DMNB photoprotecting group within the prodrug structure has a substantial impact on the photouncaging reaction. In particular, the DMNB-N_s-VFB prodrug exhibited a notable increase in ∆G^‡ for the key step of ring opining within the DMNB moiety indicative of potentially hindered kinetics of the photouncaging process compared with DMNB-N_p-VFB. Such an increase in ∆G^‡ may be attributed to the electronic influence of the N_P fragment of the prodrug. The results reported herein elaborate on the mechanism of the photoinduced release of an important anticancer drug from photoprotecting groups with the aim of enhancing our understanding of the photochemical behavior of such photosensitive pharmaceutical materials at the molecular level.     

Unveiling the Ionic Diels–Alder Reactions within the Molecular Electron Density Theory

The ionic Diels–Alder (I-DA) reactions of a series of six iminium cations with cyclopentadiene have been studied within the Molecular Electron Density Theory (MEDT). The superelectrophilic character of iminium cations, ω > 8.20 eV, accounts for the high reactivity of these species participating in I-DA reactions. The activation energies are found to be between 13 and 20 kcal·mol⁻¹ lower in energy than those associated with the corresponding Diels–Alder (DA) reactions of neutral imines. These reactions are low endo selective as a consequence of the cationic character of the TSs, but highly regioselective. Solvents have poor effects on the relative energies, and an unappreciable effect on the geometries. In acetonitrile, the activation energies increase slightly as a consequence of the better solvation of the iminium cations than the cationic TSs. Electron localization function (ELF) topological analysis of the bonding changes along the I-DA reactions shows that they are very similar to those in polar DA reactions. The present MEDT study establishes that the global electron density transfer (GEDT) taking place at the TSs of I-DA reactions, and not steric (Pauli) repulsions such as have been recently proposed, are responsible for the features of these types of DA reactions.     

Computational design of an amidase by combining the best electrostatic features of two promiscuous hydrolases

While there has been emerging interest in designing new enzymes to solve practical challenges, computer-based options to redesign catalytically active proteins are rather limited. Here, a rational QM/MM molecular dynamics strategy based on combining the best electrostatic properties of enzymes with activity in a common reaction is presented. The computational protocol has been applied to the re-design of the protein scaffold of an existing promiscuous esterase from Bacillus subtilis Bs2 to enhance its secondary amidase activity. After the alignment of Bs2 with a non-homologous amidase Candida antarctica lipase B (CALB) within rotation quaternions, a relevant spatial aspartate residue of the latter was transferred to the former as a means to favor the electrostatics of transition state formation, where a clear separation of charges takes place. Deep computational insights, however, revealed a significant conformational change caused by the amino acid replacement, provoking a shift in the pK_a of the inserted aspartate and counteracting the anticipated catalytic effect. This prediction was experimentally confirmed with a 1.3-fold increase in activity. The good agreement between theoretical and experimental results, as well as the linear correlation between the electrostatic properties and the activation energy barriers, suggest that the presented computational-based investigation can transform in an enzyme engineering approach.

A computational strategy, based on combining the best electrostatic properties of enzymes with activity on a common reaction, is presented and applied to the re-design of the protein scaffold of an promiscuous esterase to enhance its secondary amidase activity.

The application of enzymes for desired chemical transformations has been demonstrated by the report of novel and functional designed structures.^1–5 Recent advances in molecular biology and screening technologies have enabled the creation of enzymes via directed evolution. By mimicking the process of natural evolution, iterative cycles of (semi-)random mutations facilitate the improvement of proteins in the laboratory through screening and selection, and hence the identification of active variants.^6–12 Minimal structural information is needed for this strategy and distal sites critical for enzyme catalysis can also be identified. Nevertheless, directed evolution is limited by the fact that, even with the most efficient high-throughput system, only a fraction of all the possible mutants of a given enzyme can be sampled within a set timeframe.¹³ Furthermore, the development of an efficient screening system for a tailored reaction remains challenging. Recently machine-learning (ML) methods have been proposed to expedite evolution and expand the number of properties that can be optimized.^14,15 However, in order to create enzymes with novel reactivities by means of ML methods, protein engineers will have to use proteins with sequences not assigned to the designated reaction or with properties other than those of specific interest, which currently is a technical challenge. Sequence–function data from engineering experiments must be collected to catalogue the natural diversity of proteins in order to convert ML into a useful tool.¹⁵An alternative approach is a rational design, a technique that modifies selected residues at specific positions of an already existing protein scaffold through the analysis of existing mechanistic and structural data.¹⁶ To reveal the structures of the protein in the full catalytic process under physiological conditions, including metastable transition state (TS) structures, computer simulations are essential. Among all the computer-assisted design strategies, two philosophies can be identified: the redesign of the active site of an existing substrate-promiscuous enzyme and the de novo design that constructs an enzyme “from scratch”. The use of promiscuous enzymes is found to be a very promising starting point for the design of new and highly efficient biocatalysts.^17,18 However, the knowledge about the particular molecular mechanisms that allow enzymes to catalyze more than one chemical reaction is still under debate.^19–22Because both the enzyme redesign and the de novo design approach require knowledge of the TS of the reaction to be catalyzed, quantum mechanical (QM) calculations offer crucial complementary information that accelerates the development of novel designed reactions. Moreover, multiscale methods are the only tool that can offer a detailed atomistic picture of the reactions in the active site of the enzyme, which can be dramatically different from that in the gas phase or solution. In multiscale methods, electrons of the reacting fragments are explicitly described by QM methods and the large and complex interacting environment (the fully solvated protein) is described by molecular mechanics (MM) force fields. The mechanism of a reaction in the active site of an enzyme can be determined within these hybrid QM/MM methods through the extensive exploration of the Free Energy Surface (FES). This allows the determination of the rate-limiting step in a multi-step process and, within the framework of Transition State Theory (TST),²³ the prediction of rate constants directly comparable with experiments. Previous studies combining computer simulations with experimental kinetic measurements have demonstrated the good agreement that can be achieved,^24,25 which obviously depends on the quality of both simulations and experiments. In this regard, the error in the determination of activation free energies associated with the use of computational methods such as the umbrella sampling method,^26,27 employed in the present study, is usually accepted to be within 1 kcal mol⁻¹.²⁸Optimizing the secondary activity of promiscuous enzymes is a non-trivial challenge as can be illustrated by analysing the Bacillus subtilis esterase Bs2. While Bs2 is recognized as a serine hydrolase whose primary reaction is the hydrolysis of esters, it can also catalyze the hydrolysis of the amide bond of N-(4-nitrophenyl)-butyramide as a secondary reaction (Fig. 1).^22,29 Previously directed evolution experiments by Arnold and co-workers resulted in a 7-mutation variant with a 100-fold enhancement of the esterase activity (using para-nitrobenzyl butyrate as the substrate).³⁰ Bornscheuer and co-workers used a combination of directed evolution and rational design based on docking and classical energy minimization to get a 3-fold increase of the amidase activity of Bs2 after two single mutations.²² In a larger context, despite the successes of different computer-assisted designs of new enzymes, it has been argued that the high activities of the best artificial enzymes have been largely due to directed evolution and the contribution of computation was comparatively modest.³¹Open in a separate windowFig. 1Schematic representation of the reaction mechanism of the hydrolysis of N-(4-nitrophenyl)-butyramide catalyzed by Bs2. (a) Acylation step: the nucleophilic addition of Ser189 to the carbonyl followed by the breaking of the C–N bond is triggered by His399-assisted proton shuffling, and the leaving group, in this case, is 4-nitroaniline. (b) Hydrolysis step: the nucleophilic addition of a water molecule followed by the resolution of the acyl–enzyme complex is triggered by His399-assisted protein shuffling, yielding butyric acid as a product.We envisaged that creating mutations to optimize the preorganization of the protein environment will result in a variant that exhibits improved activity for the desired reaction.³² Based on our recent QM/MM studies of different enzymatic reactions, we have quantified and shown how the reactivity of different proteins can be rationalized from their electrostatic properties,^{24,25,33–36} as the pioneering studies reported by Warshel and co-workers.^37–39 The computed changes of the electrostatic potential or the electric field exerted by the studied proteins on the key atoms of the substrates reflect that there is a small reorganization of these entities when evolving from the reactant state (RS) to the TS at the lowest energy cost.^24,33–36 The electrostatic effects within the active site of the enzyme, therefore, appear to be critical for the electronic reorganization of the reactants during chemical transformations. These studies support the idea that the electrostatic properties of enzymes are the origin of their catalytic features;⁴⁰ consequently, we view that a detailed understanding of the molecular mechanism, including the evolution of electrostatic potential generated in the active site of the enzyme, could be useful in future computer-assisted protein design methods.To engineer enzymes with optimal electrostatic preorganization, comparative analysis between unrelated natural enzymes that catalyze the same chemical reactions can be a reliable strategy. In previous studies, we have shown that Candida antarctica lipase B (CALB) also displays amidase activity similar to that of Bs2, though being non-homologous with each other.^24,25 QM/MM studies of the amidase reaction catalyzed by wild-type Bs2 and CALB enzymes were previously conducted.^24,25 We expect that the favorable features of each enzyme could be isolated and combined to create a redesigned enzyme with improved catalytic activity for the secondary amidase reaction. In the present paper, based on our knowledge derived from previous comparative studies, and by applying the concept of electrostatic pre-organization,^{24,33–36,40–43} a variant with improved activity for the designated amidase reaction was generated. After overlapping the structures of both proteins in one of the located TSs, through the use of a rotation quaternion around selected atoms of the substrate, a catalytically improved Bs2 variant was delineated. In particular, residues of Bs2 with an unfavorable electrostatic effect on catalysis were substituted by those placed in an equivalent spatial position in CALB with a favorable effect, as explained in detail below. The QM/MM FES of the full catalytic reaction in the proposed variant, combined with the experimental characterization, will be used to propose a general computer-based strategy that can be potentially used to design new enzymes.     

Explorations on Thermodynamic and Kinetic Performances of Various Cationic Exchange Durations for Synthetic Clinoptilolite

Various cation–exchanged clinoptilolites (M–CPs, M = Li⁺, Cs⁺, Ca²⁺, Sr²⁺) were prepared, and their exchanged thermodynamic (and kinetic) properties and adsorption performances for CH₄, N₂, and CO₂ were investigated. The results demonstrated that the relative crystallinity of M–CP_S decreased with the increase of exchange times. Their chemisorbed water weight loss gradually increased with the increasing exchange times, except that of Cs–x–CP. The Δ_rGmθ values of exchange process of Li⁺, Cs⁺, Ca²⁺, or Sr² presented the increased trend with the enhanced exchange times, but they decreased as the temperature increased. The negative Δ_rGmθ values and the positive Δ_rHmθ and Δ_rSmθ values suggested that the exchanged procedure belonged to spontaneous, endothermic, and entropy-increasing behaviors; their kinetic performances followed a pseudo–second–order model. However, the calculated E_a values of exchange process showed the increased tendencies with the enhanced exchange times, indicating that the exchange process became more difficult. Finally, the preliminary adsorption results indicated that the maximum adsorption amount at 273 K and 1 bar was 0.51 mmol/g of CH₄ and 0.38 mmol/g of N₂ by (Na, K)–CP, and 2.32 mmol/g of CO₂ by Li–6–CP.     

Influence of Modified Terminal Segments on Dynamic Modulus and Viscosity of Dendrimer

Theoretical studies of the influence of modified terminal segments (TSs) on the relaxation spectrum of a dendrimer and dendrimer mechanical properties such as dynamic viscosity, η(ω), the elastic, G′(ω), and loss, G″(ω), moduli have been carried out by the Rouse model. Two major types of modified TS have been studied: (i) TS with an attached rigid massive group (i.e., TSs with additional friction) and (ii) TSs with a length different from the length of an inner segment. In the low‐frequency region, G′(ω), G″(ω), and η(ω) increase with the rise of friction of TS. In the high‐frequency region, dynamic moduli and viscosity depend on the length of TS. In the intermediate region, the moduli and viscosity are determined by a combined parameter: the characteristic time of TS, τ_end, which depends on the friction and length of TS. For both types of TSs, the position of the G″(ω) maximum, ω_max, depends on τ_end. In most of the considered cases, the linear dependence of ω_max on τ_end has been found. The method, which takes into account a deceleration of TS mobility with the rise in the number of generations, n, has been proposed. It was supposed that the effect of the deceleration corresponds to the forming of a dense surface shell with the rise of n, but similar behavior can also be caused by other reasons. In this case, ω_max shifts to the low‐frequency region with an increase in the number of generations. The conclusions of the theory developed in this paper are in agreement with results of the experiment, in which G′(ω) and G″(ω) were obtained for polyamidoamine dendrimers.

     

Identifying the true origins of selectivity in chiral phosphoric acid catalyzed N-acyl-azetidine desymmetrizations

The first catalytic intermolecular desymmetrization of azetidines was reported by Sun and coworkers in 2015 using a BINOL-derived phosphoric acid catalyst (J. Am. Chem. Soc. 2015, 137, 5895–5898). To uncover the mechanism of the reaction and the origins of the high enantioselectivity, Density Functional Theory (DFT) calculations were performed at the B97D3/6-311+G(2d,2p)/SMD(toluene)//B97D3/6-31G(d,p)/CPCM(toluene) level of theory. Comparison of four possible activation modes confirms that this reaction proceeds through the bifunctional activation of the azetidine nitrogen and the thione tautomer of the 2-mercaptobenzothiazole nucleophile. Upon thorough conformational sampling of the enantiodetermining transition structures (TSs), a free energy difference of 2.0 kcal mol⁻¹ is obtained, accurately reproducing the experimentally measured 88% e.e. at 80 °C. This energy difference is due to both decreased distortion and increased non-covalent interactions in the pro-(S) TS. To uncover the true origins of selectivity, the TSs optimized with the full catalyst were compared to those optimized with a model catalyst through steric maps. It is found that the arrangements displayed by the substrates are controlled by strict primary orbital interaction requirements at the transition complex, and their ability to fit into the catalyst pocket drives the selectivity. A general model of selectivity for phosphoric acid-catalyzed azetidine desymmetrizations is proposed, which is based on the preference of the nucleophile and benzoyl group to occupy empty quadrants of the chiral catalyst pocket.

The origins of selectivity in azetidine desymmetrizations have been determined computationally. Comparison of structures with model and full catalysts provided key details missed by typical analyses of the stereodetermining transition structures.     

Transition structure selectivity in enzyme catalysis: a QM/MM study of chorismate mutase

Two different transition structures (TSs) have been located and characterized for the chorismate conversion to prephenate in Bacillus subtilis chorismate mutase by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) calculations. GRACE software, combined with an AM1/CHARMM24/TIP3P potential, has been used involving full gradient relaxation of the position of ca. 3300 atoms. These TSs have been connected with their respective reactants and products by the intrinsic reaction coordinate (IRC) procedure carried out in the presence of the protein environment, thus obtaining for the first time a realistic enzymatic reaction path for this reaction. Similar QM/MM computational schemes have been applied to study the chemical reaction solvated by ca. 500 water molecules. Comparison of these results together with gas phase calculations has allowed understanding of the catalytic efficiency of the protein. The enzyme stabilizes one of the TSs (TS_OHout) by means of specific hydrogen bond interactions, while the other TS (TS_OHin) is the preferred one in vacuum and in water. The enzyme TS is effectively more polarized but less dissociative than the corresponding solvent and gas phase TSs. Electrostatic stabilization and an intramolecular charge-transfer process can explain this enzymatically induced change. Our theoretical results provide new information on an important enzymatic transformation and the key factors responsible for efficient selectivity are clarified. Received: 25 March 2000 / Accepted: 7 August 2000 / Published online: 23 November 2000