Genome-wide analysis of GDSL-type esterases/lipases in <Emphasis Type="Italic">Arabidopsis</Emphasis>

Key message

In this present study, we introduce a fundamental framework and provide information regarding the possible roles of GDSL-type esterase/lipase gene family in Arabidopsis.

Abstract

GDSL-type esterases/lipases are hydrolytic enzymes with multifunctional properties such as broad substrate specificity, regiospecificity, and stereoselectivity. In this study, we identified 105 GDSL-type esterase/lipase genes in Arabidopsis thaliana by conducting a comprehensive computational analysis. Expression studies indicated that GDSL-type lipase proteins showed varied expression patterns. Phylogenetic tree analysis indicated that AtGELP (Arabidopsis thaliana GDSL-type esterase/lipase protein) gene family was divided into four clades. The phylogenetic analysis, combined with protein motif architectures, and expression profiling were used to predict the roles AtGELP genes. To investigate the physical roles of the AtGELP gene family, we successfully screened 88 AtGELP T-DNA knockout lines for 54 AtGELP genes from 199 putative SALK T-DNA mutants. Transgenic plants of AtGELP genes were used to elucidate the phenotypic characteristics in various developmental stages or stress conditions. Our results suggest that the AtGELP genes have diverse physical functions such as affecting the germination rate and early growth of seedlings subjected to high concentrations of glucose, or being involved in biotic stress responses.

     

GDSL family of serine esterases/lipases

GDSL esterases and lipases are hydrolytic enzymes with multifunctional properties such as broad substrate specificity and regiospecificity. They have potential for use in the hydrolysis and synthesis of important ester compounds of pharmaceutical, food, biochemical, and biological interests. This new subclass of lipolytic enzymes possesses a distinct GDSL sequence motif different from the GxSxG motif found in many lipases. Unlike the common lipases, GDSL enzymes do not have the so called nucleophile elbow. Studies show that GDSL hydrolases have a flexible active site that appears to change conformation with the presence and binding of the different substrates, much like the induced fit mechanism proposed by Koshland. Some of the GDSL enzymes have thioesterase, protease, arylesterase, and lysophospholipase activity, yet they appear to be the same protein with similar molecular weight (22–60 kDa for most esterases), although some have multiple glycosylation sites with higher apparent molecular weight. GDSL enzymes have five consensus sequence (I–V) and four invariant important catalytic residues Ser, Gly, Asn, and His in blocks I, II, III, and V, respectively. The oxyanion structure led to a new designation of these enzymes as SGNH-hydrolase superfamily or subfamily. Phylogenetic analysis revealed that block IIA which belonged to the SGNH-hydrolases was found only in clade I. Therefore, this family of hydrolases represents a new example of convergent evolution of lipolytic enzymes. These enzymes have little sequence homology to true lipases. Another important differentiating feature of GDSL subfamily of lipolytic enzymes is that the serine-containing motif is closer to the N-terminus unlike other lipases where the GxSxG motif is near the center. Since the first classification of these subclass or subfamily of lipases as GDSL(S) hydrolase, progress has been made in determining the consensus sequence, crystal structure, active site and oxyanion residues, secondary structure, mechanism of catalysis, and understanding the conformational changes. Nevertheless, much still needs to be done to gain better understanding of in vivo biological function, 3-D structure, how this group of enzymes evolved to utilize many different substrates, and the mechanism of reactions. Protein engineering is needed to improve the substrate specificity, enantioselectivity, specific activity, thermostability, and heterologous expression in other hosts (especially food grade microorganisms) leading to eventual large scale production and applications. We hope that this review will rekindle interest among researchers and the industry to study and find uses for these unique enzymes.     

High-throughput assays for lipases and esterases

In the past few years a considerable number of high-throughput screening (HTS) systems have been developed, especially for lipases and esterases. In this review, a range of HTS methods for the directed evolution of these hydrolases are covered. This includes spectrophotometric and fluorimetric formats as well as other approaches to allow for fast, efficient and reliable identification of desired enzyme variants within large mutant libraries. In addition, methods for library creation and application of lipases and esterases are briefly covered.     

Methods to increase enantioselectivity of lipases and esterases

Lipases and esterases are frequently used in the synthesis of optically pure compounds; however, natural enzymes do not always show sufficiently high enantioselectivity. Variation of the structure of the substrates, modification of the reaction system or protein engineering (e.g. the expression of pure enzymes, rational design or directed evolution) are strategies that can be employed to improve the distinction between two enantiomers or enantiotopic groups.     

