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

Bacterial triacylglycerol lipase is a potential cholesterol esterase: Identification of a key determinant for sterol-binding specificity

Cholesterol esterase (Che) from Burkholderia stabilis (BsChe) is a homolog of well-characterized and industrially relevant bacterial triacylglycerol lipases (Lips). BsChe is a rare bacterial Lip enzyme that exhibits practical Che activity and is currently used in clinical applications to determine total serum cholesterol levels. To investigate the sterol specificity of BsChe, we determined the X-ray structure of BsChe. We discovered a local structural change in the active-site cleft, which might be related to substrate binding and product release. We also performed molecular docking studies by using the X-ray models of BsChe and cholesterol linoleate (CLL), the most favorable substrate for BsChe. The results showed that the sterol moieties of reasonable CLL docking poses localized to a specific active-site cleft surface formed by Leu266 and Ile287, which are unconserved among Burkholderia Lip homologs. Site-directed mutagenesis identified these residues as essential for the Che activity of BsChe, and Leu or Ile substitution conferred marked Che activity to Burkholderia Lips. In particular, Burkholderia cepacia and Burkholderia ubonensis Lips with the V266L/L287I double mutation exhibited ~50-fold and 500-fold higher Che activities than those of the wild-type enzymes, respectively. These results provide new insights into the substrate-binding mechanisms and selectivities of bacterial Lips.

Carboxylic Ester Hydrolases in Bacteria: Active Site, Structure, Function and Application

Carboxylic ester hydrolases (CEHs), which catalyze the hydrolysis of carboxylic esters to produce alcohol and acid, are identified in three domains of life. In the Protein Data Bank (PDB), 136 crystal structures of bacterial CEHs (424 PDB codes) from 52 genera and metagenome have been reported. In this review, we categorize these structures based on catalytic machinery, structure and substrate specificity to provide a comprehensive understanding of the bacterial CEHs. CEHs use Ser, Asp or water as a nucleophile to drive diverse catalytic machinery. The α/β/α sandwich architecture is most frequently found in CEHs, but 3-solenoid, β-barrel, up-down bundle, α/β/β/α 4-layer sandwich, 6 or 7 propeller and α/β barrel architectures are also found in these CEHs. Most are substrate-specific to various esters with types of head group and lengths of the acyl chain, but some CEHs exhibit peptidase or lactamase activities. CEHs are widely used in industrial applications, and are the objects of research in structure- or mutation-based protein engineering. Structural studies of CEHs are still necessary for understanding their biological roles, identifying their structure-based functions and structure-based engineering and their potential industrial applications.

A new benchmark illustrates that integration of geometric constraints inferred from enzyme reaction chemistry can increase enzyme active site modeling accuracy

Enzymes play a critical role in a wide array of industrial, medical, and research applications and with the recent explosion of genomic sequencing, we now have sequences for millions of enzymes for which there is no known structure. In order to utilize modern computational design tools for constructing inhibitors or engineering novel catalysts, the ability to accurately model enzymes is critical. A popular approach for modeling enzymes are comparative modeling techniques which can often accurately predict the global structural features. However, achieving atomic accuracy of an active site remains a challenge and is an issue when trying to utilize the molecular details for designing inhibitors or enhanced catalysts. Here we explore integrating knowledge about the required geometric orientation of conserved catalytic residues into the comparative modeling process in order to improve modeling accuracy. In order to investigate the utility of adding this information, we first carefully construct a benchmark set of reference structures to use. Consistent with previous findings, our benchmark demonstrates that the geometry between catalytic residues across an enzyme family is conserved and does not tend to deviate by more than 0.5Å. We then find that by integrating these geometric constraints during modeling, we can double the number of atomic level accuracy models (<1Å RMSD to the crystal structure ligand) within our benchmarking dataset, even for targets with templates as low as 20-30% sequence identity. Catalytic residues within an enzyme family are highly conserved and can often be readily identified through comparative sequence analysis to a known structure within the enzyme family. Therefore utilizing this readily available information has the potential to significantly improve drug design and enzyme engineering efforts for which there is no known structure for the enzyme of interest.

