Enzymatic mechanisms for catalysis of enolization: ketosteroid isomerase

Noncanonical Amino Acids in Biocatalysis

In recent years, powerful genetic code reprogramming methods have emerged that allow new functional components to be embedded into proteins as noncanonical amino acid (ncAA) side chains. In this review, we will illustrate how the availability of an expanded set of amino acid building blocks has opened a wealth of new opportunities in enzymology and biocatalysis research. Genetic code reprogramming has provided new insights into enzyme mechanisms by allowing introduction of new spectroscopic probes and the targeted replacement of individual atoms or functional groups. NcAAs have also been used to develop engineered biocatalysts with improved activity, selectivity, and stability, as well as enzymes with artificial regulatory elements that are responsive to external stimuli. Perhaps most ambitiously, the combination of genetic code reprogramming and laboratory evolution has given rise to new classes of enzymes that use ncAAs as key catalytic elements. With the framework for developing ncAA-containing biocatalysts now firmly established, we are optimistic that genetic code reprogramming will become a progressively more powerful tool in the armory of enzyme designers and engineers in the coming years.

Computational Design of an Electro-Organocatalyst for Conversion of CO2 into Formaldehyde

Density functional theory calculations employing a hybrid implicit/explicit solvation method were used to explore a new strategy for electrochemical conversion of CO2 using an electro-organocatalyst. A particular structural motif is identified that consists of an electron-rich vicinal enediamine (>N-C═C-N<) backbone, which is capable of activating CO2 by the formation of a C-C bond while subsequently facilitating the transfer of electrons from a chemically inert cathode to ultimately produce formaldehyde. Unlike transition metal-based electrocatalysts, the electro-organocatalyst is not constrained by scaling relations between the formation energies of activated CO2 and adsorbed CO, nor is it expected to be active for the competing hydrogen evolution reaction. The rate-limiting steps are found to occur during two proton-coupled electron transfer (PCET) sequences and are associated with the transfer of a proton from a proton transfer mediator to a carbon atom on the electro-organocatalyst. The difficulty of this step in the second PCET sequence necessitates an electrode potential of -0.85 V vs RHE to achieve the maximum turnover frequency. In addition, it is postulated that the electro-organocatalyst should also be capable of forming long-chain aldehydes by successively carrying out reductive aldol condensation to grow the alkyl chain one carbon at a time.

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Potent GST Ketosteroid Isomerase Activity Relevant to Ecdysteroidogenesis in the Malaria Vector Anopheles gambiae

Nobo is a glutathione transferase (GST) crucially contributing to ecdysteroid biosynthesis in insects of the orders Diptera and Lepidoptera. Ecdysone is a vital steroid hormone in insects, which governs larval molting and metamorphosis, and the suppression of its synthesis has potential as a novel approach to insect growth regulation and combatting vectors of disease. In general, GSTs catalyze detoxication, whereas the specific function of Nobo in ecdysteroidogenesis is unknown. We report that Nobo from the malaria-spreading mosquito Anopheles gambiae is a highly efficient ketosteroid isomerase catalyzing double-bond isomerization in the steroids 5-androsten-3,17-dione and 5-pregnen-3,20-dione. These mammalian ketosteroids are unknown in mosquitoes, but the discovered prominent catalytic activity of these compounds suggests that the unknown Nobo substrate in insects has a ketosteroid functionality. Aminoacid residue Asp111 in Nobo is essential for activity with the steroids, but not for conventional GST substrates. Further characterization of Nobo may guide the development of new insecticides to prevent malaria.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Carboxylate Catalyzed Isomerization of β,γ‐Unsaturated N ‐Acetylcysteamine Thioesters**

We demonstrate herein the capacity of simple carboxylate salts – tetrametylammonium and tetramethylguanidinium pivalate – to act as catalysts in the isomerization of β,γ‐unsaturated thioesters to α,β‐unsaturated thioesters. The carboxylate catalysts gave reaction rates comparable to those obtained with DBU, but with fewer side reactions. The reaction exhibits a normal secondary kinetic isotope effect (k_1H/k_1D=1.065±0.026) with a β,γ‐deuterated substrate. Computational analysis of the mechanism provides a similar value (k_1H/k_1D=1.05) with a mechanism where γ‐reprotonation of the enolate intermediate is rate determining.

Oriented internal electrostatic fields: an emerging design element in coordination chemistry and catalysis

