Moving Through Barriers in Science and Life.

This first serious attempt at an autobiographical accounting has forced me to sit still long enough to compile my thoughts about a long personal and scientific journey. I especially hope that my trajectory will be of interest and perhaps beneficial to much younger women who are just getting started in their careers. To paraphrase from Virginia Woolf's writings in A Room of One's Own at the beginning of the 20th century, "for most of history Anonymous was a Woman." However, Ms. Woolf is also quoted as saying "nothing has really happened until it has been described," a harbinger of the enormous historical changes that were about to be enacted and recorded by women in the sciences and other disciplines. The progress in my chosen field of study-the chemical basis of enzyme action-has also been remarkable, from the first description of an enzyme's 3D structure to a growing and deep understanding of the origins of enzyme catalysis.

Directed Evolution of Protein Catalysts.

Directed evolution is a powerful technique for generating tailor-made enzymes for a wide range of biocatalytic applications. Following the principles of natural evolution, iterative cycles of mutagenesis and screening or selection are applied to modify protein properties, enhance catalytic activities, or develop completely new protein catalysts for non-natural chemical transformations. This review briefly surveys the experimental methods used to generate genetic diversity and screen or select for improved enzyme variants. Emphasis is placed on a key challenge, namely how to generate novel catalytic activities that expand the scope of natural reactions. Two particularly effective strategies, exploiting catalytic promiscuity and rational design, are illustrated by representative examples of successfully evolved enzymes. Opportunities for extending these approaches to more complex biocatalytic systems are also considered.

Perturbative diffraction methods resolve a conformational switch that facilitates a two-step enzymatic mechanism.

Enzymes catalyze biochemical reactions through precise positioning of substrates, cofactors, and amino acids to modulate the transition-state free energy. However, the role of conformational dynamics remains poorly understood due to poor experimental access. This shortcoming is evident with Escherichia coli dihydrofolate reductase (DHFR), a model system for the role of protein dynamics in catalysis, for which it is unknown how the enzyme regulates the different active site environments required to facilitate proton and hydride transfer. Here, we describe ligand-, temperature-, and electric-field-based perturbations during X-ray diffraction experiments to map the conformational dynamics of the Michaelis complex of DHFR. We resolve coupled global and local motions and find that these motions are engaged by the protonated substrate to promote efficient catalysis. This result suggests a fundamental design principle for multistep enzymes in which pre-existing dynamics enable intermediates to drive rapid electrostatic reorganization to facilitate subsequent chemical steps.

Engineering substrate channeling in a bifunctional terpene synthase.

Fusicoccadiene synthase from Phomopsis amygdala (PaFS) is a bifunctional terpene synthase. It contains a prenyltransferase (PT) domain that generates geranylgeranyl diphosphate (GGPP) from dimethylallyl diphosphate and three equivalents of isopentenyl diphosphate, and a cyclase domain that converts GGPP into fusicoccadiene, a precursor of the diterpene glycoside Fusicoccin A. The two catalytic domains are connected by a flexible 69-residue linker. The PT domain mediates oligomerization to form predominantly octamers, with cyclase domains randomly splayed out around the PT core. Surprisingly, despite the random positioning of cyclase domains, substrate channeling is operative in catalysis since most of the GGPP generated by the PT remains on the enzyme for cyclization. Here, we demonstrate that covalent linkage of the PT and cyclase domains is not required for GGPP channeling, although covalent linkage may improve channeling efficiency. Moreover, GGPP competition experiments with other diterpene cyclases indicate that the PaFS PT and cyclase domains are preferential partners regardless of whether they are covalently linked or not. The cryoelectron microscopy structure of the 600-kD "linkerless" construct, in which the 69-residue linker is spliced out and replaced with the tripeptide PTQ, reveals that cyclase pairs associate with all four sides of the PT octamer and exhibit fascinating quaternary structural flexibility. These results suggest that optimal substrate channeling is achieved when a cyclase domain associates with the side of the PT octamer, regardless of whether the two domains are covalently linked and regardless of whether this interaction is transient or locked in place.

