Selectivity of lipases: Conformational analysis of suggested intermediates in ester hydrolysis of chiral primary and secondary alcohols

Directed evolution of an enantioselective epoxide hydrolase: uncovering the source of enantioselectivity at each evolutionary stage

Directed evolution of enzymes as enantioselective catalysts in organic chemistry is an alternative to traditional asymmetric catalysis using chiral transition-metal complexes or organocatalysts, the different approaches often being complementary. Moreover, directed evolution studies allow us to learn more about how enzymes perform mechanistically. The present study concerns a previously evolved highly enantioselective mutant of the epoxide hydrolase from Aspergillus niger in the hydrolytic kinetic resolution of racemic glycidyl phenyl ether. Kinetic data, molecular dynamics calculations, molecular modeling, inhibition experiments, and X-ray structural work for the wild-type (WT) enzyme and the best mutant reveal the basis of the large increase in enantioselectivity (E = 4.6 versus E = 115). The overall structures of the WT and the mutant are essentially identical, but dramatic differences are observed in the active site as revealed by the X-ray structures. All of the experimental and computational results support a model in which productive positioning of the preferred (S)-glycidyl phenyl ether, but not the (R)-enantiomer, forms the basis of enhanced enantioselectivity. Predictions regarding substrate scope and enantioselectivity of the best mutant are shown to be possible.

Combined X-ray diffraction and QM/MM study of the Burkholderia cepacia lipase-catalyzed secondary alcohol esterification

To understand the origin of high enantioselectivity of Burkholderia cepacia lipase (BCL) toward secondary alcohol, (R,S)-1-phenoxy-2-hydroxybutane (1), and its ester (E1), we determined the crystal structure of BCL complexed with phosphonate analogue of S-E1 and accomplished a series of MM, MC, and QM/MM studies. We have found that the inhibitor in the S configuration binds into the BCL active site in the same manner as the R isomer, with an important difference: while in case of the R-inhibitor the H-bond between its alcohol oxygen and catalytic His286 can be formed, in the case of the S-inhibitor this is not possible. Molecular modeling for both E1 enantiomers revealed orientations in which all hydrogen bonds characteristic of productive binding are formed. To check the possibility of chemical transformation, four different orientations of the substrate (two for each enantiomer) were chosen, and a series of ab initio QM/MM calculations were accomplished. Starting from the covalent complex, we modeled the ester (E1) hydrolysis and the alcohol (1) esterification. The calculations revealed that ester release is possible starting with all four covalent complexes. Alcohol release from the BCL-E1 complex in which the S-substrate is bound in the same manner as the S-inhibitor in the crystal structure however is not possible. These results show that the crystallographically determined binding modes should be taken with caution when modeling chemical reactions.

Enthalpie and Entropie Contributions in the Transesterification of Sucrose: Computational Study of Lipases and Subtilisin

Transesterification of sucrose with fatty acids catalyzed by subtilisin Carlsberg occurs with regioselectivity that is different from that in lipases. Thermomyces lanuginosus lipase (TlL) and Candida antarctica lipase B (CALB) catalyze synthesis at positions 6 and 6', with differing abilities, while subtilisin catalysis leads to the 1'-acylated sucrose. The catalytic machinery in lipases is approximately mirrored in subtilisins but different pocket morphologies including size, shape, and rearrangement of the catalytic elements underlies the differing regioselectivities. The thermodynamic consequences of these differences on the above reactions have been explored systematically using computational methods, determining the free energies of interaction of the putative transition-state adducts. Analysis of the conformers with the lowest transition state energies (protein-ligand interactions and vibrational entropy contributions) indicates that enthalpic factors control specificities in lipases while entropic factors are more important in subtilisin.

Expanding the substrate scope of enzymes: combining mutations obtained by CASTing

