Search

With advances in protein structure predictions, the number of available high-quality structures has increased dramatically. In light of these advances, structure-based enzyme engineering is expected to become increasingly important for optimizing biocatalysts for industrial processes. Here, we present AsiteDesign, a Monte Carlo-based protocol for structure-based engineering of active sites. AsiteDesign provides a framework for introducing new catalytic residues in a given binding pocket to either create a new catalytic activity or alter the existing one. AsiteDesign is implemented using pyRosetta and incorporates enhanced sampling techniques to efficiently explore the search space. The protocol was tested by designing an alternative catalytic triad in the active site of Pseudomonas fluorescens esterase (PFE). The designed variant was experimentally verified to be active, demonstrating that AsiteDesign can find alternative catalytic triads. Additionally, the AsiteDesign protocol was employed to enhance the hydrolysis of a bulky chiral substrate (1-phenyl-2-pentyl acetate) by PFE. The experimental verification of the designed variants demonstrated that F158L/F198A and F125A/F158L mutations increased the hydrolysis of 1-phenyl-2-pentyl acetate from 8.9 to 66.7 and 23.4%, respectively, and reversed the enantioselectivity of the enzyme from (R) to (S)-enantiopreference, with 32 and 55% enantiomeric excess (ee), respectively., Peer Reviewed, Postprint (published version)