The power of oriented electrostatic fields (ESFs) to influence chemical bonding and reactivity is a phenomenon of rapidly growing interest. The presence of strong ESFs has recently been implicated as one of the most significant contributors to the activity of select enzymes, wherein alignment of a substrate's changing dipole moment with a strong, local electrostatic field has been shown to be responsible for the majority of the enzymatic rate enhancement. Outside of enzymology, researchers have studied the impacts of "internal" electrostatic fields via the addition of ionic salts to reactions and the incorporation of charged functional groups into organic molecules (both experimentally and computationally), and "externally" via the implementation of bulk fields between electrode plates. Incorporation of charged moieties into homogeneous inorganic complexes to generate internal ESFs represents an area of high potential for novel catalyst design. This field has only begun to materialize within the past 10 years but could be an area of significant impact moving forward, since it provides a means for tuning the properties of molecular complexes via a method that is orthogonal to traditional strategies, thereby providing possibilities for improved catalytic conditions and novel reactivity. In this perspective, we highlight recent developments in this area and offer insights, obtained from our own research, on the challenges and future directions of this emerging field of research.

Glycerol-3-Phosphate Dehydrogenase: The K120 and K204 Side Chains Define an Oxyanion Hole at the Enzyme Active Site

The cationic K120 and K204 side chains lie close to the C-2 carbonyl group of substrate dihydroxyacetone phosphate (DHAP) at the active site of glycerol-3-phosphate dehydrogenase (GPDH), and the K120 side chain is also positioned to form a hydrogen bond to the C-1 hydroxyl of DHAP. The kinetic parameters for unactivated and phosphite dianion-activated GPDH-catalyzed reduction of glycolaldehyde and acetaldehyde (AcA) show that the transition state for the former reaction is stabilized by ca 5 kcal/mole by interactions of the C-1 hydroxyl group with the protein catalyst. The K120A and K204A substitutions at wild-type GPDH result in similar decreases in k_cat, but K_m is only affected by the K120A substitution. These results are consistent with 3 kcal/mol stabilizing interactions between the K120 or K204 side chains and a negative charge at the C-2 oxygen at the transition state for hydride transfer from NADH to DHAP. This stabilization resembles that observed at oxyanion holes for other enzymes. There is no detectable rescue of the K204A variant by ethylammonium cation (EtNH₃⁺), compared with the efficient rescue of the K120A variant. This is consistent with a difference in the accessibility of the variant enzyme active sites to exogenous EtNH₃⁺. The K120A/K204A substitutions cause a (6 × 10⁶)-fold increase in the promiscuity of wild-type hlGPDH for catalysis of the reduction of AcA compared to DHAP. This may reflect conservation of the active site for an ancestral alcohol dehydrogenase, whose relative activity for catalysis of reduction of AcA increases with substitutions that reduce the activity for reduction of the specific substrate DHAP.

Structure and chemistry of enzymatic active sites that play a role in the switch and conformation mechanism

Enzymes, which are biological molecules, are constructed from polypeptide chains, and these molecules are activated through reaction mechanisms. It is the role of enzymes to speed up chemical reactions that are used to build or break down cell structures. Activation energy is reduced by the enzymes' selective binding of substrates in a protected environment. In enzyme tertiary structures, the active sites are commonly situated in a "cleft," which necessitates the diffusion of substrates and products. The amino acid residues of the active site may be far apart in the primary structure owing to the folding required for tertiary structure. Due to their critical role in substrate binding and attraction, changes in amino acid structure at or near the enzyme's active site usually alter enzyme activity. At the enzyme's active site, or where the chemical reactions occur, the substrate is bound. Enzyme substrates are the primary targets of the enzyme's active site, which is designed to assist in the chemical reaction. This chapter elucidates the summary of structure and chemistry of enzymes, their active site features, charges and role of water in the structures to clarify the biochemistry of the enzymes in the depth of atomic features.

Glutathione Transferases as Efficient Ketosteroid Isomerases

In addition to their well-established role in detoxication, glutathione transferases (GSTs) have other biological functions. We are focusing on the ketosteroid isomerase activity, which appears to contribute to steroid hormone biosynthesis in mammalian tissues. A highly efficient GST A3-3 is present in some, but not all, mammals. The alpha class enzyme GST A3-3 in humans and the horse shows the highest catalytic efficiency with kcat/Km values of approximately 107 M-1s-1, ranking close to the most active enzymes known. The expression of GST A3-3 in steroidogenic tissues suggests that the enzyme has evolved to support the activity of 3β-hydroxysteroid dehydrogenase, which catalyzes the formation of 5-androsten-3,17-dione and 5-pregnen-3,20-dione that are substrates for the double-bond isomerization catalyzed by GST A3-3. The dehydrogenase also catalyzes the isomerization, but its kcat of approximately 1 s-1 is 200-fold lower than the kcat values of human and equine GST A3-3. Inhibition of GST A3-3 in progesterone-producing human cells suppress the formation of the hormone. Glutathione serves as a coenzyme contributing a thiolate as a base in the isomerase mechanism, which also involves the active-site Tyr9 and Arg15. These conserved residues are necessary but not sufficient for the ketosteroid isomerase activity. A proper assortment of H-site residues is crucial to efficient catalysis by forming the cavity binding the hydrophobic substrate. It remains to elucidate why some mammals, such as rats and mice, lack GSTs with the prominent ketosteroid isomerase activity found in certain other species. Remarkably, the fruit fly Drosophila melanogaster, expresses a GSTE14 with notable steroid isomerase activity, even though Ser14 has evolved as the active-site residue corresponding to Tyr9 in the mammalian alpha class.

