Predicting the binding mode of carbocations produced in sesquiterpene synthase enzymes is not unlike finding a piece of hay in a haystack. A new method for tackling this problem is described.

The ability to biosynthetically produce chemicals beyond what is commonly found in Nature requires the discovery of novel enzyme function. Here we utilize two approaches to discover enzymes that enable specific production of longer-chain (C5–C8) alcohols from sugar. The first approach combines bioinformatics and molecular modelling to mine sequence databases, resulting in a diverse panel of enzymes capable of catalysing the targeted reaction. The median catalytic efficiency of the computationally selected enzymes is 75-fold greater than a panel of naively selected homologues. This integrative genomic mining approach establishes a unique avenue for enzyme function discovery in the rapidly expanding sequence databases. The second approach uses computational enzyme design to reprogramme specificity. Both approaches result in enzymes with >100-fold increase in specificity for the targeted reaction. When enzymes from either approach are integrated in vivo, longer-chain alcohol production increases over 10-fold and represents >95% of the total alcohol products.

Autophagy is important for degradation and recycling of intracellular components. In a diversity of genera and species, orthologs and paralogs of the yeast Atg4 and Atg8 proteins are crucial in the biogenesis of double-membrane autophagosomes that carry the cellular cargoes to vacuoles and lysosomes. Although many plant genome sequences are available, the ATG4 and ATG8 sequence analysis is limited to some model plants. We identified 28 ATG4 and 116 ATG8 genes from the available 18 different plant genome sequences. Gene structures and protein domain sequences of ATG4 and ATG8 are conserved in plant lineages. Phylogenetic analyses classified ATG8s into 3 subgroups suggesting divergence from the common ancestor. The ATG8 expansion in plants might be attributed to whole genome duplication, segmental and dispersed duplication, and purifying selection. Our results revealed that the yeast Atg4 processes Arabidopsis ATG8 but not human LC3A (HsLC3A). In contrast, HsATG4B can process yeast and plant ATG8s in vitro but yeast and plant ATG4s cannot process HsLC3A. Interestingly, in Nicotiana benthamiana plants the yeast Atg8 is processed compared to HsLC3A. However, HsLC3A is processed when coexpressed with HsATG4B in plants. Molecular modeling indicates that lack of processing of HsLC3A by plant and yeast ATG4 is not due to lack of interaction with HsLC3A. Our in-depth analyses of ATG4 and ATG8 in the plant lineage combined with results of cross-kingdom ATG8 processing by ATG4 further support the evolutionarily conserved maturation of ATG8. Broad ATG8 processing by HsATG4B and lack of processing of HsLC3A by yeast and plant ATG4s suggest that the cross-kingdom ATG8 processing is determined by ATG8 sequence rather than ATG4.

Terpene synthases comprise a family of enzymes that convert acyclic oligo-isoprenyl diphosphates to terpene natural products with complex, polycyclic carbon backbones via the generation and protection of carbocation intermediates. To accommodate this chemistry, terpene synthase active sites generally are lined with alkyl and aromatic, i.e., nonpolar, sidechains. Predicting the correct, mechanistically relevant binding modes for entire terpene synthase reaction pathways remains an unsolved challenge. Here we describe a method for identifying such modes: , a series of protocols to predict the orientation of carbon skeletons of substrates and derived carbocations relative to the bound diphosphate group in terpene synthase active sites. Using this recipe for bornyl diphosphate synthase, we have predicted binding modesthat are consistent with all current experimental observations, including the results of isotope labeling experiments and known stereoselectivity. In addition, the predicted binding modes recapitulate key findings of a seminal study involving more computationally demanding QM/MM molecular dynamics methods on part of this pathway. This work illustrates the value of the approach as a starting point for more involved calculations and sets the stage for the rational engineering of this family of enzymes.