Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Docking”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Tetrazole-based inhibitors of the bacterial enzyme N -succinyl- L,L -2,6-diaminopimelic acid desuccinylase as potential antibiotics

Based on a hit from a high-throughput screen, a series of phenyltetrazole amides was synthesized and assayed for inhibitory potency against DapE from Haemophilus influenzae (HiDapE). The inhibitory potency was modest but confirmed, with the most potent analog containing an aminothiazole moiety displaying an IC 50 = 50.2 ± 5.0 μM. Finally, docking reveals a potential binding mode wherein the amide carbonyl bridges both zinc atoms in the active site, and the tetrazole forms key hydrogen bonds with Arg330.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. In this study, we propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Rational Enzyme Design for Controlled Functionalization of Acetylated Xylan for Cell-Free Polymer Biosynthesis

Xylan O-acetyltransferase 1 (XOAT1) is involved in O-acetylating the backbone of hemicellulose xylan. Recent structural analysis of XOAT1 showed two unequal lobes forming a cleft that is predicted to accommodate and position xylan acceptors into proximity with the catalytic triad. Here, we used docking and molecular dynamics simulations to investigate the optimal orientation of xylan in the binding cleft of XOAT1 and identify putative key residues (Gln445 and Arg444 on Minor lobe & Asn312, Met311 and Asp403 on Major lobe) involved in substrate interactions. Site-directed mutagenesis coupled with biochemical analyses revealed the major lobe of XOAT1 is important for xylan binding. Mutation of single key residues yielded XOAT1 variants with various enzymatic efficiencies that are applicable to one-pot synthesis of xylan polymers with different degrees of O-acetylation. Taken together, our results demonstrate the effectiveness of computational modeling in guiding enzyme engineering aimed at modulating xylan and redesigning plant cell walls.