Globally Optimized Molecular Embeddings for Dynamic Reaction Solvate Shell Optimization and Active Site Design

We demonstrate how a full QM/MM derivatization of the recently developed GOCAT model can be utilized in the global optimization of molecular embeddings. To this end, we provide two distinct examples: An $$\text {S}_\text {N}2$$ S N 2 reaction, and one enzymatic example of recent interest, the ketosteroid isomerase. These serve us to highlight the advantages of such an approach and sketch the roadmap for further improvements.

Structural characterization and Kemp eliminase activity of the Mycobacterium smegmatis Ketosteroid Isomerase

The Kemp elimination reaction, involving the ring-opening of benzoxazole and its derivatives under the action of natural enzymes or chemical catalysts, has been of interest to researchers since its discovery. Because this reaction does not exist in all currently known metabolic pathways, the computational design of Kemp eliminases has provided valuable insights into principles of enzymatic catalysis. However, it was discovered that the naturally occurring promiscuous enzymes ydbC, xapA and ketosteroid isomerase also can catalyze Kemp elimination. Here, we report the crystal structure of ketosteroid isomerase (KSI) from Mycobacterium smegmatis MC2 155. MsKSI crystallizes in the P2₁2₁2₁ space group with two molecules in an asymmetric unit, and ultracentrifugation data confirms that it forms a stable dimer in solution, consistent with the 1.9 Å-resolution structure. Our assays confirm that MsKSI accelerates the Kemp elimination of 5-nitrobenzoxazole (5NBI) with an optimal pH of 5.5. A 2.35 Å resolution crystal structure of the MsKSI-5NBI complex reveals that the substrate 5NBI is bound in the active pocket of the enzyme composed of hydrophobic residues. In addition, the Glu127 residue is proposed to play an important role as a general base in proton transfer and breaking weak O-N bonds to open the five-membered ring. This work provides a starting point for exploring the artificial modification of MsKSI using the natural enzyme as the backbone.

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts

As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.

Biosynthesis of the Immunosuppressant (-)-FR901483

We report characterization of the biosynthetic pathway of the potent immunosuppressant (-)-FR901483 (1) through heterologous expression and enzymatic assays. The biosynthetic logic to form the azatricyclic alkaloid is consistent with those proposed in biomimetic syntheses and involves aza-spiro annulation of dityrosyl-piperazine to form a ketoaldehyde intermediate, followed by regioselective aldol condensation, stereoselective ketoreduction, and phosphorylation. A possible target of 1 is proposed based on the biosynthetic studies.

BIMP-Catalyzed 1,3-Prototropic Shift for the Highly Enantioselective Synthesis of Conjugated Cyclohexenones

A bifunctional iminophosphorane (BIMP)-catalysed enantioselective synthesis of α,β-unsaturated cyclohexenones through a facially selective 1,3-prototropic shift of β,γ-unsaturated prochiral isomers, under mild reaction conditions and in short reaction times, on a range of structurally diverse substrates, is reported. α,β-Unsaturated cyclohexenone products primed for downstream derivatisation were obtained in high yields (up to 99 %) and consistently high enantioselectivity (up to 99 % ee). Computational studies into the reaction mechanism and origins of enantioselectivity, including multivariate linear regression of TS energy, were carried out and the obtained data were found to be in good agreement with experimental findings.

PR10/Bet v1-like Proteins as Novel Contributors to Plant Biochemical Diversity

Pathogenesis-related (PR) proteins constitute a broad class of plant proteins with analogues found throughout nature from bacteria to higher eukaryotes. PR proteins were first noted in plants as part of the hypersensitive response, but have since been assigned an array of biological roles. The PR10/Bet v1-like proteins are a subset of PR proteins characterized by an ability to bind a wide range of lipophilic ligands, uniquely positioning them as contributors to specialized biosynthetic pathways. PR10/Bet v1-like proteins participate in the production of plant alkaloids and phenolics including flavonoids, both as general binding proteins and in special cases as catalysts. Owing initially to the perceived allergenic properties of PR10/Bet v1-like proteins, many were studied at the structural level to elucidate the basis for ligand binding. These studies provided a foundation for more recent efforts to understand higher-level structural order and how PR10/Bet v1-like proteins catalyse key reactions in plant pathways. Synthetic biology aimed at reconstituting plant-specialized metabolism in microorganisms uses knowledge of these proteins to fine-tune performance in new systems.