Enhancing luciferase activity and stability through generative modeling of natural enzyme sequences.

The availability of natural protein sequences synergized with generative AI provides new paradigms to engineer enzymes. Although active enzyme variants with numerous mutations have been designed using generative models, their performance often falls short of their wild type counterparts. Additionally, in practical applications, choosing fewer mutations that can rival the efficacy of extensive sequence alterations is usually more advantageous. Pinpointing beneficial single mutations continues to be a formidable task. In this study, using the generative maximum entropy model to analyze Renilla luciferase (RLuc) homologs, and in conjunction with biochemistry experiments, we demonstrated that natural evolutionary information could be used to predictively improve enzyme activity and stability by engineering the active center and protein scaffold, respectively. The success rate to improve either luciferase activity or stability of designed single mutants is ~50%. This finding highlights nature's ingenious approach to evolving proficient enzymes, wherein diverse evolutionary pressures are preferentially applied to distinct regions of the enzyme, ultimately culminating in an overall high performance. We also reveal an evolutionary preference in RLuc toward emitting blue light that holds advantages in terms of water penetration compared to other light spectra. Taken together, our approach facilitates navigation through enzyme sequence space and offers effective strategies for computer-aided rational enzyme engineering.

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology.

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

Redox Biochemistry of the Genetic Code.

New findings on the chemistry of the amino acids, their role in protein folding, and their sequential primordial introduction have uncovered concealed causalities in genetic code evolution. The genetically encoded amino acids successively provided (i) membrane anchors, (ii) halophilic protein folds, (iii) mesophilic protein folds, (iv) metal ligation, and (v) antioxidation.

"Mix and match" auto-assembly of glycosyltransferase domains delivers biocatalysts with improved substrate promiscuity.

Glycosyltransferases (GT) catalyze the glycosylation of bioactive natural products, including peptides and proteins, flavonoids, and sterols, and have been extensively used as biocatalysts to generate glycosides. However, the often narrow substrate specificity of wild-type GTs requires engineering strategies to expand it. The GT-B structural family is constituted by GTs that share a highly conserved tertiary structure in which the sugar donor and acceptor substrates bind in dedicated domains. Here, we have used this selective binding feature to design an engineering process to generate chimeric glycosyltransferases that combine auto-assembled domains from two different GT-B enzymes. Our approach enabled the generation of a stable dimer with broader substrate promiscuity than the parent enzymes that were related to relaxed interactions between domains in the dimeric GT-B. Our findings provide a basis for the development of a novel class of heterodimeric GTs with improved substrate promiscuity for applications in biotechnology and natural product synthesis.

Structure of methyltransferase RedM that forms the dimethylpyrrolinium of the bisindole reductasporine.

Bisindoles are biologically active natural products that arise from the oxidative dimerization of two molecules of l-tryptophan. In bacterial bisindole pathways, a core set of transformations is followed by the action of diverse tailoring enzymes that catalyze reactions that lead to diverse bisindole products. Among bisindoles, reductasporine is distinct due to its dimethylpyrrolinium structure. Its previously reported biosynthetic gene cluster encodes two unique tailoring enzymes, the imine reductase RedE and the dimethyltransferase RedM, which were shown to produce reductasporine from a common bisindole intermediate in recombinant E. coli. To gain more insight into the unique tailoring enzymes in reductasporine assembly, we reconstituted the biosynthetic pathway to reductasporine in vitro and then solved the 1.7 Å resolution structure of RedM. Our work reveals RedM adopts a variety of conformational changes with distinct open and closed conformations, and site-directed mutagenesis alongside sequence analysis identifies important active site residues. Finally, our work sets the stage for understanding how RedM evolved to react with a pyrrolinium scaffold and may enable the development of new dimethyltransferase catalysts.

An imine reductase that captures reactive intermediates in the biosynthesis of the indolocarbazole reductasporine.