In a previous paper, the combinatorial active-site saturation test (CAST) was introduced as an effective strategy for the directed evolution of enzymes toward broader substrate acceptance. CASTing comprises the systematic design and screening of focused libraries around the complete binding pocket, but it is only the first step of an evolutionary process because only the initial libraries of mutants are considered. In the present study, a simple method is presented for further optimization of initial hits by combining the mutational changes obtained from two different libraries. Combined lipase mutants were screened for hydrolytic activity against six notoriously difficult substrates (bulky carboxylic acid esters) and improved mutants showing significantly higher activity were identified. The enantioselectivity of the mutants in the hydrolytic kinetic resolution of two substrates was also studied, with the best mutant-substrate combination resulting in a selectivity factor of E=49. Finally, the catalytic profile of the evolved mutants in the hydrolysis of simple nonbranched carboxylic acid esters, ranging from acetate to palmitate, was studied for theoretical reasons.

Combining regio- and enantioselectivity of lipases for the preparation of (R)-4-chloro-2-butanol

Preparation of 98% ee (R)-4-chloro-2-butanol was carried out by the enzymatic hydrolysis of chlorohydrin esters, using fungal resting cells and commercial enzymes. Hydrolyzes were carried out using lipases from Candida antarctica (Novozym 435), C. rugosa, Rhizomucor miehei (Lipozyme IM), Burkolia cepacia, and resting cells of Rhizopus oryzae and Aspergillus flavus. The influence of the enzyme, the solvent, the temperature, and the alkyl chain length on the selectivity of hydrolyzes of isomeric mixtures of chlorohydrin esters is described. Regioselectivity was higher than 95% for some of the tested lipases. Novozym 435 allowed preparation of the (R)-4-chloro-2-butanol after 15 min of reaction at 30-40 degrees C.

An Inverse Substrate Orientation for the Regioselective Acylation of 3′,5′-Diaminonucleosides Catalyzed by Candida antarctica lipase B?

Candida antarctica lipase B (CAL-B) catalyzes the regioselective acylation of natural thymidine with oxime esters and also the regioselective acylation of an analogue, 3',5'-diamino-3',5'-dideoxythymidine with nonactivated esters. In both cases, acylation favors the less hindered 5'-position over the 3'-position by upto 80-fold. Computer modeling of phosphonate transition-state analogues for the acylation of thymidine suggests that CAL-B favors acylation of the 5'-position because this orientation allows the thymine ring to bind in a hydrophobic pocket and forms stronger key hydrogen bonds than acylation of the 3'-position. On the other hand, computer modeling of phosphonamidate analogues of the transition states for acylation of either the 3'- or 5'-amino groups in 3',5'-diamino-3',5'-dideoxythymidine shows similar orientations and hydrogen bonds and, thus, does not explain the high regioselectivity. However, computer modeling of inverse structures, in which the acyl chain binds in the nucleophile pocket and vice versa, does rationalize the observed regioselectivity. The inverse structures fit the 5'-, but not the 3'-intermediate thymine ring, into the hydrophobic pocket, and form a weak new hydrogen bond between the O-2 carbonyl atom of the thymine and the nucleophile amine only for the 5'-intermediate. A water molecule might transfer a proton from the ammonium group to the active-site histidine. As a test of this inverse orientation, we compared the acylation of thymidine and 3',5'-diamino-3',5'-dideoxythymidine with butyryl acyl donors and with isosteric methoxyacetyl acyl donors. Both acyl donors reacted at equal rates with thymidine, but the methoxyacetyl acyl donor reacted four times faster than the butyryl acyl donor with 3',5'-diamino-3',5'-dideoxythymidine. This faster rate is consistent with an inverse orientation for 3',5'-diamino-3',5'-dideoxythymidine, in which the ether oxygen atom of the methoxyacetyl group can form a similar hydrogen bond to the nucleophilic amine. This combination of modeling and experiments suggests that such lipase-catalyzed reactions of apparently close substrate analogues like alcohols and amines might follow different pathways.

Mirror-Image Packing in Enantiomer Discrimination