Lipases: An Overview

Lipases are ubiquitous enzymes, widespread in nature. They were first isolated from bacteria in the early nineteenth century, and the associated research continuously increased due to the characteristics of these enzymes. This chapter reviews the main sources, structural properties, and industrial applications of these highly studied enzymes.

The Role of Solvent-Accessible Leu-208 of Cold-Active Pseudomonas fluorescens Strain AMS8 Lipase in Interfacial Activation, Substrate Accessibility and Low-Molecular Weight Esterification in the Presence of Toluene

The alkaline cold-active lipase from Pseudomonas fluorescens AMS8 undergoes major structural changes when reacted with hydrophobic organic solvents. In toluene, the AMS8 lipase catalytic region is exposed by the moving hydrophobic lid 2 (Glu-148 to Gly-167). Solvent-accessible surface area analysis revealed that Leu-208, which is located next to the nucleophilic Ser-207 has a focal function in influencing substrate accessibility and flexibility of the catalytic pocket. Based on molecular dynamic simulations, it was found that Leu-208 strongly facilitates the lid 2 opening via its side-chain. The KM and Kcat/KM of L208A mutant were substrate dependent as it preferred a smaller-chain ester (pNP-caprylate) as compared to medium (pNP-laurate) or long-chain (pNP-palmitate) esters. In esterification of ethyl hexanoate, L208A promotes a higher ester conversion rate at 20 °C but not at 30 °C, as a 27% decline was observed. Interestingly, the wild-type (WT) lipase's conversion rate was found to increase with a higher temperature. WT lipase AMS8 esterification was higher in toluene as compared to L208A. Hence, the results showed that Leu-208 of AMS8 lipase plays an important role in steering a broad range of substrates into its active site region by regulating the flexibility of this region. Leu-208 is therefore predicted to be crucial for its role in interfacial activation and catalysis in toluene.

The Lid Domain in Lipases: Structural and Functional Determinant of Enzymatic Properties

Lipases are important industrial enzymes. Most of the lipases operate at lipid–water interfaces enabled by a mobile lid domain located over the active site. Lid protects the active site and hence responsible for catalytic activity. In pure aqueous media, the lid is predominantly closed, whereas in the presence of a hydrophobic layer, it is partially opened. Hence, the lid controls the enzyme activity. In the present review, we have classified lipases into different groups based on the structure of lid domains. It has been observed that thermostable lipases contain larger lid domains with two or more helices, whereas mesophilic lipases tend to have smaller lids in the form of a loop or a helix. Recent developments in lipase engineering addressing the lid regions are critically reviewed here. After on, the dramatic changes in substrate selectivity, activity, and thermostability have been reported. Furthermore, improved computational models can now rationalize these observations by relating it to the mobility of the lid domain. In this contribution, we summarized and critically evaluated the most recent developments in experimental and computational research on lipase lids.

How the Same Core Catalytic Machinery Catalyzes 17 Different Reactions: the Serine-Histidine-Aspartate Catalytic Triad of α/β-Hydrolase Fold Enzymes

Enzymes within a family often catalyze different reactions. In some cases, this variety stems from different catalytic machinery, but in other cases the machinery is identical; nevertheless, the enzymes catalyze different reactions. In this review, we examine the subset of α/β-hydrolase fold enzymes that contain the serine-histidine-aspartate catalytic triad. In spite of having the same protein fold and the same core catalytic machinery, these enzymes catalyze seventeen different reaction mechanisms. The most common reactions are hydrolysis of C-O, C-N and C-C bonds (Enzyme Classification (EC) group 3), but other enzymes are oxidoreductases (EC group 1), acyl transferases (EC group 2), lyases (EC group 4) or isomerases (EC group 5). Hydrolysis reactions often follow the canonical esterase mechanism, but eight variations occur where either the formation or cleavage of the acyl enzyme intermediate differs. The remaining eight mechanisms are lyase-type elimination reactions, which do not have an acyl enzyme intermediate and, in four cases, do not even require the catalytic serine. This diversity of mechanisms from the same catalytic triad stems from the ability of the enzymes to bind different substrates, from the requirements for different chemical steps imposed by these new substrates and, only in about half of the cases, from additional hydrogen bond partners or additional general acids/bases in the active site. This detailed analysis shows that binding differences and non-catalytic residues create new mechanisms and are essential for understanding and designing efficient enzymes.