Imine reductases (IREDs) and reductive aminases have been used in the synthesis of chiral amine products for drug manufacturing; however, little is known about their biological contexts. Here we employ structural studies and site-directed mutagenesis to interrogate the mechanism of the IRED RedE from the biosynthetic pathway to the indolocarbazole natural product reductasporine. Cocrystal structures with the substrate-mimic arcyriaflavin A reveal an extended active site cleft capable of binding two indolocarbazole molecules. Site-directed mutagenesis of a conserved aspartate in the primary binding site reveals a new role for this residue in anchoring the substrate above the NADPH cofactor. Variants targeting the secondary binding site greatly reduce catalytic efficiency, while accumulating oxidized side-products. As indolocarbazole biosynthetic intermediates are susceptible to spontaneous oxidation, we propose the secondary site acts to protect against autooxidation, and the primary site drives catalysis through precise substrate orientation and desolvation effects. The structure of RedE with its extended active site can be the starting point as a new scaffold for engineering IREDs and reductive aminases to intercept large substrates relevant to industrial applications.

Enzyme-cargo encapsulation peptides bind between tessellating tiles of the bacterial microcompartment shell.

Bacterial microcompartments are prokaryotic organelles comprising encapsulated enzymes within a thin protein shell. They facilitate metabolic processing including propanediol, choline, glycerol, and ethanolamine utilization, and they accelerate carbon fixation in cyanobacteria. Enzymes targeted to the inside of the microcompartment frequently possess a cargo-encapsulation peptide, but the site to which the peptide binds is unclear. We provide evidence that the encapsulation peptides bind to the hydrophobic groove formed between tessellating subunits of the shell proteins. In silico docking studies provide a compelling model of peptide binding to this prominent hydrophobic groove. This result is consistent with the now widely accepted view that the convex side of the shell oligomers faces the lumen of the microcompartment. The binding of the encapsulation peptide to the groove between tessellating shell protein tiles explains why it has been difficult to define the peptide binding site using other methods, provides a mechanism by which encapsulation-peptide bearing enzymes can promote shell assembly, and explains how the presence of cargo affects the size and shape of the bacterial microcompartment. This knowledge may be exploited in engineering microcompartments or disease prevention by hampering cargo encapsulation.

Cofactor Binding Triggers Rapid Conformational Remodelling of the Active Site of a Methyltransferase Ribozyme.

The methyltransferase ribozyme SMRZ-1 utilizes S-adenosyl-methionine (SAM) and Cu (II) ions to methylate RNA. Comparison of the SAM bound and unbound RNA structures has shown a conformational change in the RNA. However, the contribution of specific interactions and the role of a pseudo-triplex motif in the catalytic centre on the methylation reaction is not completely understood. In this study, we have used atomic substitutions and mutational analysis to investigate the reaction specificity and the key interactions required for catalysis. Substitution of the fluorescent nucleotide 2-aminopurine within the active ribozyme enabled the conformational dynamics of the RNA upon co-factor binding to be explored using fluorescence spectroscopy. We show that fast co-factor binding (t1/2 ∼ 0.7 seconds) drives a conformational change in the RNA to facilitate methyl group transfer. The importance of stacking interactions at the pseudo-triplex motif and chelation of the Cu (II) ion were shown to be essential for SAM binding.

Klebsiella pneumoniae carbapenemase variant 44 acquires ceftazidime-avibactam resistance by altering the conformation of active-site loops.