Selectivity of lipases and esterases towards phenol esters

The selectivity of 28 lipases and esterases in the hydrolysis of butanoates of o-, m- or p-substituted phenols was investigated in a microtiterplate format. The phenols released during enzyme-catalyzed hydrolysis were converted in situ with Gibbs’ reagent to form a blue indophenol complex, which was quantified spectrophotometrically at 600 nm. Substantial differences in rates were found, which exhibits that the type and position of the substituent at the alkyl group has a strong influence on the selectivity of the enzymes. For various enzymes, the p-nitro derivative was the best substrate, whereas for other enzymes the m-Cl-derivative was preferentially hydrolyzed. Analysis of the data using the Hammett equation showed that sometimes the observed changes followed a predictable trend, but in several cases the result is very unexpected.     

Synthesis activity-based zymography for detection of lipases and esterases

A new zymography method for lipases and esterases was developed on the basis of the esterification reaction between fatty acids and alcohols. The enzymes were separated by SDS-PAGE and native PAGE. The gel was washed and then incubated in an aqueous solution containing fatty acids (oleic acid 18:1 or caprylic acid 8:0) and dodecanol. Synthesis was visualized by in situ precipitation of water-insoluble and non-diffusible fatty acid esters, such as dodecyl oleate and dodecyl octanoate. The synthesis activity-based zymography was confirmed with different enzyme samples, including commercial lipase preparations, purified recombinant lipase and cutinase, and crude culture supernatants of lipolytic enzyme-producing soil bacteria.     

The compositional patterns of the avian genomes and their evolutionary implications

The compositional distributions of large (main-band) DNA fragments from eight birds belonging to eight different orders (including both paleognathous and neognathous species) are very broad and extremely close to each other. These findings, which are paralleled by the compositional similarity of homologous coding sequences and their codon positions, support the idea that birds are a monophyletic group.The compositional distribution of third-codon positions of genes from chicken, the only avian species for which a relatively large number of coding sequences is known, is very broad and bimodal, the minor GC-richer peak reaching 100% GC. The very high compositional heterogeneity of avian genomes is accompanied (as in the case of mammalian genomes) by a very high speciation rate compared to cold-blooded vertebrates which are characterized by genomes that are much less heterogeneous. The higher GC levels attained by avian compared to mammalian genomes might be correlated with the higher body temperature (41–43°C) of birds compared to mammals (37°C).A comparison of GC levels of coding sequences and codon positions from man and chicken revealed very close average GC levels and standard deviations. Homologous coding sequences and codon positions from man and chicken showed a surprisingly high degree of compositional similarity which was, however, higher for GC-poor than for GC-rich sequences. This indicates that GC-poor isochores of warm-blooded vertebrates reflect the composition of the isochores of the genome of the common reptilian ancestor of mammals and birds, which underwent only a small compositional change at the transition from cold- to warm-blooded vertebrates. In contrast, the GC-rich isochores of birds and mammals are the result of large compositional changes at the same evolutionary transition, where were in part different in the two classes of warm-blooded vertebrates.Correspondence to: G. Bernaadi     

MELDB: a database for microbial esterases and lipases

MELDB is a comprehensive protein database of microbial esterases and lipases which are hydrolytic enzymes important in the modern industry. Proteins in MELDB are clustered into groups according to their sequence similarities based on a local pairwise alignment algorithm and a graph clustering algorithm (TribeMCL). This differs from traditional approaches that use global pairwise alignment and joining methods. Our procedure was able to reduce the noise caused by dubious alignment in the distantly related or unrelated regions in the sequences. In the database, 883 esterase and lipase sequences derived from microbial sources are deposited and conserved parts of each protein are identified. HMM profiles of each cluster were generated to classify unknown sequences. Contents of the database can be keyword-searched and query sequences can be aligned to sequence profiles and sequences themselves.     

A simple activity staining protocol for lipases and esterases

A simple activity staining protocol for rapid detection and differentiation of lipases and esterases was developed based on pH drop due to fatty acids released following lipolysis. Though the detection of lipolysis as a function of drop in pH is not new, the present method has been made more sensitive by the judicious selection of the initial pH of the chromogenic substrate, which has been set near the end point of the dye so that even a slight drop in pH results in immediate color change. In the present case, the dye phenol red was taken, which has the end point at pH 7.3–7.4 where the color is pink. A slight drop due to fatty acid release results in yellow coloration. The assay has high reproducibility and can detect as low as 0.5 p-NPP enzyme units within 15 min. In addition, this method can be used for various lipidic substrates such as oils and tributyrin, making it suitable for both lipases and esterases.     