Lipases: an overview

Lipases are ubiquitous enzymes, widespread in nature. They were first isolated from bacteria in the early nineteenth century and the associated research continuously increased due to the particular characteristics of these enzymes. This chapter reviews the main sources, structural properties, and industrial applications of these highly studied enzymes.

Solvent dielectric effect and side chain mutation on the structural stability of Burkholderia cepacia lipase active site: a quantum mechanical/molecular mechanics study

Quantum mechanical and molecular dynamics methods were used to analyze the structure and stability of neutral and zwitterionic configurations of the extracted active site sequence from a Burkholderia cepacia lipase, histidyl-seryl-glutamin (His86-Ser87-Gln88) and its mutated form, histidyl-cysteyl-glutamin (His86-Cys87-Gln88) in vacuum and different solvents. The effects of solvent dielectric constant, explicit and implicit water molecules and side chain mutation on the structure and stability of this sequence in both neutral and zwitterionic forms are represented. The quantum mechanics computations represent that the relative stability of zwitterionic and neutral configurations depends on the solvent structure and its dielectric constant. Therefore, in vacuum and the considered non-polar solvents, the neutral form of the interested sequences is more stable than the zwitterionic form, while their zwitterionic form is more stable than the neutral form in the aqueous solution and the investigated polar solvents in most cases. However, on the potential energy surfaces calculated, there is a barrier to proton transfer from the positively charged ammonium group to the negatively charged carboxylat group or from the ammonium group to the adjacent carbonyl oxygen and or from side chain oxygen and sulfur to negatively charged carboxylat group. Molecular dynamics simulations (MD) were also performed by using periodic boundary conditions for the zwitterionic configuration of the hydrated molecules in a box of water molecules. The obtained results demonstrated that the presence of explicit water molecules provides the more compact structures of the studied molecules. These simulations also indicated that side chain mutation and replacement of sulfur with oxygen leads to reduction of molecular flexibility and packing.

Modeling of solvent-dependent conformational transitions in Burkholderia cepacia lipase

BACKGROUND

The characteristic of most lipases is the interfacial activation at a lipid interface or in non-polar solvents. Interfacial activation is linked to a large conformational change of a lid, from a closed to an open conformation which makes the active site accessible for substrates. While for many lipases crystal structures of the closed and open conformation have been determined, the pathway of the conformational transition and possible bottlenecks are unknown. Therefore, molecular dynamics simulations of a closed homology model and an open crystal structure of Burkholderia cepacia lipase in water and toluene were performed to investigate the influence of solvents on structure, dynamics, and the conformational transition of the lid.

RESULTS

The conformational transition of B. cepacia lipase was dependent on the solvent. In simulations of closed B. cepacia lipase in water no conformational transition was observed, while in three independent simulations of the closed lipase in toluene the lid gradually opened during the first 10-15 ns. The pathway of conformational transition was accessible and a barrier was identified, where a helix prevented the lid from opening to the completely open conformation. The open structure in toluene was stabilized by the formation of hydrogen bonds.In simulations of open lipase in water, the lid closed slowly during 30 ns nearly reaching its position in the closed crystal structure, while a further lid opening compared to the crystal structure was observed in toluene. While the helical structure of the lid was intact during opening in toluene, it partially unfolded upon closing in water. The closing of the lid in water was also observed, when with eight intermediate structures between the closed and the open conformation as derived from the simulations in toluene were taken as starting structures. A hydrophobic beta-hairpin was moving away from the lid in all simulations in water, which was not observed in simulations in toluene. The conformational transition of the lid was not correlated to the motions of the beta-hairpin structure.

