Boron catalysis in a designer enzyme.

Enzymes play an increasingly important role in improving the benignity and efficiency of chemical production, yet the diversity of their applications lags heavily behind chemical catalysts as a result of the relatively narrow range of reaction mechanisms of enzymes. The creation of enzymes containing non-biological functionalities facilitates reaction mechanisms outside nature's canon and paves the way towards fully programmable biocatalysis1-3. Here we present a completely genetically encoded boronic-acid-containing designer enzyme with organocatalytic reactivity not achievable with natural or engineered biocatalysts4,5. This boron enzyme catalyses the kinetic resolution of hydroxyketones by oxime formation, in which crucial interactions with the protein scaffold assist in the catalysis. A directed evolution campaign led to a variant with natural-enzyme-like enantioselectivities for several different substrates. The unique activation mode of the boron enzyme was confirmed using X-ray crystallography, high-resolution mass spectrometry (HRMS) and 11B NMR spectroscopy. Our study demonstrates that genetic-code expansion can be used to create evolvable enantioselective enzymes that rely on xenobiotic catalytic moieties such as boronic acids and access reaction mechanisms not reachable through catalytic promiscuity of natural or engineered enzymes.

Structural Elucidation and Engineering of a Bacterial Carbohydrate Oxidase.

Flavin-dependent carbohydrate oxidases are valuable tools in biotechnological applications due to their high selectivity in the oxidation of carbohydrates. In this study, we report the biochemical and structural characterization of a recently discovered carbohydrate oxidase from the bacterium Ralstonia solanacearum, which is a member of the vanillyl alcohol oxidase flavoprotein family. Due to its exceptionally high activity toward N-acetyl-d-galactosamine and N-acetyl-d-glucosamine, the enzyme was named N-acetyl-glucosamine oxidase (NagOx). In contrast to most known (fungal) carbohydrate oxidases, NagOx could be overexpressed in a bacterial host, which facilitated detailed biochemical and enzyme engineering studies. Steady state kinetic analyses revealed that non-acetylated hexoses were also accepted as substrates albeit with lower efficiency. Upon determination of the crystal structure, structural insights into NagOx were obtained. A large cavity containing a bicovalently bound FAD, tethered via histidyl and cysteinyl linkages, was observed. Substrate docking highlighted how a single residue (Leu251) plays a key role in the accommodation of N-acetylated sugars in the active site. Upon replacement of Leu251 (L251R mutant), an enzyme variant was generated with a drastically modified substrate acceptance profile, tuned toward non-N-acetylated monosaccharides and disaccharides. Furthermore, the activity toward bulkier substrates such as the trisaccharide maltotriose was introduced by this mutation. Due to its advantage of being overexpressed in a bacterial host, NagOx can be considered a promising alternative engineerable biocatalyst for selective oxidation of monosaccharides and oligosaccharides.

Fixing Flavins: Hijacking a Flavin Transferase for Equipping Flavoproteins with a Covalent Flavin Cofactor.

Most flavin-dependent enzymes contain a dissociable flavin cofactor. We present a new approach for installing in vivo a covalent bond between a flavin cofactor and its host protein. By using a flavin transferase and carving a flavinylation motif in target proteins, we demonstrate that "dissociable" flavoproteins can be turned into covalent flavoproteins. Specifically, four different flavin mononucleotide-containing proteins were engineered to undergo covalent flavinylation: a light-oxygen-voltage domain protein, a mini singlet oxygen generator, a nitroreductase, and an old yellow enzyme-type ene reductase. Optimizing the flavinylation motif and expression conditions led to the covalent flavinylation of all four flavoproteins. The engineered covalent flavoproteins retained function and often exhibited improved performance, such as higher thermostability or catalytic performance. The crystal structures of the designed covalent flavoproteins confirmed the designed threonyl-phosphate linkage. The targeted flavoproteins differ in fold and function, indicating that this method of introducing a covalent flavin-protein bond is a powerful new method to create flavoproteins that cannot lose their cofactor, boosting their performance.

Selecting Better Biocatalysts by Complementing Recoded Bacteria.

