Esterases receive special attention because of 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 rank (classify) the 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 data sets.

No abstract

Iterative saturation mutagenesis was performed over six residues delimiting the substrate binding pocket of a high redox potential chimeric laccase with the aim of enhancing its activity over sinapic acid, a ligninrelated phenol of industrial interest. In total, more than 15000 clones were screened and two selected variants, together with the parent-type laccase, were purified and characterized. The new variants presented shifted pH activity profiles and enhanced turnover rates on sinapic acid and its methyl ester, whereas the oxidation of related phenols was not significantly enhanced. Neither the enzyme's redox potential nor the mechanism of the reaction was affected. Quantum mechanics and molecular dynamics calculations were done to rationalize the effect of the selected mutations, revealing the critical role of the residues of the enzyme pocket to provide the precise binding of the substrate that enables an efficient electron transfer to the T1 copper. The results presented highlight the power of combining directed evolution and computational approaches to give novel solutions in enzyme engineering and to understand the mechanistic reasons behind them, offering new insights for further rational design towards specific targets.This work was funded by INDOX (KBBE-2013-7-613549) European project and NOESIS (BIO2014-56388-R) and CTQ2013-\ud 48287-R Spanish National Projects. I. P. and G. S. acknowledge the Spanish Research Council (CSIC) and MINECO for their respective predoctoral fellowships.Peer ReviewedPostprint (published version

From laccase design to application of the overexpressed biocatalyst in an industrial environment for eco-friendly synthesis of polyaniline and dyes.

Effects of altering the properties of an active site in an enzymatic homogeneous catalyst have been extensively reported. However, the possibility of increasing the number of such sites, as commonly done in heterogeneous catalytic materials, remains unexplored, particularly because those have to accommodate appropriate residues in specific configurations. This possibility was investigated by using a serine ester hydrolase as the target enzyme. By using the Protein Energy Landscape Exploration software, which maps ligand diffusion and binding, we found a potential binding pocket capable of holding an extra catalytic triad and oxyanion hole contacts. By introducing two mutations, this binding pocket became a catalytic site. Its substrate specificity, substrate preference, and catalytic activity were different from those of the native site of the wild type ester hydrolase and other hydrolases, due to the differences in the active site architecture. Converting the binding pocket into an extra catalytic active site was proven to be a successful approach to create a serine ester hydrolase with two functional reactive groups. Our results illustrate the accuracy and predictive nature of modern modeling techniques, opening novel catalytic opportunities coming from the presence of different catalytic environments in single enzymes.