Esterases receive special attention because their wide distribution in biological systems and environments and their importance for physiology and chemical synthesis. The prediction of esterases substrate promiscuity level from sequence data and the molecular reasons why certain such enzymes are more promiscuous than others, remain to be elucidated. This limits the surveillance of the sequence space for esterases potentially leading to new versatile biocatalysts and new insights into their role in cellular function. Here we performed an extensive analysis of the substrate spectra of 145 phylogenetically and environmentally diverse microbial esterases, when tested with 96 diverse esters. We determined the primary factors shaping their substrate range by analyzing substrate range patterns in combination with structural analysis and protein-ligand simulations. We found a structural parameter that helps ranking (classifying) promiscuity level of esterases from sequence data at 94% accuracy. This parameter, the active site effective volume, exemplifies the topology of the catalytic environment by measuring the active site cavity volume corrected by the relative solvent accessible surface area (SASA) of the catalytic triad. Sequences encoding esterases with active site effective volumes (cavity volume/SASA) above a threshold show greater substrate spectra, which can be further extended in combination with phylogenetic data. This measure provides also a valuable tool for interrogating substrates capable of being converted. This measure, found to be transferred to phosphatases of the haloalkanoic acid dehalogenase superfamily and possibly other enzymatic systems, represents a powerful tool for low-cost bioprospecting for esterases with broad substrate ranges, in large scale sequence datasets.
Parabens (alkyl esters of 4-hydroxybenzoic acid) are widely used as preservatives in drugs, cosmetic products, and foodstuffs. Safety concerns have recently increased due to the potential health risks associated to exposure to large amounts of these substances. Biotransformation of parabens mainly includes hydrolysis of the ester bond and glucuronidation reactions. The hydrolysis and glucuronidation of a series of six parabens differing by the nature of the alkyl group were investigated in human liver microsomes and plasma, and the major human UDP-glucuronosyltransferase (UGT) isoforms involved in the reaction were identified. Methyl- and ethylparaben were stable in human plasma, with 95% of the initial concentration remaining after 24 h. On the other hand, propyl-, butyl- and benzylparaben concentrations decreased by 50% under similar conditions. In contrast, rapid hydrolysis was measured with human liver microsomes depending on the alkyl chain length, with t(1/2) varying from 22 min for methylparaben to 87 min for butylparaben. All parabens were actively glucuronidated by liver microsomes, in comparison to 4-hydroxybenzoic acid. They were mainly substrates of human recombinant UGT1A1, UGT1A8, UGT1A9, UGT2B7, UGT2B15 and UGT2B17. In conclusion, the parabens were readily metabolized in human liver through esterase hydrolysis and glucuronidation by several UGT isoforms. These results suggest that these parabens do not accumulate in human tissue.
Human liver has numerous hydrolytic enzymes involved in metabolism of endogenous and exogenous esters. Of these enzymes, carboxylesterases (EC 3.1.1.1) form an important group which hydrolyses many diverse ester substrates, including pro-ester drugs. Carboxylesterase activity was investigated in liver subcellular fractions from 22 individuals using the general carboxylesterase substrate phenylvalerate and the homologous series of esters methyl-, ethyl-, propyl-, butyl- and benzylparaben. The intra- and inter-individual variation in phenylvalerate and paraben metabolism was compared. Rates of hydrolysis were higher in microsomal fractions than cytosolic fractions for all compounds. The rate of paraben hydrolysis varied depending on the size of the paraben alcohol leaving group, showing a decrease with increasing leaving group size. Comparisons showed that individuals with high rates of hydrolysis towards methyl paraben also showed high rates of hydrolysis to the other parabens and phenylvalerate. Phenylvalerate as a non-specific carboxylesterase substrate was hydrolysed mainly by hCE1 in human livers and there was good correlation with small alcohol leaving group parabens, suggesting hCE1 involvement. Lower correlations with larger alcohol leaving group parabens are consistent with more hCE2 involvement.