In vivo selections are powerful tools for the directed evolution of enzymes. However, the need to link enzymatic activity to cellular survival makes selections for enzymes that do not fulfill a metabolic function challenging. Here, we present an in vivo selection strategy that leverages recoded organisms addicted to non-canonical amino acids (ncAAs) to evolve biocatalysts that can provide these building blocks from synthetic precursors. We exemplify our platform by engineering carbamoylases that display catalytic efficiencies more than five orders of magnitude higher than those observed for the wild-type enzyme for ncAA-precursors. As growth rates of bacteria under selective conditions correlate with enzymatic activities, we were able to elicit improved variants from populations by performing serial passaging. By requiring minimal human intervention and no specialized equipment, we surmise that our strategy will become a versatile tool for the in vivo directed evolution of diverse biocatalysts.

Gene Fusion and Directed Evolution to Break Structural Symmetry and Boost Catalysis by an Oligomeric C-C Bond-Forming Enzyme.

Gene duplication and fusion are among the primary natural processes that generate new proteins from simpler ancestors. Here we adopted this strategy to evolve a promiscuous homohexameric 4-oxalocrotonate tautomerase (4-OT) into an efficient biocatalyst for enantioselective Michael reactions. We first designed a tandem-fused 4-OT to allow independent sequence diversification of adjacent subunits by directed evolution. This fused 4-OT was then subjected to eleven rounds of directed evolution to give variant 4-OT(F11), which showed an up to 320-fold enhanced activity for the Michael addition of nitromethane to cinnamaldehydes. Crystallographic analysis revealed that 4-OT(F11) has an unusual asymmetric trimeric architecture in which one of the monomers is flipped 180° relative to the others. This gene duplication and fusion strategy to break structural symmetry is likely to become an indispensable asset of the enzyme engineering toolbox, finding wide use in engineering oligomeric proteins.

Catalytic and structural properties of ATP-dependent caprolactamase from Pseudomonas jessenii.