A Preorganized Electric Field Leads to Minimal Geometrical Reorientation in the Catalytic Reaction of Ketosteroid Isomerase

Electrostatic interactions play a pivotal role in enzymatic catalysis and are increasingly modeled explicitly in computational enzyme design; nevertheless, they are challenging to measure experimentally. Using vibrational Stark effect (VSE) spectroscopy, we have measured electric fields inside the active site of the enzyme ketosteroid isomerase (KSI). These studies have shown that these fields can be unusually large, but it has been unclear to what extent they specifically stabilize the transition state (TS) relative to a ground state (GS). In the following, we use crystallography and computational modeling to show that KSI's intrinsic electric field is nearly perfectly oriented to stabilize the geometry of its reaction's TS. Moreover, we find that this electric field adjusts the orientation of its substrate in the ground state so that the substrate needs to only undergo minimal structural changes upon activation to its TS. This work provides evidence that the active site electric field in KSI is preorganized to facilitate catalysis and provides a template for how electrostatic preorganization can be measured in enzymatic systems.

On the Role of Water as a Catalyst in Prebiotic Chemistry

In this manuscript we provide computational support to the catalytic role of water in all kinds of pseudopericyclic reactions operating in the reductive acid cycle, as well as in other metabolic processes. Water catalysis is not limited to those reactions where simple translocation of hydrogen atoms occurs, such as those represented by NuH+E→Nu-EH general equation. Indeed, water catalysis is more general and extremely important in tautomerization reactions of the type HX-Y=Z→X=Y-ZH, which operate in the reductive citric acid cycle and metabolic processes. Moreover, the comprehensive theoretical study reported herein illustrates that these reactions appear to behave as authentic enzyme-catalyzed reactions showing Michaelis-Menten behavior, however with the abnormal singularity that the concentration of the catalytic "water clusters" of different length and nature must be taken as a huge number. Overall, the results presented are suggestive of the workability of the so-called "metabolism first" proposal in a hot water world, as water catalysis eliminates the dilution problem frequently associated to this proposal.

The role of Brønsted base basicity in estimating carbon acidity at enzyme active sites: a caveat

Many enzymes catalyze the abstraction of a proton from a carbon acid substrate to initiate a variety of reactions; however, the development of a complete quantitative description of enzyme-catalyzed heterolytic cleavage of a C-H bond remains a challenge to enzymologists. To determine the pK value for such substrates bound at the active site, recent studies have estimated the equilibrium for formation of the deprotonated intermediate at the active site, however, accurate knowledge of the pK of the conjugate acid of the Brønsted base catalyst (BH⁺) is also required. Herein, it is shown that using the value of pK of the enzyme-substrate complex can underestimate the value of pK by an amount between zero and pδ, where pδ is the change in basicity of BH⁺ upon going from the enzyme-substrate complex to the enzyme-intermediate complex.

Structural basis for human sterol isomerase in cholesterol biosynthesis and multidrug recognition

3-β-hydroxysteroid-Δ8, Δ7-isomerase, known as Emopamil-Binding Protein (EBP), is an endoplasmic reticulum membrane protein involved in cholesterol biosynthesis, autophagy, oligodendrocyte formation. The mutation on EBP can cause Conradi-Hunermann syndrome, an inborn error. Interestingly, EBP binds an abundance of structurally diverse pharmacologically active compounds, causing drug resistance. Here, we report two crystal structures of human EBP, one in complex with the anti-breast cancer drug tamoxifen and the other in complex with the cholesterol biosynthesis inhibitor U18666A. EBP adopts an unreported fold involving five transmembrane-helices (TMs) that creates a membrane cavity presenting a pharmacological binding site that accommodates multiple different ligands. The compounds exploit their positively-charged amine group to mimic the carbocationic sterol intermediate. Mutagenesis studies on specific residues abolish the isomerase activity and decrease the multidrug binding capacity. This work reveals the catalytic mechanism of EBP-mediated isomerization in cholesterol biosynthesis and how this protein may act as a multi-drug binder.

Elucidation of the biosynthetic pathways of boldenone in the equine testis

Boldenone is an anabolic-androgenic steroid that is prohibited in equine sports. Urine from the uncastrated male horse contains boldenone that is thought to be of endogenous origin and thus a threshold ('cut-off') concentration has been adopted internationally for free and conjugated boldenone to help distinguish cases of doping from its natural production. The testis is likely to be a source of boldenone. Qualitative analysis was performed on extracts of equine testicular homogenates (n = 3 horses) incubated non-spiked and in the presence of its potential precursors using liquid chromatography tandem mass spectrometry (LC-MS/MS) and LC high resolution mass spectrometry (LC-HRMS). Samples were analysed both underivatised and derivatised to increase the certainty of identification. In addition to previously reported endogenous steroids, analysis of non-spiked testicular tissue samples demonstrated the presence of boldenone and boldienone at trace levels in the equine testis. Incubation of homogenates with deuterium or carbon isotope labelled testosterone and androstenedione resulted in the matching stable isotope analogues of boldenone and boldienone being formed. Additionally, deuterium and carbon labelled 2-hydroxyandrostenedione was detected, raising the possibility that this steroid is a biosynthetic intermediate. In conclusion, boldenone and boldienone are naturally present in the equine testis, with the biosynthesis of these steroids arising from the conversion of testosterone and androstenedione. However, additional work employing larger numbers of animals, further enzyme kinetic experiments and pure reference standards for 2-OH androstenedione isomers would be required to better characterize the pathways involved in these transformations.

Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis

The enormous rate accelerations observed for many enzyme catalysts are due to strong stabilizing interactions between the protein and reaction transition state. The defining property of these catalysts is their specificity for binding the transition state with a much higher affinity than substrate. Experimental results are presented which show that the phosphodianion-binding energy of phosphate monoester substrates is used to drive conversion of their protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes. These results are generalized to other enzyme-catalyzed reactions. The existence of many enzymes in flexible, entropically rich, and inactive ground states provides a mechanism for utilization of ligand-binding energy to mold these catalysts into stiff and active forms. This reduces the substrate-binding energy expressed at the Michaelis complex, while enabling the full and specific expression of large transition-state binding energies. Evidence is presented that the complexity of enzyme conformational changes increases with increases in the enzymatic rate acceleration. The requirement that a large fraction of the total substrate-binding energy be utilized to drive conformational changes of floppy enzymes is proposed to favor the selection and evolution of protein folds with multiple flexible unstructured loops, such as the TIM-barrel fold. The effect of protein motions on the kinetic parameters for enzymes that undergo ligand-driven conformational changes is considered. The results of computational studies to model the complex ligand-driven conformational change in catalysis by triosephosphate isomerase are presented.

An Ab Initio QM/MM Study of the Electrostatic Contribution to Catalysis in the Active Site of Ketosteroid Isomerase

The electric field in the hydrogen-bond network of the active site of ketosteroid isomerase (KSI) has been experimentally measured using vibrational Stark effect (VSE) spectroscopy, and utilized to study the electrostatic contribution to catalysis. A large gap was found in the electric field between the computational simulation based on the Amber force field and the experimental measurement. In this work, quantum mechanical (QM) calculations of the electric field were performed using an ab initio QM/MM molecular dynamics (MD) simulation and electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method. Our results demonstrate that the QM-derived electric field based on the snapshots from QM/MM MD simulation could give quantitative agreement with the experiment. The accurate calculation of the electric field inside the protein requires both the rigorous sampling of configurations, and a QM description of the electrostatic field. Based on the direct QM calculation of the electric field, we theoretically confirmed that there is a linear correlation relationship between the activation free energy and the electric field in the active site of wild-type KSI and its mutants (namely, D103N, Y16S, and D103L). Our study presents a computational protocol for the accurate simulation of the electric field in the active site of the protein, and provides a theoretical foundation that supports the link between electric fields and enzyme catalysis.

Monitoring dehydroepiandrosterone (DHEA) in the urine of Thoroughbred geldings for doping control purposes

The use of testosterone and its pro-drugs, such as dehydroepiandrosterone (DHEA), is currently regulated in horseracing by the application of international testosterone thresholds. However, additional steroidomic approaches, such as steroid ratios, to distinguish overall adrenal stimulation from drug administrations and an equine biological passport for longitudinal steroid profiling of individual animals could be advantageous in equine doping testing. Thus, DHEA concentrations and related ratios (testosterone [T] to DHEA and DHEA to epitestosterone [E]) were assessed in the reference population by quantitative analysis of 200 post-race gelding urine samples using liquid chromatography-tandem mass spectrometry. DHEA concentrations ranged between 0.9 and 136.6 ng/mL (mean 12.8 ng/mL), T:DHEA ratios between 0.06 and 1.85 (mean 0.43), and DHEA:E ratios between 0.21 and 13.56 (mean 2.20). Based on the reference population statistical upper limits of 5.4 for T:DHEA ratio and 48.1 for DHEA:E ratio are proposed with a risk of 1 in 10 000 for a normal outlier exceeding the value. Analysis of post-administration urine samples collected following administrations of DHEA, Equi-Bolic® (a mix of DHEA and pregnenolone) and testosterone propionate to geldings showed that the upper limit for T:DHEA ratio was exceeded following testosterone propionate administration and DHEA:E ratio following DHEA administrations and thus these ratios could be used as additional biomarkers when determining the cause of an atypical testosterone concentration. Additionally, DHEA concentrations and ratios can be used as a starting point to establish reference ranges for an equine biological passport.