Synthetic chemists often exploit the high enantioselectivity of lipases to prepare pure enantiomers of primary alcohols, but the molecular basis for this enantioselectivity is unknown. The crystal structures of two phosphonate transition-state analogs bound to Burkholderia cepacia lipase reveal this molecular basis for a typical primary alcohol: 2-methyl-3-phenyl-1-propanol. The enantiomeric alcohol moieties adopt surprisingly similar orientations, with only subtle differences that make it difficult to predict how to alter enantioselectivity. These structures, along with a survey of previous structures of enzyme bound enantiomers, reveal that binding of enantiomers does not involve an exchange of two substituent positions as most researchers assumed. Instead, the enantiomers adopt mirror-image packing, where three of the four substituents at the stereocenter lie in similar positions. The fourth substituent, hydrogen, points in opposite directions.

Learning from Directed Evolution: Theoretical Investigations into Cooperative Mutations in Lipase Enantioselectivity

Molecular modeling with classical force-fields has been used to study the reactant complex and the tetrahedral intermediate in lipase-catalyzed ester hydrolysis in 20 enzyme/substrate combinations. The R and S enantiomers of alpha-methyldecanoic acid ester served as substrates for the wild-type lipase from Pseudomonas aeruginosa and nine selected mutants. After suitable preparation of initial structures from an available wild-type crystal structure, each system was subjected to 1 ns CHARMM force-field molecular dynamics simulations. The resulting geometric and energetic changes allow interpretation of some experimentally observed effects of mutations, particularly with regard to the "hot spots" at residues 155 and 162. The replacement S155F enhances S enantiopreference through a steric relay involving Leu162. The double mutation S53P + L162G improves S enantioselectivity by creating a new binding pocket for the S enantiomer with an additional stabilizing hydrogen bond to His83. The simulations provide insight into remote and cooperative effects of mutations.

A quantitative model for predicting enzyme enantioselectivity: application to Burkholderia cepacia lipase and 3-(aryloxy)-1,2-propanediol derivatives

We describe a new approach for predicting the enantioselectivity of enzymes towards racemic compounds. It is based on comparative binding energy (COMBINE) analysis. The approach is used to rationalise the enantioselectivity of Burkholderia cepacia lipase (BCL) towards thirteen racemic 3-(aryloxy)-1,2-propanediols in the process of acylation. According to our molecular modelling study the two 3-(aryloxy)-1,2-propanediols enantiomers bind in the BCL active site in different orientations. To derive a quantitative structure-activity relationship (QSAR), the difference in the interaction energy between two enantiomers with each amino acid residue was computed. These residue-based energy differences were then subjected to chemometric analysis and 3D QSAR models were derived. The models were able to unambiguously predict the fast-reacting enantiomer and the approximate magnitude of the enantioselectivity. The study enabled identification of interactions between the substrate and the lipase amino acid residues that play key roles in secondary alcohol enantiodifferentiation. From the results, it was possible to propose modifications of both, substrate and protein, which would directionally modify enantioselectivity of BCL towards secondary aryl-alcohols.

Complex of Burkholderia cepacia lipase with transition state analogue of 1-phenoxy-2-acetoxybutane: biocatalytic, structural and modelling study

In a series of four racemic phenoxyalkyl-alkyl carbinols, 1-phenoxy-2-hydroxybutane (1) is enantioselectively acetylated by Burkholderia cepacia (formerly Pseudomonas cepacia) lipase with an E value > or = 200, whereas for the other three racemates E was found to be < or = 4. To explain the high preference of B. cepacia lipase for (R)-(+)-1, a precursor of its transition state analogue with a tetrahedral P-atom, (R(P),S(P))-O-(2R)-(1-phenoxybut-2-yl)methylphosphonic acid chloride was prepared and crystallized in complex with B. cepacia lipase. The X-ray structure of the complex was determined, allowing to compare the conformation of the inhibitor with results of molecular modelling.

Novel microbial lipases: catalytic activity in reactions in organic media

