Acyltransferases isolated from Pseudomonas protegens (PpATase) and Pseudomonas fluorescens (PfATase) have recently been reported to catalyze the Friedel-Crafts acylation, providing a biological version of this classical organic reaction. These enzymes catalyze the cofactor-independent acylation of monoacetylphloroglucinol (MAPG) to diacetylphloroglucinol (DAPG) and phloroglucinol (PG) and have been demonstrated to have a wide substrate scope, making them valuable for potential applications in biocatalysis. Herein, we present a detailed reaction mechanism of PpATase on the basis of quantum chemical calculations, employing a large model of the active site. The proposed mechanism is consistent with available kinetics, mutagenesis, and structural data. The roles of various active site residues are analyzed. Very importantly, the Asp137 residue, located more than 10 angstrom from the substrate, is predicted to be the proton source for the protonation of the substrate in the rate-determining step. This key prediction is corroborated by site-directed mutagenesis experiments. Based on the current calculations, the regioselectivity of PpATase and its specificity toward non-natural substrates can be rationalized.

The computational study of enantioselective reactions is a challenging task that requires high accuracy, as very small energy differences have to be reproduced. Quantum chemical methods, most commonly density functional theory, are today an important tool in this pursuit. This Perspective describes recent efforts in modeling asymmetric reactions in enzymes by means of the quantum chemical cluster approach. The methodology is described briefly and a number of illustrative case studies performed recently at our laboratory are presented. The reviewed enzymes are limonene epoxide hydrolase, soluble epoxide hydrolase, arylmalonate decarboxylase, phenolic acid decarboxylase, benzoylformate decarboxylase, secondary alcohol dehydrogenase, acyl transferase, and norcoclaurine synthase. The challenges encountered in each example are discussed, and the modeling lessons learned are highlighted.

The acyltransferase from Mycobacterium smegmatis (MsAcT) complements the well-established acylation activity of hydrolases in organic solvents with its activity to perform acylation reactions (among other reactions) in an aqueous environment. The enzyme’s potential is however limited, due to its poor regio- and stereoselectivity with enantioselectivities (E-values) below 20 for bulky (aromatic) substrates. By applying computer-guided rational design, a library of single variants was designed that allowed conversion of a set of previously challenging substrates with good activity and E-values up to >200. The computational predictions were found to be in agreement with experimental data, which in turn allowed for the generation of even more active and selective double variants. Overall, the produced set of variants provides a toolbox for the enantioselective acylation of challenging alcohols in water, effectively contributing to an alternative to reactions in organic solvents.

Acyl transferase from Mycobacterium smegmatis (MsAcT) is a promising biocatalyst because it catalyzes an acyl transfer reaction in aqueous solution, thereby accepting many primary and secondary alcohols as substrates. MsAcT also exhibits high enantioselectivity for a selected number of secondary alcohols. To increase the applicability of this enzyme for the production of optically active compounds, a detailed understanding of the reaction mechanism and the factors that affect enantioselectivity is essential. Herein, quantum chemical calculations are employed to study the reactions of two secondary alcohols, 1-isopropyl propargyl alcohol and 2-hydroxy propanenitrile, for which the enzyme displays opposite enantiopreference, favoring the S enantiomer in the former case and R enantiomer in the latter. A model of the active site has been designed and for both substrates various binding modes are evaluated and the intermediates and transition states along the reaction path are then located. The calculated energy profiles agree with the experimental observations, and reproduce the selectivity outcome. Through a detailed analysis of the geometries of key transition states, insights into the origins of the enantiopreference are obtained.

The acyl transferase from Mycobacterium smegmatis (MsAcT) catalyzes the acyl transfer between a range of primary and secondary alcohols, whereby its outstanding ability is to perform this reaction in aqueous solution. Therefore, MsAcT opens different options for acylation reactions enabling alternatives for many conventionally hydrolytic enzymes used in biocatalysis. Nevertheless, hydrolysis is still a major side reaction of this enzyme. To provide a detailed understanding of the competition between hydrolysis and transesterification reactions, a combination of density functional theory and free energy perturbation methods have been employed. The relative binding free energies and the energy profiles of the chemical steps involved in the reaction were calculated for a number of substrates. The calculations show that the enzyme active site exhibits a higher affinity for substrates with an aromatic ring. The rate-determining step corresponds to the collapse of a negatively charged tetrahedral intermediate in the substrate acylation half-reaction. The intrinsic barriers of the transesterification and hydrolysis half-reactions are calculated to be of similar heights, suggesting that the determining factor in the MsAcT specificity is the higher binding affinity of the active site for the alcohol substrates relative to water. Finally, the influence of the acyl donor on the MsAcT-catalyzed reaction is also investigated by considering different esters in the calculations.

The peptidyl transfer reaction on the large ribosomal subunit depends on the protonation state of the amine nucleophile and exhibits a large kinetic solvent isotope effect (KSIE similar to 8). In contrast, the related peptidyl-tRNA hydrolysis reaction involved in termination shows a KSIE of similar to 4 and a pH-rate profile indicative of base catalysis. It is, however, unclear why these reactions should proceed with different mechanisms, as the experimental data suggests. One explanation is that two competing mechanisms may be operational in the peptidyl transferase center (PTC). Herein, we explored this possibility by re-examining the previously proposed proton shuttle mechanism and testing the feasibility of general base catalysis also for peptide bond formation. We employed a large cluster model of the active site and different reaction mechanisms were evaluated by density functional theory calculations. In these calculations, the proton shuttle and general base mechanisms both yield activation energies comparable to the experimental values. However, only the proton shuttle mechanism is found to be consistent with the experimentally observed pH-rate profile and the KSIE. This suggests that the PTC promotes the proton shuttle mechanism for peptide bond formation, while prohibiting general base catalysis, although the detailed mechanism by which general base catalysis is excluded remains unclear.