Functional proteomic analysis of lipases and esterases in cultured human adipocytes

This study reports on the analysis of the lipolytic proteome of cultured human fat cells. We used specific affinity tags to detect and identify the lipolytic and esterolytic enzymes in human subcutaneous (Sc) and visceral (Visc) adipocytes. For this purpose, differentiated fat cells were incubated with a fluorescent suicide inhibitor followed by protein separation using one- or two-dimensional gel electrophoresis. After detection by fluorescence laser scanning, the labeled proteins were tryptically digested and peptides were identified by mass spectrometry. In addition, a biotinylated probe was used for specific enzyme labeling with subsequent avidin affinity isolation of the tagged proteins. Finally, we determined the quantitative differences in protein expression levels between subcutaneous and visceral adipocytes using differential activity-based gel electrophoresis (DABGE). We found that the lipase/esterase patterns of both cell types are very similar, except for some proteins that were only found in Sc cells. Two novel enzyme candidates identified in this study were overexpressed and characterized using biologically relevant glycerolipid substrates in vitro. Both of them showed pronounced hydrolytic activities on hydrophobic acylglycerols and therefore may be considered lipases. The physiological functions of the novel lipolytic proteins in vivo are currently subject to investigation.     

How do lipases and esterases work: the electrostatic contribution

This work explores the role of one of the factors explaining lipase/esterase activity: the contribution of electrostatic interactions to lipase/esterase activity. The electrostatic potential distribution on the molecular surface of an enzyme as a function of pH determines, to a large extent, the enzyme's pH activity profile. Other important factors include the presence and distribution of polar and hydrophobic residues in the active cleft. We have mapped the electrostatic potential distribution as a function of pH on the molecular surface of nine lipases/esterases for which the 3D structure is experimentally known. A comparison of these potential maps at different pH values with the corresponding pH-activity profile, pH optimum or pH range where the activity displayed by the enzyme is maximum, has revealed a considerable correlation. A negative potential in the active site appears correlated with maximum activity towards triglycerides, which has prompted us to propose a model for product release ('The electrostatic catapult model') after cleavage of an ester bond. At the same time as the bottom of the active site cleft becomes negatively charged, other nearby regions also titrate and become negatively charged when pH becomes more alkaline, for some of the studied lipases. If such lipases also show phospholipase activity (such as guinea pig lipase-related proteins 2 chimera) we raise the hypothesis that such other titratable regions after becoming negatively charged might stabilise the positive charge present in the polar head of phospholipids, such as phosphatidylcholine, phosphatidylethanolamine and phosphatidylserine. The distribution of polar, weak polar and non-polar residues on the molecular surface of each studied lipase, in particular the active site region, was compared for all the lipases studied. The combination of graphical visualisation of the electrostatic potential maps and the polarity maps combined with knowledge about the location of key residues on the protein surface allows us to envision atomic models for lipolytic activity.     

Gel filtration applied to the study of lipases and other esterases

1. Sephadex G-100 and G-200 gel-filtration columns were calibrated for molecular-weight estimation with proteins of known molecular weights, and used to study the composition of several lipase or esterase preparations. 2. Enzymes from cow's milk, rat adipose tissue and pig pancreas were detected in the column effluents by their ability to liberate free acid from emulsified tributyrin at pH 8·5. 3. Four tributyrinases were detected in preparations from individual cow's milks. Molecular weights 62000, 75000 and 112000 were estimated for three of them, but although the fourth may be of unusually low molecular weight an estimate was not possible. 4. Extracts of rat adipose tissue apparently contained six tributyrinases (molecular weights 39000, 47000, 55000, 68000, 75000 and 200000) but the relative amounts of these enzymes varied widely from rat to rat. 5. Tributyrinase activity in juice expressed from pig pancreatic tissue was due mainly to one enzyme (molecular weight 42000). On the other hand, activity in extracts of acetone-dried pancreas was confined to material of molecular weight > 10⁶, which may be an aggregated form of the lower-molecular-weight enzyme. 6. Activity in fractionated wheat-germ extracts was assayed with emulsified triacetin substrate, and was evidently due to one enzyme (molecular weight 51000). 7. Some problems arising in the application of gel filtration to the study of lipase–esterase systems were indicated.     

