Structure and evolutionary trace-assisted screening of a residue swapping the substrate ambiguity and chiral specificity in an esterase

Transforming an esterase into an enantioselective catecholase through bioconjugation of a versatile metal-chelating inhibitor

Metal complexes introduced into protein scaffolds can generate versatile biomimetic catalysts endowed with a variety of catalytic properties. Here, we synthesized and covalently bound a bipyridinyl derivative to the active centre of an esterase to generate a biomimetic catalyst that shows catecholase activity and enantioselective catalytic oxidation of (+)-catechin.

The Mobility of the Cap Domain Is Essential for the Substrate Promiscuity of a Family IV Esterase from Sorghum Rhizosphere Microbiome

Metagenomics offers the possibility to screen for versatile biocatalysts. In this study, the microbial community of the Sorghum bicolor rhizosphere was spiked with technical cashew nut shell liquid, and after incubation, the environmental DNA (eDNA) was extracted and subsequently used to build a metagenomic library. We report the biochemical features and crystal structure of a novel esterase from the family IV, EH₀, retrieved from an uncultured sphingomonad after a functional screen in tributyrin agar plates. EH₀ (optimum temperature [T_opt], 50°C; melting temperature [T_m], 55.7°C; optimum pH [pH_opt], 9.5) was stable in the presence of 10 to 20% (vol/vol) organic solvents and exhibited hydrolytic activity against p-nitrophenyl esters from acetate to palmitate, preferably butyrate (496 U mg^-1), and a large battery of 69 structurally different esters (up to 30.2 U mg^-1), including bis(2-hydroxyethyl)-terephthalate (0.16 ± 0.06 U mg^-1). This broad substrate specificity contrasts with the fact that EH₀ showed a long and narrow catalytic tunnel, whose access appears to be hindered by a tight folding of its cap domain. We propose that this cap domain is a highly flexible structure whose opening is mediated by unique structural elements, one of which is the presence of two contiguous proline residues likely acting as possible hinges, which together allow for the entrance of the substrates. Therefore, this work provides a new role for the cap domain, which until now was thought to be an immobile element that contained hydrophobic patches involved in substrate prerecognition and in turn substrate specificity within family IV esterases. IMPORTANCE A better understanding of structure-function relationships of enzymes allows revelation of key structural motifs or elements. Here, we studied the structural basis of the substrate promiscuity of EH₀, a family IV esterase, isolated from a sample of the Sorghum bicolor rhizosphere microbiome exposed to technical cashew nut shell liquid. The analysis of EH₀ revealed the potential of the sorghum rhizosphere microbiome as a source of enzymes with interesting properties, such as pH and solvent tolerance and remarkably broad substrate promiscuity. Its structure resembled those of homologous proteins from mesophilic Parvibaculum and Erythrobacter spp. and hyperthermophilic Pyrobaculum and Sulfolobus spp. and had a very narrow, single-entry access tunnel to the active site, with access controlled by a capping domain that includes a number of nonconserved proline residues. These structural markers, distinct from those of other substrate-promiscuous esterases, can help in tuning substrate profiles beyond tunnel and active site engineering.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Odd one out? Functional tuning of Zymomonas mobilis pyruvate kinase is narrower than its allosteric, human counterpart

Various protein properties are often illuminated using sequence comparisons of protein homologs. For example, in analyses of the pyruvate kinase multiple sequence alignment, the set of positions that changed during speciation ("phylogenetic" positions) were enriched for "rheostat" positions in human liver pyruvate kinase (hLPYK). (Rheostat positions are those which, when substituted with various amino acids, yield a range of functional outcomes). However, the correlation was moderate, which could result from multiple biophysical constraints acting on the same position during evolution and/or various sources of noise. To further examine this correlation, we here tested Zymomonas mobilis PYK (ZmPYK), which has <65% sequence identity to any other PYK sequence. Twenty-six ZmPYK positions were selected based on their phylogenetic scores, substituted with multiple amino acids, and assessed for changes in Kapp-PEP . Although we expected to identify multiple, strong rheostat positions, only one moderate rheostat position was detected. Instead, nearly half of the 271 ZmPYK variants were inactive and most others showed near wild-type function. Indeed, for the active ZmPYK variants, the total range of Kapp,PEP values ("tunability") was 40-fold less than that observed for hLPYK variants. The combined functional studies and sequence comparisons suggest that ZmPYK has evolved functional and/or structural attributes that differ from the rest of the family. We hypothesize that including such "orphan" sequences in MSA analyses obscures the correlations used to predict rheostat positions. Finally, results raise the intriguing biophysical question as to how the same protein fold can support rheostat positions in one homolog but not another.

EPIMUTESTR: a nearest neighbor machine learning approach to predict cancer driver genes from the evolutionary action of coding variants

Discovering rare cancer driver genes is difficult because their mutational frequency is too low for statistical detection by computational methods. EPIMUTESTR is an integrative nearest-neighbor machine learning algorithm that identifies such marginal genes by modeling the fitness of their mutations with the phylogenetic Evolutionary Action (EA) score. Over cohorts of sequenced patients from The Cancer Genome Atlas representing 33 tumor types, EPIMUTESTR detected 214 previously inferred cancer driver genes and 137 new candidates never identified computationally before of which seven genes are supported in the COSMIC Cancer Gene Census. EPIMUTESTR achieved better robustness and specificity than existing methods in a number of benchmark methods and datasets.