Klebsiella pneumoniae carbapenemase 2 (KPC-2) is an important source of drug resistance as it can hydrolyze and inactivate virtually all ß-lactam antibiotics. KPC-2 is potently inhibited by avibactam via formation of a reversible carbamyl linkage of the inhibitor with the catalytic serine of the enzyme. However, the use of avibactam in combination with ceftazidime (CAZ-AVI) has led to the emergence of CAZ-AVI-resistant variants of KPC-2 in clinical settings. One such variant, KPC-44, bears a 15 amino acid duplication in one of the active-site loops (270-loop). Here, we show that the KPC-44 variant exhibits higher catalytic efficiency in hydrolyzing ceftazidime, lower efficiency toward imipenem and meropenem, and a similar efficiency in hydrolyzing ampicillin, than the WT KPC-2 enzyme. In addition, the KPC-44 variant enzyme exhibits 12-fold lower AVI carbamylation efficiency than the KPC-2 enzyme. An X-ray crystal structure of KPC-44 showed that the 15 amino acid duplication results in an extended and partially disordered 270-loop and also changes the conformation of the adjacent 240-loop, which in turn has altered interactions with the active-site omega loop. Furthermore, a structure of KPC-44 with avibactam revealed that formation of the covalent complex results in further disorder in the 270-loop, suggesting that rearrangement of the 270-loop of KPC-44 facilitates AVI carbamylation. These results suggest that the duplication of 15 amino acids in the KPC-44 enzyme leads to resistance to CAZ-AVI by modulating the stability and conformation of the 270-, 240-, and omega-loops.

Komagataella phaffii Erp41 is a protein disulfide isomerase with unprecedented disulfide bond catalyzing activity when coupled to glutathione.