Modelling substrate specificity and enantioselectivity for lipases and esterases by substrate-imprinted docking

Background  

Previously, ways to adapt docking programs that were developed for modelling inhibitor-receptor interaction have been explored. Two main issues were discussed. First, when trying to model catalysis a reaction intermediate of the substrate is expected to provide more valid information than the ground state of the substrate. Second, the incorporation of protein flexibility is essential for reliable predictions.     

Phylogenetic and evolutionary analysis of NBS-encoding genes in Rosaceae fruit crops

Phylogenetic relationships of the nucleotide binding site (NBS)-encoding resistance gene homologues (RGHs) among 12 species in five genera of Rosaceae fruit crops were evaluated. A total of 228 Rosaceous RGHs were deeply separated into two distinct clades, designated as TIR (sequences within this clade containing a Toll Interleukin-1 Receptor domain) and NonTIR (sequences lacking a TIR domain). Most Rosaceous RGH genes were phylogenetically distinct from Arabidopsis, Rice or Pine genes, except for a few Rosaceous members which grouped closely with Arabidopsis genes. Within Rosaceae, sequences from multiple species were often phylogenetically clustered together, forming heterogenous groups, however, apple- and chestnut rose-specific groups really exist. Gene duplication followed by sequence divergence were proposed as the mode for the evolution of a large number of distantly or closely related RGH genes in Rosaceae, and this mode may play a role in the generation of new resistance specificity. Positively selected sites within NBS-coding region were detected and thus nucleotide variation within NBS domain may function in determining disease resistance specificity. This study also discusses the synteny of a genomic region that encompass powdery mildew resistance locus among Malus, Prunus and Rosa, which may have potential use for fruit tree disease breeding and important gene cloning.     

Polymer membrane ion-selective electrodes as a convenient tool for lipases and esterases assays

We have proposed a novel assay for lipases and esterases activity determination based on potentiometry with ion-selective electrodes (ISEs). Enzyme preparations, obtained from the living cells, are complex mixtures of various proteins, short peptides, lipids, carbohydrates, and other compounds. The most commonly used quantitative methods in enzyme studies are based on spectrophotometric or spectroflourimetric protocols which has significant limitations. They are not valid for samples that are turbid or strongly colored. To overcome those drawbacks we have proposed an assay based on potentiometry with ISEs for lipases and esterases activity determination. This electrochemical methodology represents an attractive tool for enzyme analysis, because of its low detection limit, independence from sample volume and from sample turbidity. The usefulness of this assay has been proven by the determination of the activity of various raw enzymes “acetone powders” isolated from animal tissues. Moreover, activities of fractions obtained during purification of one of those raw biocatalysts were also determined that way. The reliability of determination enzyme activity with ISE assay was proven by comparison with a classical spectrophotometric method.     

Evaluation of droplet-based microfluidic platforms as a convenient tool for lipases and esterases assays

Abstract

The accurate estimation of kinetic parameters is of fundamental importance for biochemical studies for research and industry. In this paper, we demonstrate the application of a modular microfluidic system for execution of enzyme assays that allow determining the kinetic parameters of the enzymatic reactions such as V_max – the maximum rate of reaction and K_M – the Michaelis constant. For experiments, the fluorogenic carbonate as a probe for a rapid determination of the kinetic parameters of hydrolases, such as lipases and esterases, was used. The microfluidic system together with the method described yields the kinetic constants calculated from the concentration of enzymatic product changes via a Michaelis–Menten model using the Lambert function W(x). This modular microfluidic system was validated on three selected enzymes (hydrolases).     

Release of short chain fatty acids from cream lipids by commercial lipases and esterases

Lipases and esterases are frequently used in dairy production processes to enhance the buttery flavour of the end product. Short chain fatty acids, and especially butanoic acid, play a key role in this and different enzymes with specificity towards short chain fatty acids are commercially available as potent flavouring tools. We have compared six lipases/esterases associated with buttery flavour production. Although specificity to short chain fatty acids was ascribed to each enzyme, clear differences in free fatty acid profiles were found when these enzymes were applied on cream. Candida cylindraceae lipase was the most useful enzyme for buttery flavour production in cream with the highest yield of free fatty acids (57 g oleic acid 100 g⁻¹ fat), no release of long chain fatty acids and specificity towards butanoic acid.