Biotechnological solutions will be a key aspect in our immediate future society, where optimized enzymatic processes through enzyme engineering might be an important solution for waste transformation, clean energy production, biodegradable materials, and green chemistry, for example. Here we advocate the importance of structural-based bioinformatics and molecular modeling tools in such developments. We summarize our recent experiences indicating a great prediction/success ratio, and we suggest that an early in silico phase should be performed in enzyme engineering studies. Moreover, we demonstrate the potential of a new technique combining Rosetta and PELE, which could provide a faster and more automated procedure, an essential aspect for a broader use., This work has also been supported by predoctoral fellowships FPU19/00608 and PRE2020-091825, and the PID2019-106370RBI00/AEI/10.13039/501100011033 grant from the Spanish Ministry of Science and Innovation., Peer Reviewed, Postprint (author's final draft)

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.

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.

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.

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 class I pancreatic elastase from Atlantic salmon is considered to be a cold-adapted enzyme in view of the cold habitat, the reduced thermostability of the enzyme, and the fact that it is faster than its mesophilic porcine counterpart at room temperature. However, no experimental characterization of its catalytic properties at lower temperatures has actually been reported. Here we use extensive computer simulations of its catalytic reaction, at different temperatures and with different peptide substrates, to compare its characteristics with those of porcine pancreatic elastase, with which it shares 67% sequence identity. We find that both enzymes have a preference for smaller aliphatic residues at the P1 position, while the reaction rate with phenylalanine at P1 is predicted to be substantially lower. With the former class of substrates, the calculated reaction rates for salmon enzyme are consistently higher than those of the porcine ortholog at all temperatures examined, and the difference is most pronounced at the lowest temperature. As observed for other cold-adapted enzymes, this is caused by redistribution of the activation free energy in terms of enthalpy and entropy and can be linked to differences in the mobility of surface-exposed loops in the two enzymes. Such mobility changes are found to be reflected by characteristic sequence conservation patterns in psychrophilic and mesophilic species. Hence, calculations of mutations in a single surface loop show that the temperature dependence of the catalytic reaction is altered in a predictable way.

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.

CONSPECTUS: The role played by entropy for the enormous rate enhancement achieved by enzymes has been debated for many decades. There are, for example, several confirmed cases where the activation free energy is reduced by around 10 kcal/mol due to entropic effects, corresponding to a rate enhancement of similar to 10(7) compared to the uncatalyzed reaction. However, despite substantial efforts from both the experimental and theoretical side, no real consensus has been reached regarding the origin of such large entropic contributions to enzyme catalysis. Another remarkable instance of entropic effects is found in enzymes that are adapted by evolution to work at low temperatures, near the freezing point of water. These cold-adapted enzymes invariably show a more negative entropy and a lower enthalpy of activation than their mesophilic orthologs, which counteracts the exponential damping of reaction rates at lower temperature. The structural origin of this universal phenomenon has, however, remained elusive. The basic problem with connecting macroscopic thermodynamic quantities, such as activation entropy and enthalpy derived from Arrhenius plots, to the 3D protein structure is that the underlying detailed (microscopic) energetics is essentially inaccessible to experiment. Moreover, attempts to calculate entropy contributions by computer simulations have mostly focused only on substrate entropies, which do not provide the full picture. We have recently devised a new approach for accessing thermodynamic activation parameters of both enzyme and solution reactions from computer simulations, which turns out to be very successful. This method is analogous to the experimental Arrhenius plots and directly evaluates the temperature dependence of calculated reaction free energy profiles. Hence, by extensive molecular dynamics simulations and calculations of up to thousands of independent free energy profiles, we are able to extract activation parameters with sufficient precisi

Decomposition of activation free energies of chemical reactions, into enthalpic and entropic components, can provide invaluable signatures of mechanistic pathways both in solution and in enzymes. Owing to the large number of degrees of freedom involved in such wcondensed-phase reactions, the extensive configurational sampling needed for reliable entropy estimates is still beyond the scope of quantum chemical calculations. Here we show, for the hydrolytic deamination of cytidine and dihydrocytidine in water, how direct computer simulations of the temperature dependence of free energy profiles can be used to extract very accurate thermodynamic activation parameters. The simulations are based on empirical valence bond models, and we demonstrate that the energetics obtained is insensitive to whether these are calibrated by quantum mechanical calculations or experimental data. The thermodynamic activation parameters are in remarkable agreement with experiment results and allow discrimination among alternative mechanisms, as well as rationalization of their different activation enthalpies and entropies.

Nonsteroidal anti-inflammatory drugs (NSAIDs) are widely used for the treatment of pain, fever, inflammation, and some types of cancers. Their mechanism of action is the inhibition of isoforms 1 and 2 of the enzyme cyclooxygenase (COX-1 and COX-2, respectively). However, both nonselective and selective NSAIDs may have side effects that include gastric intestinal bleeding, peptic ulcer formation, kidney problems, and occurrences of myocardial infarction. The search for selective high-affinity COX inhibitors resulted in a number of compounds characterized by a slow, tight-binding inhibition that occurs in a two-step manner. It has been suggested that the final, only very slowly reversible, tight-binding event is the result of conformational changes in the enzyme. However, the nature of these conformational changes has remained elusive. Here we explore the structural determinants of the tight-binding phenomenon in COX-1 with molecular dynamics and free energy simulations. The calculations reveal how different classes of inhibitors affect the equilibrium between two conformational substates of the enzyme in distinctly different ways. The class of tight-binding inhibitors is found to exclusively stabilize an otherwise unfavorable enzyme conformation and bind significantly stronger to this state than to that normally observed in crystal structures. By also computing free energies of binding to the two enzyme conformations for 16 different NSAIDs, we identify an induced-fit mechanism and the key structural features associated with high-affinity tight binding. These results may facilitate the rational development of new COX inhibitors with improved selectivity profiles.

The aim of this study is pharmacochemistry of Iranian flora Artemisia sieberi and its antimalarial effects on Plasmodium berghei in vivo. This is the first application of A. sieberi for treatment of murine malaria. A. sieberi were collected at flowering stage from the Khorassan and Semnan provinces of Iran; the aerial parts were air-dried at room temperature and then powdered. The powder was macerated in methanol, filtered with Bokhner hopper and solvent was separated in rotary evaporator. Total herbal extract was subsequently processed for ether and chloroform extracts preparation. The toxicity of herbal extract was assessed on naive NMRI mice with high, average and low doses; then pathophysiological signs were assessed. Finally, the antimalarial efficacy was investigated on two groups of Plasmodium berghei infected mice. Percentage of parasitaemia and pathophysiology were also evaluated. The results of this assessment showed no toxicity even by high concentration of herbal extract. A significant reduction in percentage of parasitaemia was observed; no alterations of hepatosplenomegaly and body weight were indicated in study group. A. sieberi extracts showed antimalarial effects against murine malaria with some efficacies on reducing pathophysiology. However, there is requirement to find the major component of this herbal extract by further studies.