2,000 microbial strains were isolated from soil samples and tested to determine their lipolytic activity by employing screening techniques on solid and in liquid media. Culture broths were initially tested with 1,2-O-dilauryl-rac-glycero-3-glutaric acid-resorufinyl ester, a chromogenic substrate specific for lipases. Fourteen lipase-producing microorganisms were selected and their taxonomic identification was carried out. Hydrolysis of tributyrin or olive oil and the esterification of oleic acid with heptanol were selected to preliminary evaluate the catalytic activity of these lipases. All the selected lipases catalysed this esterification reaction with good yields. Resolution of (R,S)-2-(4-isobutylphenyl) propionic acid, (R,S)-1-phenylethanol, (R,S) 1-phenylethylamine and of (R) or (S) glycidol were performed to evaluate the stereoselectivity of these novel enzymes as biocatalysts in reactions in organic media. Lipases from the fungi Fusarium oxysporum and Ovadendron sulphureo-ochraceum gave the best yields and enantioselectivities in the resolution of racemic ibuprofen and 1-phenylethanol. Several lipases displayed a high stereoselectivity in the resolution of chiral amines by an alcoxycarbonylation reaction.

Stereoselectivity of Pseudomonas cepacia lipase toward secondary alcohols: a quantitative model

The lipase from Pseudomonas cepacia represents a widely applied catalyst for highly enantioselective resolution of chiral secondary alcohols. While its stereopreference is determined predominantly by the substrate structure, stereoselectivity depends on atomic details of interactions between substrate and lipase. Thirty secondary alcohols with published E values using P. cepacia lipase in hydrolysis or esterification reactions were selected, and models of their octanoic acid esters were docked to the open conformation of P. cepacia lipase. The two enantiomers of 27 substrates bound preferentially in either of two binding modes: the fast-reacting enantiomer in a productive mode and the slow-reacting enantiomer in a nonproductive mode. Nonproductive mode of fast-reacting enantiomers was prohibited by repulsive interactions. For the slow-reacting enantiomers in the productive binding mode, the substrate pushes the active site histidine away from its proper orientation, and the distance d(H(N epsilon) - O(alc)) between the histidine side chain and the alcohol oxygen increases, d(H(N epsilon) - O(alc)) was correlated to experimentally observed enantioselectivity: in substrates for which P. cepacia lipase has high enantioselectivity (E > 100), d(H(N epsilon) - O(alc)) is >2.2 A for slow-reacting enantiomers, thus preventing efficient catalysis of this enantiomer. In substrates of low enantioselectivity (E < 20), the distance d(H(N epsilon) - O(alc)) is less than 2.0 A, and slow- and fast-reacting enantiomers are catalyzed at similar rates. For substrates of medium enantioselectivity (20 < E < 100), d(H(N epsilon) - O(alc)) is around 2.1 A. This simple model can be applied to predict enantioselectivity of P. cepacia lipase toward a broad range of secondary alcohols.

Molecular modeling and biocatalysis: explanations, predictions, limitations, and opportunities

Rapid advances in structural biology have revealed the three-dimensional structures of many biocatalysts. Molecular modeling is the tool that links these structures with experimental observations. As a qualitative tool, current modeling methods are extremely useful. They can explain, on a molecular level, unusual features of reactions. They can predict how to increase the selectivity either by substrate modification or by site-directed mutagenesis. Quantitative predictions, for example the degree of enantioselectivity, are still not reliable, however. Modeling is limited also by the availability of three-dimensional structures. Most current modeling involves hydrolases, especially proteases and lipases, but structures for other types of enzymes are starting to appear.

Influence of water activity on the enantioselective esterification of (R,S)-ibuprofen by Candida antarctica lipase B in solventless media

The lipase-catalyzed enantioselective esterification of ibuprofen has been studied in a media, composed only of substrates. When racemic ibuprofen is used, the alcohol-chain length affects the esterification rates of individual enantiomers, but it does not affect the enantioselectivity. Water activity affects the esterification rates of (R)- and (S)-ibuprofen differently, leading to higher enantioselectivity at lower water activities. Experiments were also conducted at various (R)- to (S)-ibuprofen ratios. It appears that the esterification rate of (R)-ibuprofen is always proportional to its concentration, whereas at low water activity the esterification rate of (S)-ibuprofen shows a saturation at higher concentrations. Other 2-phenyl carboxylic acids were studied, and the increase in apparent enantioselectivity at low-water activity was not observed for the molecules tested.