Caprolactamase is the first enzyme in the caprolactam degradation pathway of Pseudomonas jessenii. It is composed of two subunits (CapA and CapB) and sequence-related to other ATP-dependent enzymes involved in lactam hydrolysis, like 5-oxoprolinases and hydantoinases. Low sequence similarity also exists with ATP-dependent acetone- and acetophenone carboxylases. The caprolactamase was produced in Escherichia coli, isolated by His-tag affinity chromatography, and subjected to functional and structural studies. Activity toward caprolactam required ATP and was dependent on the presence of bicarbonate in the assay buffer. The hydrolysis product was identified as 6-aminocaproic acid. Quantum mechanical modeling indicated that the hydrolysis of caprolactam was highly disfavored (ΔG0 '= 23 kJ/mol), which explained the ATP dependence. A crystal structure showed that the enzyme exists as an (αß)2 tetramer and revealed an ATP-binding site in CapA and a Zn-coordinating site in CapB. Mutations in the ATP-binding site of CapA (D11A and D295A) significantly reduced product formation. Mutants with substitutions in the metal binding site of CapB (D41A, H99A, D101A, and H124A) were inactive and less thermostable than the wild-type enzyme. These residues proved to be essential for activity and on basis of the experimental findings we propose possible mechanisms for ATP-dependent lactam hydrolysis.

Structure of a robust bacterial protein cage and its application as a versatile biocatalytic platform through enzyme encapsulation.

Using a newly discovered encapsulin from Mycolicibacterium hassiacum, several biocatalysts were packaged in this robust protein cage. The encapsulin was found to be easy to produce as recombinant protein. Elucidation of its crystal structure revealed that it is a spherical protein cage of 60 protomers (diameter of 23 nm) with narrow pores. By developing an effective coexpression and isolation procedure, the effect of packaging a variety of biocatalysts could be evaluated. It was shown that encapsulation results in a significantly higher stability of the biocatalysts. Most of the targeted cofactor-containing biocatalysts remained active in the encapsulin. Due to the restricted diameters of the encapsulin pores (5-9 Å), the protein cage protects the encapsulated enzymes from bulky compounds. The work shows that encapsulins may be valuable tools to tune the properties of biocatalysts such as stability and substrate specificity.

Characterization of a New DyP-Peroxidase from the Alkaliphilic Cellulomonad, Cellulomonas bogoriensis.

DyP-type peroxidases are heme-containing enzymes that have received increasing attention over recent years with regards to their potential as biocatalysts. A novel DyP-type peroxidase (CboDyP) was discovered from the alkaliphilic cellulomonad, Cellulomonas bogoriensis, which could be overexpressed in Escherichia coli. The biochemical characterization of the recombinant enzyme showed that it is a heme-containing enzyme capable to act as a peroxidase on several dyes. With the tested substrates, the enzyme is most active at acidic pH values and is quite tolerant towards solvents. The crystal structure of CboDyP was solved which revealed atomic details of the dimeric heme-containing enzyme. A peculiar feature of CboDyP is the presence of a glutamate in the active site which in most other DyPs is an aspartate, being part of the DyP-typifying sequence motif GXXDG. The E201D CboDyP mutant was prepared and analyzed which revealed that the mutant enzyme shows a significantly higher activity on several dyes when compared with the wild-type enzyme.

Design of a substrate-tailored peptiligase variant for the efficient synthesis of thymosin-α1.

The synthesis of thymosin-α1, an acetylated 28 amino acid long therapeutic peptide, via conventional chemical methods is exceptionally challenging. The enzymatic coupling of unprotected peptide segments in water offers great potential for a more efficient synthesis of peptides that are difficult to synthesize. Based on the design of a highly engineered peptide ligase, we developed a fully convergent chemo-enzymatic peptide synthesis (CEPS) process for the production of thymosin-α1via a 14-mer + 14-mer segment condensation strategy. Using structure-inspired enzyme engineering, the thiol-subtilisin variant peptiligase was tailored to recognize the respective 14-mer thymosin-α1 segments in order to create a clearly improved biocatalyst, termed thymoligase. Thymoligase catalyzes peptide bond formation between both segments with a very high efficiency (>94% yield) and is expected to be well applicable to many other ligations in which residues with similar characteristics (e.g. Arg and Glu) are present in the respective positions P1 and P1'. The crystal structure of thymoligase was determined and shown to be in good agreement with the model used for the engineering studies. The combination of the solid phase peptide synthesis (SPPS) of the 14-mer segments and their thymoligase-catalyzed ligation on a gram scale resulted in a significantly increased, two-fold higher overall yield (55%) of thymosin-α1 compared to those typical of existing industrial processes.

Characterization of Two VAO-Type Flavoprotein Oxidases from Myceliophthora thermophila.

The VAO flavoprotein family consists mostly of oxidoreductases harboring a covalently linked flavin cofactor. The linkage can be either monocovalent at position 8 with a histidine or tyrosine or bicovalent at position 8 with a histidine and at position 6 with a cysteine. Bicovalently bound flavoproteins show a preference for bulkier substrates such as oligosaccharides or secondary metabolites. The genome of the thermophilic fungus Myceliophthora thermophila C1 was found to be rich in genes encoding putative covalent VAO-type flavoproteins. Enzymes from this fungus have the advantage of being rather thermostable and homologous overexpression in M. thermophila C1 is feasible. Recently we discovered a new and VAO-type carbohydrate oxidase from this fungus: xylooligosaccharide oxidase. In this study, two other putative VAO-type oxidases, protein sequence XP_003663615 (MtVAO615) and XP_003665713 (MtVAO713), were expressed in M. thermophila C1, purified and characterized. Enzyme MtVAO615 was found to contain a bicovalently bound FAD, while enzyme MtVAO713 contained a monocovalent histidyl-bound FAD. The crystal structures of both proteins were obtained which revealed atypical active site architectures. It could be experimentally verified that both proteins, when reduced, rapidly react with molecular oxygen, a hallmark of flavoprotein oxidases. A large panel of alcohols, including carbohydrates, steroids and secondary alcohols were tested as potential substrates. For enzyme MtVAO713 low oxidase activity was discovered towards ricinoleic acid.

Metal Dependence of the Xylose Isomerase from Piromyces sp. E2 Explored by Activity Profiling and Protein Crystallography.

Xylose isomerase from Piromyces sp. E2 (PirXI) can be used to equip Saccharomyces cerevisiae with the capacity to ferment xylose to ethanol. The biochemical properties and structure of the enzyme have not been described even though its metal content, catalytic parameters, and expression level are critical for rapid xylose utilization. We have isolated the enzyme after high-level expression in Escherichia coli, analyzed the metal dependence of its catalytic properties, and determined 12 crystal structures in the presence of different metals, substrates, and substrate analogues. The activity assays revealed that various bivalent metals can activate PirXI for xylose isomerization. Among these metals, Mn2+ is the most favorable for catalytic activity. Furthermore, the enzyme shows the highest affinity for Mn2+, which was established by measuring the activation constants (Kact) for different metals. Metal analysis of the purified enzyme showed that in vivo the enzyme binds a mixture of metals that is determined by metal availability as well as affinity, indicating that the native metal composition can influence activity. The crystal structures show the presence of an active site similar to that of other xylose isomerases, with a d-xylose binding site containing two tryptophans and a catalytic histidine, as well as two metal binding sites that are formed by carboxylate groups of conserved aspartates and glutamates. The binding positions and conformations of the metal-coordinating residues varied slightly for different metals, which is hypothesized to contribute to the observed metal dependence of the isomerase activity.

Discovery of a Xylooligosaccharide Oxidase from Myceliophthora thermophila C1.

By inspection of the predicted proteome of the fungus Myceliophthora thermophila C1 for vanillyl-alcohol oxidase (VAO)-type flavoprotein oxidases, a putative oligosaccharide oxidase was identified. By homologous expression and subsequent purification, the respective protein could be obtained. The protein was found to contain a bicovalently bound FAD cofactor. By screening a large number of carbohydrates, several mono- and oligosaccharides could be identified as substrates. The enzyme exhibits a strong substrate preference toward xylooligosaccharides; hence it is named xylooligosaccharide oxidase (XylO). Chemical analyses of the product formed upon oxidation of xylobiose revealed that the oxidation occurs at C1, yielding xylobionate as product. By elucidation of several XylO crystal structures (in complex with a substrate mimic, xylose, and xylobiose), the residues that tune the unique substrate specificity and regioselectivity could be identified. The discovery of this novel oligosaccharide oxidase reveals that the VAO-type flavoprotein family harbors oxidases tuned for specific oligosaccharides. The unique substrate profile of XylO hints at a role in the degradation of xylan-derived oligosaccharides by the fungus M. thermophila C1.

Functional and structural characterization of an unusual cofactor-independent oxygenase.

The vast majority of characterized oxygenases use bound cofactors to activate molecular oxygen to carry out oxidation chemistry. Here, we show that an enzyme of unknown activity, RhCC from Rhodococcus jostii RHA1, functions as an oxygenase, using 4-hydroxyphenylenolpyruvate as a substrate. This unique and complex reaction yields 3-hydroxy-3-(4-hydroxyphenyl)-pyruvate, 4-hydroxybenzaldehyde, and oxalic acid as major products. Incubations with H2(18)O, (18)O2, and a substrate analogue suggest that this enzymatic oxygenation reaction likely involves a peroxide anion intermediate. Analysis of sequence similarity and the crystal structure of RhCC (solved at 1.78 Å resolution) reveal that this enzyme belongs to the tautomerase superfamily. Members of this superfamily typically catalyze tautomerization, dehalogenation, or decarboxylation reactions rather than oxygenation reactions. The structure shows the absence of cofactors, establishing RhCC as a rare example of a redox-metal- and coenzyme-free oxygenase. This sets the stage to study the mechanistic details of cofactor-independent oxygen activation in the unusual context of the tautomerase superfamily.

Crystal structures of two Bacillus carboxylesterases with different enantioselectivities.

Naproxen esterase (NP) from Bacillus subtilis Thai I-8 is a carboxylesterase that catalyzes the enantioselective hydrolysis of naproxenmethylester to produce S-naproxen (E>200). It is a homolog of CesA (98% sequence identity) and CesB (64% identity), both produced by B. subtilis strain 168. CesB can be used for the enantioselective hydrolysis of 1,2-O-isopropylideneglycerol (solketal) esters (E>200 for IPG-caprylate). Crystal structures of NP and CesB, determined to a resolution of 1.75Å and 2.04Å, respectively, showed that both proteins have a canonical α/ß hydrolase fold with an extra N-terminal helix stabilizing the cap subdomain. The active site in both enzymes is located in a deep hydrophobic groove and includes the catalytic triad residues Ser130, His274, and Glu245. A product analog, presumably 2-(2-hydroxyethoxy)acetic acid, was bound in the NP active site. The enzymes have different enantioselectivities, which previously were shown to result from only a few amino acid substitutions in the cap domain. Modeling of a substrate in the active site of NP allowed explaining the different enantioselectivities. In addition, Ala156 may be a determinant of enantioselectivity as well, since its side chain appears to interfere with the binding of certain R-enantiomers in the active site of NP. However, the exchange route for substrate and product between the active site and the solvent is not obvious from the structures. Flexibility of the cap domain might facilitate such exchange. Interestingly, both carboxylesterases show higher structural similarity to meta-cleavage compound (MCP) hydrolases than to other α/ß hydrolase fold esterases.

Crystal structure of α-1,4-glucan lyase, a unique glycoside hydrolase family member with a novel catalytic mechanism.

α-1,4-Glucan lyase (EC 4.2.2.13) from the red seaweed Gracilariopsis lemaneiformis cleaves α-1,4-glucosidic linkages in glycogen, starch, and malto-oligosaccharides, yielding the keto-monosaccharide 1,5-anhydro-D-fructose. The enzyme belongs to glycoside hydrolase family 31 (GH31) but degrades starch via an elimination reaction instead of hydrolysis. The crystal structure shows that the enzyme, like GH31 hydrolases, contains a (ß/α)8-barrel catalytic domain with B and B' subdomains, an N-terminal domain N, and the C-terminal domains C and D. The N-terminal domain N of the lyase was found to bind a trisaccharide. Complexes of the enzyme with acarbose and 1-dexoynojirimycin and two different covalent glycosyl-enzyme intermediates obtained with fluorinated sugar analogues show that, like GH31 hydrolases, the aspartic acid residues Asp(553) and Asp(665) are the catalytic nucleophile and acid, respectively. However, as a unique feature, the catalytic nucleophile is in a position to act also as a base that abstracts a proton from the C2 carbon atom of the covalently bound subsite -1 glucosyl residue, thus explaining the unique lyase activity of the enzyme. One Glu to Val mutation in the active site of the homologous α-glucosidase from Sulfolobus solfataricus resulted in a shift from hydrolytic to lyase activity, demonstrating that a subtle amino acid difference can promote lyase activity in a GH31 hydrolase.

Biochemical, kinetic, and structural characterization of a Bacillus tequilensis nitroreductase.

Nitroreductases (NRs) are NAD(P)H-dependent flavoenzymes that reduce nitro aromatic compounds to their corresponding arylamines via the nitroso and hydroxylamine intermediates. Because of their broad substrate scope and versatility, NRs have found application in multiple fields such as biocatalysis, bioremediation, cell-imaging and prodrug activation. However, only a limited number of members of the broad NR superfamily (> 24 000 sequences) have been experimentally characterized. Within this group of enzymes, only few are capable of amine synthesis, which is a fundamental chemical transformation for the pharmaceutical, agricultural, and textile industries. Herein, we provide a comprehensive description of a recently discovered NR from Bacillus tequilensis, named BtNR. This enzyme has previously been demonstrated to have the capability to fully convert nitro aromatic and heterocyclic compounds to their respective primary amines. In this study, we determined its biochemical, kinetic and structural properties, including its apparent melting temperature (Tm) of 59 °C, broad pH activity range (from pH 3 to 10) and a notably low redox potential (-236 ± 1 mV) in comparison to other well-known NRs. We also determined its steady-state and pre-steady-state kinetic parameters, which are consistent with other NRs. Additionally, we elucidated the crystal structure of BtNR, which resembles the well-characterized Escherichia coli oxygen-insensitive NAD(P)H nitroreductase (NfsB), and investigated the substrate binding in its active site through docking and molecular dynamics studies with four nitro aromatic substrates. Guided by these structural analyses, we probed the functional roles of active site residues by site-directed mutagenesis. Our findings provide valuable insights into the biochemical and structural properties of BtNR, as well as its potential applications in biotechnology.

Evolution of the catalytic mechanism at the dawn of the Baeyer-Villiger monooxygenases.

Enzymes are crucial for the emergence and sustenance of life on earth. How they became catalytically active during their evolution is still an open question. Two opposite explanations are plausible: acquiring a mechanism in a series of discrete steps or all at once in a single evolutionary event. Here, we use molecular phylogeny, ancestral sequence reconstruction, and biochemical characterization to follow the evolution of a specialized group of flavoprotein monooxygenases, the bacterial Baeyer-Villiger monooxygenases (BVMOs). These enzymes catalyze an intricate chemical reaction relying on three different elements: a reduced nicotinamide cofactor, dioxygen, and a substrate. Characterization of ancestral BVMOs shows that the catalytic mechanism evolved in a series of steps starting from a FAD-binding protein and further acquiring reactivity and specificity toward each of the elements participating in the reaction. Together, the results of our work portray how an intrinsically complex catalytic mechanism emerged during evolution.

Characterization of two bacterial multi-flavinylated proteins harboring multiple covalent flavin cofactors.

In recent years, studies have shown that a large number of bacteria secrete multi-flavinylated proteins. The exact roles and properties, of these extracellular flavoproteins that contain multiple covalently anchored FMN cofactors, are still largely unknown. Herein, we describe the biochemical and structural characterization of two multi-FMN-containing covalent flavoproteins, SaFMN3 from Streptomyces azureus and CbFMN4 from Clostridiaceae bacterium. Based on their primary structure, these proteins were predicted to contain three and four covalently tethered FMN cofactors, respectively. The genes encoding SaFMN3 and CbFMN4 were heterologously coexpressed with a flavin transferase (ApbE) in Escherichia coli, and could be purified by affinity chromatography in good yields. Both proteins were found to be soluble and to contain covalently bound FMN molecules. The SaFMN3 protein was studied in more detail and found to display a single redox potential (-184 mV) while harboring three covalently attached flavins. This is in line with the high sequence similarity when the domains of each flavoprotein are compared. The fully reduced form of SaFMN3 is able to use dioxygen as electron acceptor. Single domains from both proteins were expressed, purified and crystallized. The crystal structures were elucidated, which confirmed that the flavin cofactor is covalently attached to a threonine. Comparison of both crystal structures revealed a high similarity, even in the flavin binding pocket. Based on the crystal structure, mutants of the SaFMN3-D2 domain were designed to improve its fluorescence quantum yield by changing the microenvironment of the isoalloxazine moiety of the flavin cofactor. Residues that quench the flavin fluorescence were successfully identified. Our study reveals biochemical details of multi-FMN-containing proteins, contributing to a better understanding of their role in bacteria and providing leads to future utilization of these flavoprotein in biotechnology.

Identification and characterization of archaeal and bacterial F420 -dependent thioredoxin reductases.

The thioredoxin pathway is an antioxidant system present in most organisms. Electrons flow from a thioredoxin reductase to thioredoxin at the expense of a specific electron donor. Most known thioredoxin reductases rely on NADPH as a reducing cofactor. Yet, in 2016, a new type of thioredoxin reductase was discovered in Archaea which utilize instead a reduced deazaflavin cofactor (F420 H2 ). For this reason, the respective enzyme was named deazaflavin-dependent flavin-containing thioredoxin reductase (DFTR). To have a broader understanding of the biochemistry of DFTRs, we identified and characterized two other archaeal representatives. A detailed kinetic study, which included pre-steady state kinetic analyses, revealed that these two DFTRs are highly specific for F420 H2 while displaying marginal activity with NADPH. Nevertheless, they share mechanistic features with the canonical thioredoxin reductases that are dependent on NADPH (NTRs). A detailed structural analysis led to the identification of two key residues that tune cofactor specificity of DFTRs. This allowed us to propose a DFTR-specific sequence motif that enabled for the first time the identification and experimental characterization of a bacterial DFTR.

Regio- and stereoselective steroid hydroxylation by CYP109A2 from Bacillus megaterium explored by X-ray crystallography and computational modeling.

The P450 monooxygenase CYP109A2 from Bacillus megaterium DSM319 was previously found to convert vitamin D3 (VD3) to 25-hydroxyvitamin D3. Here, we show that this enzyme is also able to convert testosterone in a highly regio- and stereoselective manner to 16ß-hydroxytestosterone. To reveal the structural determinants governing the regio- and stereoselective steroid hydroxylation reactions catalyzed by CYP109A2, two crystal structures of CYP109A2 were solved in similar closed conformations, one revealing a bound testosterone in the active site pocket, albeit at a nonproductive site away from the heme-iron. To examine whether the closed crystal structures nevertheless correspond to a reactive conformation of CYP109A2, docking and molecular dynamics (MD) simulations were performed with testosterone and vitamin D3 (VD3) present in the active site. These MD simulations were analyzed for catalytically productive conformations, the relative occurrences of which were in agreement with the experimentally determined stereoselectivities if the predicted stability of each carbon-hydrogen bond was taken into account. Overall, the first-time determination and analysis of the catalytically relevant 3D conformation of CYP109A2 will allow for future small molecule ligand screening in silico, as well as enabling site-directed mutagenesis toward improved enzymatic properties of this enzyme.