CONCLUSION

Conformational transitions between the experimentally observed closed and open conformation of the lid were observed by multiple molecular dynamics simulations of B. cepacia lipase. Transitions in both directions occurred without applying restraints or external forces. The opening and closing were driven by the solvent and independent of a bound substrate molecule.

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.

Crystal structures of TM0549 and NE1324--two orthologs of E. coli AHAS isozyme III small regulatory subunit

Crystal structures of two orthologs of the regulatory subunit of acetohydroxyacid synthase III (AHAS, EC 2.2.1.6) from Thermotoga maritima (TM0549) and Nitrosomonas europea (NE1324) were determined by single-wavelength anomalous diffraction methods with the use of selenomethionine derivatives at 2.3 A and 2.5 A, respectively. TM0549 and NE1324 share the same fold, and in both proteins the polypeptide chain contains two separate domains of a similar size. Each protein contains a C-terminal domain with ferredoxin-type fold and an N-terminal ACT domain, of which the latter is characteristic for several proteins involved in amino acid metabolism. The ferredoxin domain is stabilized by a calcium ion in the crystal structure of NE1324 and by a Mg(H2O)(6)2+ ion in TM0549. Both TM0549 and NE1324 form dimeric assemblies in the crystal lattice.

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.

Inhibition or activation of Pseudomonas species lipase by 1,2-ethylene-di-N-alkylcarbamates in detergents

1,2-Ethylene-di-N-n-propylcarbamate (1) is characterized as an essential activator of Pseudomonas species lipase while 1,2-ethylene-di-N-n-butyl-, t-butyl-, n-heptyl-, and n-octyl-carbamates (2-5) are characterized as the pseudo substrate inhibitors of the enzyme in the presence of the detergent taurocholate or triton X-100. The inhibition and activation reactions are more sensitive in taurocholate than in triton X-100. From CD studies, the enzyme changes conformations in the presence of the detergent and further alters conformations by addition of the carbamate activator or inhibitor into the enzyme-detergent adduct. Therefore, this study suggests that the conformational change of lipase during interfacial activation is a continuous process to expose the active site of the enzyme to substrate. From 600 MHz (1)H NMR studies, the conformations of the alpha- and beta-methylene moieties of the activator 1,2-ethylene-di-N-n-propylcarbamate in the presence of substrate change after adding taurocholate into the mixture, and the conformations of the beta-methylene moieties of the inhibitor 1,2-ethylene-di-N-n-butylcarbamate in the presence of substrate alter after adding taurocholate into the mixture.

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.

Cloning and expression of a novel esterase gene cpoA from Burkholderia cepacia

AIMS

To screen and clone a novel enzyme with specific activity for the resolution of (R)-beta-acetylmercaptoisobutyrate (RAM) from (R,S)-beta-acetylmercaptoisobutyrate [(R,S)-ester].

METHODS AND RESULTS

A micro-organism that produces a novel esterase was isolated and identified as the bacterium Burkholderia cepacia by using the analysis of cellular fatty acids, Biolog automated microbial identification/characterization system, and 16S rRNA gene sequence analysis. A novel esterase gene was cloned from the chromosomal DNA of B. cepacia and was designated as cpoA. The cpoA encodes a polypeptide of 273 amino acids which shows a strong sequence homology with many bacterial nonhaeme chloroperoxidases. In addition, a typical serine-hydrolase motif, Gly-X-Ser-X-Gly, and the highly conserved catalytic triad, Ser95, Asp224, and His253, were identified in the deduced amino acid sequence of cpoA by multiple sequence alignment.

CONCLUSION

The cpoA cloned from B. cepacia encodes a novel esterase which is highly related to the nonhaeme chloroperoxidases.

SIGNIFICANCE AND IMPACT OF THE STUDY

This is the first report that describes the isolation and cloning of a serine esterase gene from B. cepacia, which is useful in the chiral resolution of (R,S)-ester. The cloned gene will allow additional research on the bifunctionality of the enzyme with esterase and chloroperoxidase activity at the structural and functional levels.

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.