Engineering the "Missing Link" in Biosynthetic (-)-Menthol Production: Bacterial Isopulegone Isomerase

The realization of a synthetic biology approach to microbial (1R,2S,5R)-(−)-menthol (1) production relies on the identification of a gene encoding an isopulegone isomerase (IPGI), the only enzyme in the Mentha piperita biosynthetic pathway as yet unidentified. We demonstrate that Δ5-3-ketosteroid isomerase (KSI) from Pseudomonas putida can act as an IPGI, producing (R)-(+)-pulegone ((R)-2) from (+)-cis-isopulegone (3). Using a robotics-driven semirational design strategy, we identified a key KSI variant encoding four active site mutations, which confer a 4.3-fold increase in activity over the wild-type enzyme. This was assisted by the generation of crystal structures of four KSI variants, combined with molecular modeling of 3 binding to identify key active site residue targets. The KSI variant was demonstrated to function efficiently within cascade biocatalytic reactions with downstream Mentha enzymes pulegone reductase and (−)-menthone:(−)-menthol reductase to generate 1 from 3. This study introduces the use of a recombinant IPGI, engineered to function efficiently within a biosynthetic pathway for the production of 1 in microorganisms.

Electric Fields and Fast Protein Dynamics in Enzymes

In recent years, there has been much discussion regarding the origin of enzymatic catalysis and whether including protein dynamics is necessary for understanding catalytic enhancement. An important contribution in this debate was made with the application of the vibrational Stark effect spectroscopy to measure electric fields in the active site. This provided a window on electric fields at the transition state in enzymatic reactions. We performed computational studies on two enzymes where we have shown that fast dynamics is part of the reaction mechanism and calculated the electric field near the bond-breaking event. We found that the fast motions that we had identified lead to an increase of the electric field, thus preparing an enzymatic configuration that is electrostatically favorable for the catalytic chemical step. We also studied the enzyme that has been the subject of Stark spectroscopy, ketosteroid isomerase, and found electric fields of a similar magnitude to the two previous examples.

Enormous Hydrogen Bond Strength Enhancement through π-Conjugation Gain: Implications for Enzyme Catalysis

Surprisingly large resonance-assistance effects may explain how some enzymes form extremely short, strong hydrogen bonds to stabilize reactive oxyanion intermediates and facilitate catalysis. Computational models for several enzymic residue-substrate interactions reveal that when a π-conjugated, hydrogen bond donor (XH) forms a hydrogen bond to a charged substrate (Y^-), XH can become significantly more π-electron delocalized, and this "extra" stabilization may boost the [XH···Y^-] hydrogen bond strength by ≥15 kcal/mol. This reciprocal relationship departs from the widespread pK_a concept (i.e., the idea that short, strong hydrogen bonds form when the interacting moieties have matching pK_a values), which has been the rationale for enzymic acid-base reactions. The findings presented here provide new insight into how short, strong hydrogen bonds could form in enzymes.

Electric Fields and Enzyme Catalysis

What happens inside an enzyme's active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists' attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme's active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

Methyl and p-Bromobenzyl Esters of Hydrogenated Kaurenoic Acid for Controlling Anthracnose in Common Bean Plants

Kaurenoic acid derivatives were prepared and submitted to in vitro assays with the fungus Colletotrichum lindemuthianum, which causes anthracnose disease in the common bean. The most active substances were found to be methyl and p-bromobenzylesters, 7 and 9, respectively, of the hydrogenated kaurenoic acid, which presented a minimum inhibitory concentration (MIC) of 0.097 and 0.131 mM, respectively, while the commercial fungicide methyl thiophanate (MT) presented a MIC of 0.143 mM. Substances 7 (1.401 mM) and 9 (1.886 mM) reduced the severity of anthracnose in common bean to values statistically comparable to MT (2.044 mM). According to an in silico study, both compounds 7 and 9 are inhibitors of the ketosteroid isomerase (KSI) enzyme produced by other organisms, the amino acid sequence of which could be detected in fungal genomes. These substances appeared to act against C. lindemuthianum by inhibiting its KSI. Therefore, substances 7 and 9 are promising for the development of new fungicides.

Kemp Eliminase Activity of Ketosteroid Isomerase

Kemp eliminases represent the most successful class of computationally designed enzymes, with rate accelerations of up to 10⁹-fold relative to the rate of the same reaction in aqueous solution. Nevertheless, several other systems such as micelles, catalytic antibodies, and cavitands are known to accelerate the Kemp elimination by several orders of magnitude. We found that the naturally occurring enzyme ketosteroid isomerase (KSI) also catalyzes the Kemp elimination. Surprisingly, mutations of D38, the residue that acts as a general base for its natural substrate, produced variants that catalyze the Kemp elimination up to 7000-fold better than wild-type KSI does, and some of these variants accelerate the Kemp elimination more than the computationally designed Kemp eliminases. Analysis of the D38N general base KSI variant suggests that a different active site carboxylate residue, D99, performs the proton abstraction. Docking simulations and analysis of inhibition by active site binders suggest that the Kemp elimination takes place in the active site of KSI and that KSI uses the same catalytic strategies of the computationally designed enzymes. In agreement with prior observations, our results strengthen the conclusion that significant rate accelerations of the Kemp elimination can be achieved with very few, nonspecific interactions with the substrate if a suitable catalytic base is present in a hydrophobic environment. Computational design can fulfill these requirements, and the design of more complex and precise environments represents the next level of challenges for protein design.