93 SC - Biological and Environmental Research (BER↗
DOE OSTI2023

Identification of inhibitors against SARS-CoV-2 variants of concern using virtual screening and metadynamics-based enhanced sampling

Among the variants of SARS-CoV-2, some are more infectious than the Wild-type. Interestingly, these mutations enable the virus to evade the therapeutic efforts. Hence, there is a need for candidate drug molecules that can potently bind with all the variants. Here we have adopted a strategy combining virtual screening, molecular docking followed by rigorous sampling by metadynamics simulations to find candidate molecules. From our results we found four highly potent drug candidates that can bind to the Spike-RBD of all the variants of the virus. Additionally, we also found that certain signature residues on the RBM region commonly bind to each of these inhibitors. Thus, our study not only gives information on the chemical compounds, but also residues on the proteins which could be targeted for future drug and vaccine development studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Engineering a photoenzyme to use red light

Photoenzymatic reactions involving flavin-dependent “ene”-reductases (EREDs) rely on protein-templated charge transfer (CT) complexes between the cofactor and substrate for radical initiation. These complexes typically absorb in the blue region of the electromagnetic spectrum. Here, we engineered an ERED to form CT complexes that absorb red light. Mechanistic studies indicate that red-light activity is due to the growth of a red-absorbing shoulder off the previously identified cyan absorption feature. Molecular dynamics simulations, docking, and excited-state calculations suggest that the cyan feature involves a π→π* transition on flavin, whereas the red-light absorption is a π→π* transition between flavin and the substrate. Furthermore, differences in the electronic transition are due to changes in the substrate-binding conformation and allosteric tuning of the electronic structure of the cofactor-substrate complex. Microenvironment tuning of the CT complex for red-light activity is observed with other engineered photoenzymatic reactions, highlighting this effect’s generality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Multisource Rotating Shield Brachytherapy Apparatus for Prostate Cancer

Purpose: Our purpose is to present a novel multisource rotating shield brachytherapy (RSBT) apparatus for the simultaneous precise angular and linear positioning of partially shielded {sup 153}Gd brachytherapy sources in interstitial needles for the treatment of locally advanced prostate cancer. It is designed to lower the dose to nearby healthy tissues, the urethra in particular, relative to conventional high-dose-rate brachytherapy techniques. Methods and Materials: Following needle implantation through the patient template, an angular drive mechanism is docked to the patient template. Each needle is coupled to a multisource afterloader catheter by a connector passing through a shaft. The shafts are rotated about their axes by translating a moving template between 2 stationary templates. The shafts' surfaces and moving template holes are helically threaded with the same pattern such that translation of the moving template causes simultaneous rotation of the shafts. The rotation of each shaft is mechanically transmitted to the catheter-source-shield combination, inside the needles, via several key-keyway pairs. The catheter angles are simultaneously incremented throughout treatment, and only a single 360° rotation of all catheters is needed for a full treatment. For each rotation angle, source depth in each needle is controlled by a multisource afterloader, which is proposed as an array of belt-driven linear actuators, each of which drives a wire that controls catheter depth in a needle. Results: Treatment plans demonstrated that RSBT with the proposed apparatus reduced urethral D{sub 0.1cm{sup 3}} (the minimum dose delivered to the hottest 0.1cm{sup 3} of the urethra) below that of conventional high-dose-rate brachytherapy by 31% for urethral dose gradient volume within 3 mm of the urethra surface. Treatment time to deliver 20 Gy with the proposed multisource RSBT apparatus by use of nineteen 62.4-GBq {sup 153}Gd sources was 122 minutes. Conclusions: The proposed RSBT delivery apparatus enables a mechanically feasible urethra-sparing treatment technique for prostate cancer in a clinically reasonable time frame.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

Engineered Biosynthesis of Alkyne-Tagged Polyketides by Type I PKSs

Polyketides produced by modular polyketide synthases (PKSs) are important small molecules widely used as drugs, pesticides, and biological probes. Tagging these polyketides with a clickable functionality enables the visualization, diversification, and mode of action study through bio-orthogonal chemistry. We report the de novo biosynthesis of alkyne-tagged polyketides by modular type I PKSs through starter unit engineering. Specifically, we use JamABC, a terminal alkyne biosynthetic machinery from the jamaicamide B biosynthetic pathway, in combination with representative modular PKSs. We demonstrate that JamABC works as a trans loading system for engineered type I PKSs to produce alkyne-tagged polyketides. In addition, the production efficiency can be improved by enhancing the interactions between the carrier protein (JamC) and PKSs using docking domains and site-directed mutagenesis of JamC. This work thus provides engineering guidelines and strategies that are applicable to additional modular type I PKSs to produce targeted alkyne-tagged metabolites for chemical and biological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Structure of 3-mercaptopropionic acid dioxygenase with a substrate analog reveals bidentate substrate binding at the iron center

Thiol dioxygenases are a subset of nonheme iron oxygenases that catalyze the formation of sulfinic acids from sulfhydryl-containing substrates and dioxygen. Among this class, cysteine dioxygenases (CDOs) and 3-mercaptopropionic acid dioxygenases (3MDOs) are the best characterized, and the mode of substrate binding for CDOs is well understood. However, the manner in which 3-mercaptopropionic acid (3MPA) coordinates to the nonheme iron site in 3MDO remains a matter of debate. A model for bidentate 3MPA coordination at the 3MDO Fe-site has been proposed on the basis of computational docking, whereas steady-state kinetics and EPR spectroscopic measurements suggest a thiolate-only coordination of the substrate. To address this gap in knowledge, we determined the structure of Azobacter vinelandii 3MDO (Av3MDO) in complex with the substrate analog and competitive inhibitor, 3-hydroxypropionic acid (3HPA). The structure together with DFT computational modeling demonstrates that 3HPA and 3MPA associate with iron as chelate complexes with the substrate-carboxylate group forming an additional interaction with Arg168 and the thiol bound at the same position as in CDO. A chloride ligand was bound to iron in the coordination site assigned as the O 2 -binding site. Supporting HYSCORE spectroscopic experiments were performed on the (3MPA/NO)-bound Av3MDO iron nitrosyl (S = 3/2) site. In combination with spectroscopic simulations and optimized DFT models, this work provides an experimentally verified model of the Av3MDO enzyme–substrate complex, effectively resolving a debate in the literature regarding the preferred substrate-binding denticity. These results elegantly explain the observed 3MDO substrate specificity, but leave unanswered questions regarding the mechanism of substrate-gated reactivity with dioxygen.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular chaperone RAP interacts with LRP1 in a dynamic bivalent mode and enhances folding of ligand-binding regions of other LDLR family receptors

The low-density lipoprotein receptor (LDLR) family of receptors are cell-surface receptors that internalize numerous ligands and play crucial role in various processes, such as lipoprotein metabolism, hemostasis, fetal development, etc. Previously, receptor-associated protein (RAP) was described as a molecular chaperone for LDLR-related protein 1 (LRP1), a prominent member of the LDLR family. We aimed to verify this role of RAP for LRP1 and two other LDLR family receptors, LDLR and vLDLR, and to investigate the mechanisms of respective interactions using a cell culture model system, purified system, and in silico modelling. Upon coexpression of RAP with clusters of the ligand-binding complement repeats (CRs) of the receptors in secreted form in insect cells culture, the isolated proteins had increased yield, enhanced folding, and improved binding properties compared with proteins expressed without RAP, as determined by circular dichroism and surface plasmon resonance. Within LRP1 CR-clusters II and IV, we identified multiple sites comprised of adjacent CR doublets, which provide alternative bivalent binding combinations with specific pairs of lysines on RAP. Mutational analysis of these lysines within each of isolated RAP D1/D2 and D3 domains having high affinity to LRP1 and of conserved tryptophans on selected CR-doublets of LRP1, as well as in silico docking of a model LRP1 CR-triplet with RAP, indicated a universal role for these residues in interaction of RAP and LRP1. Consequently, we propose a new model of RAP interaction with LDLR family receptors based on switching of the bivalent contacts between molecules over time in a dynamic mode.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Identification of structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Characterization of lignin-degrading enzyme PmdC, which catalyzes a key step in the synthesis of polymer precursor 2-pyrone-4,6-dicarboxylic acid

Pyrone-2,4-dicarboxylic acid (PDC) is a valuable polymer precursor that can be derived from the microbial degradation of lignin. The key enzyme in the microbial production of PDC is 4-carboxy-2-hydroxymuconate-6-semialdehyde (CHMS) dehydrogenase, which acts on the substrate CHMS. We present the crystal structure of CHMS dehydrogenase (PmdC from Comamonas testosteroni) bound to the cofactor NADP, shedding light on its three-dimensional architecture, and revealing residues responsible for binding NADP. Using a combination of structural homology, molecular docking, and quantum chemistry calculations, we have predicted the binding site of CHMS. Key histidine residues in a conserved sequence are identified as crucial for binding the hydroxyl group of CHMS and facilitating dehydrogenation with NADP. Mutating these histidine residues results in a loss of enzyme activity, leading to a proposed model for the enzyme's mechanism. These findings are expected to help guide efforts in protein and metabolic engineering to enhance PDC yields in biological routes to polymer feedstock synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The crystal structure of Grindelia robusta 7,13-copalyl diphosphate synthase reveals active site features controlling catalytic specificity

Diterpenoid natural products serve critical functions in plant development and ecological adaptation and many diterpenoids have economic value as bioproducts. The family of class II diterpene synthases catalyzes the committed reactions in diterpenoid biosynthesis, converting a common geranylgeranyl diphosphate precursor into different bicyclic prenyl diphosphate scaffolds. Enzymatic rearrangement and modification of these precursors generate the diversity of bioactive diterpenoids. We report the crystal structure of Grindelia robusta 7,13-copalyl diphosphate synthase, GrTPS2, at 2.1 Å of resolution. GrTPS2 catalyzes the committed reaction in the biosynthesis of grindelic acid, which represents the signature metabolite in species of gumweed (Grindelia spp., Asteraceae). Grindelic acid has been explored as a potential source for drug leads and biofuel production. The GrTPS2 crystal structure adopts the conserved three-domain fold of class II diterpene synthases featuring a functional active site in the γβ-domain and a vestigial α-domain. Substrate docking into the active site of the GrTPS2 apo protein structure predicted catalytic amino acids. Biochemical characterization of protein variants identified residues with impact on enzyme activity and catalytic specificity. Specifically, mutagenesis of Y457 provided mechanistic insight into the position-specific deprotonation of the intermediary carbocation to form the characteristic 7,13 double bond of 7,13-copalyl diphosphate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A characterization of recombinant Arabidopsis FRIABLE1 (FRB1) reveals robust rhamnogalacturonan-I rhamnosyltransferase activity and critical catalytic residues

Plant cell walls are glycan-rich extracellular matrices that fundamentally impact essential cellular processes, such as growth, adhesion, and cell shape acquisition. Understanding plant cell wall glycans requires the identification and characterization of the biosynthetic enzymes that produce these polymers. Most successful in vitro protein expression studies of plant cell wall glycosyltransferases have relied on insect, fungal/yeast, or human cell expression systems, whereas prokaryotic expression systems have been generally unsuccessful. Here, we show that Arabidopsis FRIABLE1 (FRB1)/rhamnogalacturonan-I rhamnosyltransferase 8 (RRT8) can be produced in Escherichia coli RosettaGami2 cells as N-terminal maltose-binding protein fusion proteins containing C-terminal 6X-His-tags. We also report the catalytic constants of FRB1/RRT8 with apparent K M and K cat values of 226 μM and 33 min -1 for UDP-Rhamnose and 117 μM and 28.7 min -1 for rhamnogalacturonan-I (RG-I), respectively. We examine the catalytic activities of mutated FRB1/RRT8 proteins based on an AlphaFold 3-generated FRB1/RRT8 protein structural model with a virtually docked UDP-Rha donor. Enzymatic characterization of the mutated and wildtype FRB1/RRT8 protein confirmed that mutation of predicted catalytic site amino acid residues resulted in a 20-fold reduction in RRT activity. FRB1 also robustly polymerizes RG-I in combination with RG-I galacturonosyltransferase 1. These results show how a robust E. coli expression system combined with artificial intelligence tools can be used to increase understanding of plant cell wall glycosyltransferase structure and function.

glycosyltransferase↗
DOE OSTI2022

BIPSPI+: Mining Type-Specific Datasets of Protein Complexes to Improve Protein Binding Site Prediction

Computational approaches for predicting protein-protein interfaces are extremely useful for understanding and modelling the quaternary structure of protein assemblies. In particular, partner-specific binding site prediction methods allow delineating the specific residues that compose the interface of protein complexes. In recent years, new machine learning and other algorithmic approaches have been proposed to solve this problem. However, little effort has been made in finding better training datasets to improve the performance of these methods. With the aim of vindicating the importance of the training set compilation procedure, in this work we present BIPSPI+, a new version of our original server trained on carefully curated datasets that outperforms our original predictor. We show how prediction performance can be improved by selecting specific datasets that better describe particular types of protein interactions and interfaces (e.g. homo/hetero). In addition, our upgraded web server offers a new set of functionalities such as the sequence-structure prediction mode, hetero- or homo-complex specialization and the guided docking tool that allows to compute 3D quaternary structure poses using the predicted interfaces. BIPSPI+ is freely available at https://bipspi.cnb.csic.es.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

AlphaFold and Structural Mass Spectrometry Enable Interrogations on the Intrinsically Disordered Regions in Cyanobacterial Light-harvesting Complex Phycobilisome

Intrinsically disordered proteins/regions (IDPRs) are a very large and functionally important class of proteins that participate in weak multivalent interactions in protein complexes. They are recalcitrant for interrogations using X-ray crystallography and cryo-EM. The IDPRs observed at the interface of the photosynthetic pigment protein complexes (PPCs) remain much less clear, e.g., the major cyanobacterial light-harvesting complex (PBS) contains an unstructured PB-loop insertion in the phycocyanobilin domain (PB domain) of ApcE (the largest polypeptide in PBS). Here, a joint platform is built to probe such structural domains. This platform is characterized by two-round progressive justifications of in silico models by using the structural mass spectrometry data. First, the AlphaFold-generated 3D structure of the PB domain (containing PB-loop) was justified in the context of PBS. Second, docking the AlphaFold-generated ApcG (a ligand) into the first-step justified structure (a receptor). The final ligand-receptor complex was then subjected to a second-round justification, again, by using unequivocal isotopically-encoded cross-links identified in LC-MS/MS. This work reveals a full-length PB-loop structure modelled in the PBS basal cylinder, free from any spatial conflicts against the other subunits in PBS. The structure of PB domain highlights the close associations of the intrinsically disordered PB-loop with its binding partners in PBS, including ApcG, another IDPR. The PB-loop region involved in the binding of photosystem II (PSII) is also discussed in the context of excitation energy transfer regulation. Finally, this work calls attention to the highly disordered, yet interrogatable interface between the light-harvesting antenna complexes and the reaction centers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Efficient discrimination of natural stereoisomers of chicoric acid, an HIV-1 integrase inhibitor

Plants from the Asteraceae family are known to contain a wide spectrum of phytochemicals with various nutraceutical properties. One important phytochemical, chicoric acid (CA), is reported to exist in plants, such as Sonchus oleraceus and Bidens pilosa, as stereoisomers. They occur either as the naturally abundant RR-chicoric acid (RR-CA), or the less abundant RS-chicoric acid (RS-CA), also known as meso-chicoric acid. To date, little is known about the biological activity of RS-CA, but there is evidence of its anti-HIV (human immunodeficiency virus) properties. In this study, a reliable analytical method was developed to distinguish between the two stereoisomers detected in S. oleraceus and B. pilosa. Structure identification used liquid chromatography–mass spectrometry (LC-MS) in combination with ultraviolet radiation (UV) -induced geometrical isomerization, molecular dynamics (MD) simulations, and density functional theory (DFT) models. Optimized structures from DFT calculations were used for docking studies against the HIV-1 integrase enzyme. Different retention times on the reverse phase chromatograms revealed the plants produce two different CA stereoisomers: S. oleraceus produced the RR-CA isomer, while B. pilosa produced the RS-CA isomer. DFT results demonstrated the RR-CA molecule was more stable than RS-CA due to the stabilizing force of intra-molecular hydrogen bonding. Differences in the HIV-1 integrase enzyme binding modes were observed, with the RR-CA being a more potent inhibitor than the RS-CA molecule. The results highlight the significance of plant metabolite structural complexity from both chemical and biological perspectives. Artificial formation of geometrical isomers, in combination with the predictive ability of DFT models and the resolving power of the LC-MS, can be exploited to distinguish closely related compounds, in particular stereoisomers.

chicoric acid, RR-chicoric acid, RS-chicoric acid,↗
DOE OSTI2021

An effective human uracil-DNA glycosylase inhibitor targets the open pre-catalytic active site conformation

Human uracil DNA-glycosylase (UDG) is the prototypic and first identified DNA glycosylase with a vital role in removing deaminated cytosine and incorporated uracil and 5-fluorouracil (5-FU) from DNA. UDG depletion sensitizes cells to high APOBEC3B deaminase and to pemetrexed (PEM) and floxuridine (5-FdU), which are toxic to tumor cells through incorporation of uracil and 5-FU into DNA. To identify small-molecule UDG inhibitors for pre-clinical evaluation, we optimized biochemical screening of a selected diversity collection of >3,000 small-molecules. We found aurintricarboxylic acid (ATA) as an inhibitor of purified UDG at an initial calculated IC 50 < 100 nM. Subsequent enzymatic assays confirmed effective ATA inhibition but with an IC50 of 700 nM and showed direct binding to the human UDG with a K D of <700 nM. ATA displays preferential, dose-dependent binding to purified human UDG compared to human 8-oxoguanine DNA glycosylase. ATA did not bind uracil-containing DNA at these concentrations. Yet, combined crystal structure and in silico docking results unveil ATA interactions with the DNA binding channel and uracil-binding pocket in an open, destabilized UDG conformation. Biologically relevant ATA inhibition of UDG was measured in cell lysates from human DLD1 colon cancer cells and in MCF-7 breast cancer cells using a host cell reactivation assay. Collective findings provide proof-of-principle for development of an ATA-based chemotype and “door stopper” strategy targeting inhibitor binding to a destabilized, open pre-catalytic glycosylase conformation that prevents active site closing for functional DNA binding and nucleotide flipping needed to excise altered bases in DNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Selectivity of enzymes involved in the formation of opposite enantiomeric series of p -menthane monoterpenoids in peppermint and Japanese catnip

Peppermint (Mentha x piperita L.) and Japanese catnip (Schizonepeta tenuifolia (Benth.) Briq.) accumulate p-menthane monoterpenoids with identical functionalization patterns but opposite stereochemistry. In the present study, we investigate the enantioselectivity of multiple enzymes involved in monoterpenoid biosynthesis in these species. Based on kinetic assays, mint limonene synthase, limonene 3-hydroxylase, isopiperitenol dehydrogenase, isopiperitenone reductase, and menthone reductase exhibited significant enantioselectivity toward intermediates of the pathway that proceeds through (-)-4S-limonene. Limonene synthase, isopiperitenol dehydrogenase and isopiperitenone reductase of Japanese catnip preferred intermediates of the pathway that involves (+)-4R-limonene, whereas limonene 3-hydroxylase was not enantioselective, and the activities of pulegone reductase and menthone reductase were too low to acquire meaningful kinetic data. Molecular modeling studies with docked ligands generally supported the experimental data obtained with peppermint enzymes, indicating that the preferred enantiomer was aligned well with the requisite cofactor and amino acid residues implicated in catalysis. A striking example for enantioselectivity was peppermint (-)-menthone reductase, which binds (-)-menthone with exquisite affinity but was predicted to bind (+)-menthone in a non-productive orientation that positions its carbonyl functional group at considerable distance to the NADPH cofactor. Here, the work presented here lays the groundwork for structure-function studies aimed at unraveling how enantioselectivity evolved in closely related species of the Lamiaceae and beyond.

59 BASIC BIOLOGICAL SCIENCES↗