In the methylotrophic yeast Komagataella phaffii, we identified an endoplasmic reticulum-resident protein disulfide isomerase (PDI) family member, Erp41, with a peculiar combination of active site motifs. Like fungal ERp38, it has two thioredoxin-like domains which contain active site motifs (a and a'), followed by an alpha-helical ERp29c C-terminal domain (c domain). However, while the a domain has a typical PDI-like active site motif (CGHC), the a' domain instead has CGYC, a glutaredoxin-like motif which confers to the protein an exceptional affinity for GSH/GSSG. This combination of active site motifs has so far been unreported in PDI-family members. Homology searches revealed ERp41 is present in the genome of some plants, fungal parasites, and a few nonconventional yeasts, among which are Komagataella spp. and Yarrowia lipolytica. These yeasts are both used for the production of secreted recombinant proteins. Here, we analyzed the activity of K. phaffii Erp41. We report that it is nonessential in K. phaffii, and that it can catalyze disulfide bond formation in partnership with the sulfhydryl oxidase Ero1 in vitro with higher turnover rates than the canonical PDI from K. phaffii, Pdi1, but slower activation times. We show how Erp41 has unusually fast glutathione-coupled oxidation activity and relate it to its unusual combination of active sites in its thioredoxin-like domains. We further describe how this determines its unusually efficient catalysis of dithiol oxidation in peptide and protein substrates.

Structural and biochemical basis for regiospecificity of the flavonoid glycosyltransferase UGT95A1.

Glycosylation is a predominant strategy plants use to fine-tune the properties of small molecule metabolites to affect their bioactivity, transport, and storage. It is also important in biotechnology and medicine as many glycosides are utilized in human health. Small molecule glycosylation is largely carried out by family 1 glycosyltransferases. Here, we report a structural and biochemical investigation of UGT95A1, a family 1 GT enzyme from Pilosella officinarum that exhibits a strong, unusual regiospecificity for the 3'-O position of flavonoid acceptor substrate luteolin. We obtained an apo crystal structure to help drive the analyses of a series of binding site mutants, revealing that while most residues are tolerant to mutations, key residues M145 and D464 are important for overall glycosylation activity. Interestingly, E347 is crucial for maintaining the strong preference for 3'-O glycosylation, while R462 can be mutated to increase regioselectivity. The structural determinants of regioselectivity were further confirmed in homologous enzymes. Our study also suggests that the enzyme contains large, highly dynamic, disordered regions. We showed that while most disordered regions of the protein have little to no implication in catalysis, the disordered regions conserved among investigated homologs are important to both the overall efficiency and regiospecificity of the enzyme. This report represents a comprehensive in-depth analysis of a family 1 GT enzyme with a unique substrate regiospecificity and may provide a basis for enzyme functional prediction and engineering.

Defining the functional properties of cyclopropane fatty acid synthase from Pseudomonas aeruginosa PAO1.

Cyclopropane fatty acid synthases (CFAS) catalyze the conversion of unsaturated fatty acids to cyclopropane fatty acids (CFAs) within bacterial membranes. This modification alters the biophysical properties of membranes and has been correlated with virulence in several human pathogens. Despite the central role played by CFAS enzymes in regulating bacterial stress responses, the mechanistic properties of the CFAS enzyme family and the consequences of CFA biosynthesis remain largely uncharacterized in most bacteria. We report the first characterization of the CFAS enzyme from Pseudomonas aeruginosa (PA), an opportunistic human pathogen with complex membrane biology that is frequently associated with antimicrobial resistance and high tolerance to various external stressors. We demonstrate that CFAs are produced by a single enzyme in PA and that cfas gene expression is upregulated during the transition to stationary phase and in response to oxidative stress. Analysis of PA lipid extracts reveal a massive increase in CFA production as PA cells enter stationary phase and help define the optimal membrane composition for in vitro assays. The purified PA-CFAS enzyme forms a stable homodimer and preferentially modifies phosphatidylglycerol lipid substrates and membranes with a higher content of unsaturated acyl chains. Bioinformatic analysis across bacterial phyla shows highly divergent amino acid sequences within the lipid-binding domain of CFAS enzymes, perhaps suggesting distinct membrane-binding properties among different orthologs. This work lays an important foundation for further characterization of CFAS in P. aeruginosa and for examining the functional differences between CFAS enzymes from different bacteria.

Cdc14 phosphatases use an intramolecular pseudosubstrate motif to stimulate and regulate catalysis.

Cdc14 phosphatases are related structurally and mechanistically to protein tyrosine phosphatases (PTPs) but evolved a unique specificity for phosphoSer-Pro-X-Lys/Arg sites primarily deposited by cyclin-dependent kinases. This specialization is widely conserved in eukaryotes. The evolutionary reconfiguration of the Cdc14 active site to selectively accommodate phosphoSer-Pro likely required modification to the canonical PTP catalytic cycle. While studying Saccharomyces cerevisiae Cdc14, we discovered a short sequence in the disordered C terminus, distal to the catalytic domain, which mimics an optimal substrate. Kinetic analyses demonstrated this pseudosubstrate binds the active site and strongly stimulates rate-limiting phosphoenzyme hydrolysis, and we named it "substrate-like catalytic enhancer" (SLiCE). The SLiCE motif is found in all Dikarya fungal Cdc14 orthologs and contains an invariant glutamine, which we propose is positioned via substrate-like contacts to assist orientation of the hydrolytic water, similar to a conserved active site glutamine in other PTPs that Cdc14 lacks. AlphaFold2 predictions revealed vertebrate Cdc14 orthologs contain a conserved C-terminal alpha helix bound to the active site. Although apparently unrelated to the fungal sequence, this motif also makes substrate-like contacts and has an invariant glutamine in the catalytic pocket. Altering these residues in human Cdc14A and Cdc14B demonstrated that it functions by the same mechanism as the fungal motif. However, the fungal and vertebrate SLiCE motifs were not functionally interchangeable, illuminating potential active site differences during catalysis. Finally, we show that the fungal SLiCE motif is a target for phosphoregulation of Cdc14 activity. Our study uncovered evolution of an unusual stimulatory pseudosubstrate motif in Cdc14 phosphatases.

Catalytic mechanism of trans-2-enoyl-CoA reductases in the fatty acid elongation cycle and its cooperative action with fatty acid elongases.

The fatty acid (FA) elongation cycle produces very-long-chain FAs with ≥C21, which have unique physiological functions. Trans-2-enoyl-CoA reductases (yeast, Tsc13; mammals, TECR) catalyze the reduction reactions in the fourth step of the FA elongation cycle and in the sphingosine degradation pathway. However, their catalytic residues and coordinated action in the FA elongation cycle complex are unknown. To reveal these, we generated and analyzed Ala-substituted mutants of 15 residues of Tsc13. An in vitro FA elongation assay showed that nine of these mutants were less active than WT protein, with E91A and Y256A being the least active. Growth complementation analysis, measurement of ceramide levels, and deuterium-sphingosine labeling revealed that the function of the E91A mutant was substantially impaired in vivo. In addition, we found that the activity of FA elongases, which catalyze the first step of the FA elongation cycle, were reduced in the absence of Tsc13. Similar results were observed in Tsc13 E91A-expressing cells, which is attributable to reduced interaction between the Tsc13 E91A mutant and the FA elongases Elo2/Elo3. Finally, we found that E94A and Y248A mutants of human TECR, which correspond to E91A and Y256A mutants of Tsc13, showed reduced and almost no activity, respectively. Based on these results and the predicted three-dimensional structure of Tsc13, we speculate that Tyr256/Tyr248 of Tsc13/TECR is the catalytic residue that supplies a proton to trans-2-enoyl-CoAs. Our findings provide a clue concerning the catalytic mechanism of Tsc13/TECR and the coordinated action in the FA elongation cycle complex.

G-site residue S67 is involved in the fungicide-degrading activity of a tau class glutathione S-transferase from Carica papaya.

Thiram is a toxic fungicide extensively used for the management of pathogens in fruits. Although it is known that thiram degrades in plant tissues, the key enzymes involved in this process remain unexplored. In this study, we report that a tau class glutathione S-transferase (GST) from Carica papaya can degrade thiram. This enzyme was easily obtained by heterologous expression in Escherichia coli, showed low promiscuity toward other thiuram disulfides, and catalyzed thiram degradation under physiological reaction conditions. Site-directed mutagenesis indicated that G-site residue S67 shows a key influence for the enzymatic activity toward thiram, while mutation of residue S13, which reduced the GSH oxidase activity, did not significantly affect the thiram-degrading activity. The formation of dimethyl dithiocarbamate, which was subsequently converted into carbon disulfide, and dimethyl dithiocarbamoylsulfenic acid as the thiram degradation products suggested that thiram undergoes an alkaline hydrolysis that involves the rupture of the disulfide bond. Application of the GST selective inhibitor 4-chloro-7-nitro-2,1,3-benzoxadiazole reduced papaya peel thiram-degrading activity by 95%, indicating that this is the main degradation route of thiram in papaya. GST from Carica papaya also catalyzed the degradation of the fungicides chlorothalonil and thiabendazole, with residue S67 showing again a key influence for the enzymatic activity. These results fill an important knowledge gap in understanding the catalytic promiscuity of plant GSTs and reveal new insights into the fate and degradation products of thiram in fruits.

Differential molecular mechanisms of substrate recognition by selenium methyltransferases, INMT and TPMT, in selenium detoxification and excretion.

It is known that the recommended dietary allowance of selenium (Se) is dangerously close to its tolerable upper intake level. Se is detoxified and excreted in urine as trimethylselenonium ion (TMSe) when the amount ingested exceeds the nutritional level. Recently, we demonstrated that the production of TMSe requires two methyltransferases: thiopurine S-methyltransferase (TPMT) and indolethylamine N-methyltransferase (INMT). In this study, we investigated the substrate recognition mechanisms of INMT and TPMT in the Se-methylation reaction. Examination of the Se-methyltransferase activities of two paralogs of INMT, namely, nicotinamide N-methyltransferase and phenylethanolamine N-methyltransferase, revealed that only INMT exhibited Se-methyltransferase activity. Consistently, molecular dynamics simulations demonstrated that dimethylselenide was preferentially associated with the active center of INMT. Using the fragment molecular orbital method, we identified hydrophobic residues involved in the binding of dimethylselenide to the active center of INMT. The INMT-L164R mutation resulted in a deficiency in Se- and N-methyltransferase activities. Similarly, TPMT-R152, which occupies the same position as INMT-L164, played a crucial role in the Se-methyltransferase activity of TPMT. Our findings suggest that TPMT recognizes negatively charged substrates, whereas INMT recognizes electrically neutral substrates in the hydrophobic active center embedded within the protein. These observations explain the sequential requirement of the two methyltransferases in producing TMSe.