Molecular Seesaw: How Increased Hydrogen Bonding Can Hinder Excited-State Proton Transfer

A previously unexplained effect in the relative rate of excited-state intramolecular proton transfer (ESIPT) in related indole derivatives is investigated using both theory and experiment. Ultrafast spectroscopy [ J. Phys. Chem. A, 2015, 119, 5618-5625 ] found that although the diol 1,3-bis(2-pyridylimino)-4,7-dihydroxyisoindole exhibits two equivalent intramolecular hydrogen bonds, the ESIPT rate associated with tautomerization of either hydrogen bond is a factor of 2 slower than that of the single intramolecular hydrogen bond in the ethoxy-ol 1,3-bis(2-pyridylimino)-4-ethoxy-7-hydroxyisoindole. Excited-state electronic structure calculations suggest a resolution to this puzzle by revealing a seesaw effect in which the two hydrogen bonds of the diol are both longer than the single hydrogen bond in the ethoxy-ol. Semiclassical rate theory recovers the previously unexplained trends and leads to clear predictions regarding the relative H/D kinetic isotope effect (KIE) for ESIPT in the two systems. The theoretical KIE predictions are tested using ultrafast spectroscopy, confirming the seesaw effect.

A Critical Test of the Electrostatic Contribution to Catalysis with Noncanonical Amino Acids in Ketosteroid Isomerase

The vibrational Stark effect (VSE) has been used to measure the electric field in the active site of ketosteroid isomerase (KSI). These measured fields correlate with ΔG(⧧) in a series of conventional mutants, yielding an estimate for the electrostatic contribution to catalysis (Fried et al. Science 2014, 346, 1510-1513). In this work we test this result with much more conservative variants in which individual Tyr residues in the active site are replaced by 3-chlorotyrosine via amber suppression. The electric fields sensed at the position of the carbonyl bond involved in charge displacement during catalysis were characterized using the VSE, where the field sensitivity has been calibrated by vibrational Stark spectroscopy, solvatochromism, and MD simulations. A linear relationship is observed between the electric field and ΔG(⧧) that interpolates between wild-type and more drastic conventional mutations, reinforcing the evaluation of the electrostatic contribution to catalysis in KSI. A simplified model and calculation are developed to estimate changes in the electric field accompanying changes in the extended hydrogen-bond network in the active site. The results are consistent with a model in which the O-H group of a key active site tyrosine functions by imposing a static electrostatic potential onto the carbonyl bond. The model suggests that the contribution to catalysis from the active site hydrogen bonds is of similar weight to the distal interactions from the rest of the protein. A similar linear correlation was also observed between the proton affinity of KSI's active site and the catalytic rate, suggesting a direct connection between the strength of the H-bond and the electric field it exerts.

Role of conserved Met112 residue in the catalytic activity and stability of ketosteroid isomerase

Ketosteroid isomerase (3-oxosteroid Δ(5)-Δ(4)-isomerase, KSI) from Pseudomonas putida catalyzes allylic rearrangement of the 5,6-double bond of Δ(5)-3-ketosteroid to 4,5-position by stereospecific intramolecular transfer of a proton. The active site of KSI is formed by several hydrophobic residues and three catalytic residues (Tyr14, Asp38, and Asp99). In this study, we investigated the role of a hydrophobic Met112 residue near the active site in the catalysis, steroid binding, and stability of KSI. Replacing Met112 with alanine (yields M112A) or leucine (M112L) decreased the kcat by 20- and 4-fold, respectively. Compared with the wild type (WT), M112A and M112L KSIs showed increased KD values for equilenin, an intermediate analogue; these changes suggest that loss of packing at position 112 might lead to unfavorable steroid binding, thereby resulting in decreased catalytic activity. Furthermore, M112A and M112L mutations reduced melting temperature (Tm) by 6.4°C and 2.5°C, respectively. These changes suggest that favorable packing in the core is important for the maintenance of stability in KSI. The M112K mutation decreased kcat by 2000-fold, compared with the WT. In M112K KSI structure, a new salt bridge was formed between Asp38 and Lys112. This bridge could change the electrostatic potential of Asp38, and thereby contribute to the decreased catalytic activity. The M112K mutation also decreased the stability by reducing Tm by 4.1°C. Our data suggest that the Met112 residue may contribute to the catalytic activity and stability of KSI by providing favorable hydrophobic environments and compact packing in the catalytic core.

Asymmetric Catalysis with CO2 : The Direct α-Allylation of Ketones

Quaternary stereocenters are found in numerous bioactive molecules. The Tsuji-Trost reaction has proven to be a powerful C-C bond forming process, and, at least in principle, should be well suited to access quaternary stereocenters via the α-allylation of ketones. However, while indirect approaches are known, the direct, catalytic asymmetric α-allylation of branched ketones has been elusive until today. By combining "enol catalysis" with the use of CO2 as a formal catalyst for asymmetric catalysis, we have now developed a solution to this problem: we report a direct, highly enantioselective and highly atom-economic Tsuji-Trost allylation of branched ketones with allylic alcohol. Our reaction delivers products bearing quaternary stereocenters with high enantioselectivity and water as the sole by-product. We expect our methodology to be of utility in asymmetric catalysis and inspire the design of other highly atom-economic transformations.

Theoretical study of enzymatically catalyzed tautomerization of carbon acids in aqueous solution: quantum calculations and steered molecular dynamics simulations

The thermodynamics and kinetics of enzymatically assisted reactions of carbon acids were studied theoretically in this work. Quantum electronic (QE) structure calculations and steered molecular dynamics (SMD) simulations were carried out. Three 3-butenal tautomerization reactions that proceed from the β,γ-unsaturated reactant (R) to the α,β-unsaturated carbon acid product (P) and occur in two elementary steps through an intermediate (I) were studied, ignoring or including the surrounding aqueous medium in the calculations. The Gibbs free energies of activation of the R ⇆ I enolization and I ⇆ P ketonization steps were found to decrease considerably when residues simulating enzymes were introduced into these processes. Although the processes became slightly more favorable thermodynamically when the solution was included in the simulations, they became less favorable kinetically. The results from SMD simulations of these reactions were qualitatively consistent with the values we obtained using QE as well as those found by other authors in similar studies. Our simulations also allowed us to perform a detailed study of these reactions in solution.

Role of Loop-Clamping Side Chains in Catalysis by Triosephosphate Isomerase

The side chains of Y208 and S211 from loop 7 of triosephosphate isomerase (TIM) form hydrogen bonds to backbone amides and carbonyls from loop 6 to stabilize the caged enzyme–substrate complex. The effect of seven mutations [Y208T, Y208S, Y208A, Y208F, S211G, S211A, Y208T/S211G] on the kinetic parameters for TIM catalyzed reactions of the whole substrates dihydroxyacetone phosphate and d-glyceraldehyde 3-phosphate [(k_cat/K_m)_GAP and (k_cat/K_m)_DHAP] and of the substrate pieces glycolaldehyde and phosphite dianion (k_cat/K_HPiK_GA) are reported. The linear logarithmic correlation between these kinetic parameters, with slope of 1.04 ± 0.03, shows that most mutations of TIM result in an identical change in the activation barriers for the catalyzed reactions of whole substrate and substrate pieces, so that the transition states for these reactions are stabilized by similar interactions with the protein catalyst. The second linear logarithmic correlation [slope = 0.53 ± 0.16] between k_cat for isomerization of GAP and K_d^⧧ for phosphite dianion binding to the transition state for wildtype and many mutant TIM-catalyzed reactions of substrate pieces shows that ca. 50% of the wildtype TIM dianion binding energy, eliminated by these mutations, is expressed at the wildtype Michaelis complex, and ca. 50% is only expressed at the wildtype transition state. Negative deviations from this correlation are observed when the mutation results in a decrease in enzyme reactivity at the catalytic site. The main effect of Y208T, Y208S, and Y208A mutations is to cause a reduction in the total intrinsic dianion binding energy, but the effect of Y208F extends to the catalytic site.

Steroidogenesis in plants--Biosynthesis and conversions of progesterone and other pregnane derivatives

In plants androstanes, estranes, pregnanes and corticoids have been described. Sometimes 17β-estradiol, androsterone, testosterone or progesterone were summarized as sex hormones. These steroids influence plant development: cell divisions, root and shoot growth, embryo growth, flowering, pollen tube growth and callus proliferation. First reports on the effect of applicated substances and of their endogenous occurrence date from the early twenties of the last century. This caused later on doubts on the identity of the compounds. Best investigated is the effect of progesterone. Main steps of the progesterone biosynthetic pathway have been analyzed in Digitalis. Cholesterol-side-chain-cleavage, pregnenolone and progesterone formation as well as the stereospecific reduction of progesterone are described and the corresponding enzymes are presented. Biosynthesis of androstanes, estranes and corticoids is discussed. Possible progesterone receptors and physiological reactions on progesterone application are reviewed.