Dušan Petrović

Shuffling Active Site Substate Populations Affects Catalytic Activity: The Case of Glucose Oxidase

Glucose oxidase has wide applications in the pharmaceutical, chemical, and food industries. Many recent studies have enhanced key properties of this enzyme using directed evolution, yet without being able to reveal why these mutations are actually beneficial. This work presents a synergistic combination of experimental and computational methods, indicating how mutations, even when distant from the active site, positively affect glucose oxidase catalysis. We have determined the crystal structures of glucose oxidase mutants containing molecular oxygen in the active site. The catalytically important His516 residue has been previously shown to be flexible in the wild-type enzyme. The molecular dynamics simulations performed in this work allow us to quantify this floppiness, revealing that His516 exists in two states: catalytic and noncatalytic. The relative populations of these two substates are almost identical in the wild-type enzyme, with His516 readily shuffling between them. In the glucose oxidase mutants, on the other hand, the mutations enrich the catalytic His516 conformation and reduce the flexibility of this residue, leading to an enhancement in their catalytic efficiency. This study stresses the benefit of active site preorganization with respect to enzyme conversion rates by reducing molecular reorientation needs. We further suggest that the computational approach based on Hamiltonian replica exchange molecular dynamics, used in this study, may be a general approach to screening in silico for improved enzyme variants involving flexible catalytic residues.

Enzyme Architecture: Modeling the Operation of a Hydrophobic Clamp in Catalysis by Triosephosphate Isomerase

Triosephosphate isomerase (TIM) is a proficient catalyst of the reversible isomerization of dihydroxyacetone phosphate (DHAP) to d-glyceraldehyde phosphate (GAP), via general base catalysis by E165. Historically, this enzyme has been an extremely important model system for understanding the fundamentals of biological catalysis. TIM is activated through an energetically demanding conformational change, which helps position the side chains of two key hydrophobic residues (I170 and L230), over the carboxylate side chain of E165. This is critical both for creating a hydrophobic pocket for the catalytic base and for maintaining correct active site architecture. Truncation of these residues to alanine causes significant falloffs in TIM’s catalytic activity, but experiments have failed to provide a full description of the action of this clamp in promoting substrate deprotonation. We perform here detailed empirical valence bond calculations of the TIM-catalyzed deprotonation of DHAP and GAP by both wild-type TIM and its I170A, L230A, and I170A/L230A mutants, obtaining exceptional quantitative agreement with experiment. Our calculations provide a linear free energy relationship, with slope 0.8, between the activation barriers and Gibbs free energies for these TIM-catalyzed reactions. We conclude that these clamping side chains minimize the Gibbs free energy for substrate deprotonation, and that the effects on reaction driving force are largely expressed at the transition state for proton transfer. Our combined analysis of previous experimental and current computational results allows us to provide an overview of the breakdown of ground-state and transition state effects in enzyme catalysis in unprecedented detail, providing a molecular description of the operation of a hydrophobic clamp in triosephosphate isomerase.

Conformational Sampling of the Intrinsically Disordered C-terminal Tail of DERA is Important for Enzyme Catalysis

2-Deoxyribose-5-phosphate aldolase (DERA) catalyzes the reversible conversion of acetaldehyde and glyceraldehyde-3-phosphate into deoxyribose-5-phosphate. DERA is used as a biocatalyst for the synthesis of drugs such as statins and is a promising pharmaceutical target due to its involvement in nucleotide catabolism. Despite previous biochemical studies suggesting the catalytic importance of the C-terminal tyrosine residue found in several bacterial DERAs, the structural and functional basis of its participation in catalysis remains elusive because the electron density for the last eight to nine residues (i.e., the C-terminal tail) is absent in all available crystal structures. Using a combination of NMR spectroscopy and molecular dynamics simulations, we conclusively show that the rarely studied C-terminal tail of E. coli DERA (ecDERA) is intrinsically disordered and exists in equilibrium between open and catalytically relevant closed states, where the C-terminal tyrosine (Y259) enters the active site. Nuclear Overhauser effect distance restraints, obtained due to the presence of a substantial closed state population, were used to derive the solution-state structure of the ecDERA closed state. Real-time NMR hydrogen/deuterium exchange experiments reveal that Y259 is required for efficiency of the proton abstraction step of the catalytic reaction. Phosphate titration experiments show that, in addition to the phosphate-binding residues located near the active site, as observed in the available crystal structures, ecDERA contains previously unknown auxiliary phosphate-binding residues on the C-terminal tail which could facilitate in orienting Y259 in an optimal position for catalysis. Thus, we present significant insights into the structural and mechanistic importance of the ecDERA C-terminal tail and illustrate the role of conformational sampling in enzyme catalysis.

Simulation-Guided Design of Cytochrome P450 for Chemo- and Regioselective Macrocyclic Oxidation

Engineering high chemo-, regio-, and stereoselectivity is a prerequisite for enzyme usage in organic synthesis. Cytochromes P450 can oxidize a broad range of substrates, including macrocycles, which are becoming popular scaffolds for therapeutic agents. However, a large conformational space explored by macrocycles not only reduces the selectivity of oxidation but also impairs computational enzyme design strategies based on docking and molecular dynamics (MD) simulations. We present a novel design workflow that uses enhanced-sampling Hamiltonian replica exchange (HREX) MD and focuses on quantifying the substrate binding for suggesting the mutations to be made. This computational approach is applied to P450 BM3 with the aim to shift regioselectively toward one of the numerous possible positions during β-cembrenediol oxidation. The predictions are experimentally tested and the resulting product distributions validate our design strategy, as single mutations led up to 5-fold regioselectivity increases. We thus conclude that the HREX-MD-based workflow is a promising tool for the identification of positions for mutagenesis aiming at P450 enzymes with improved regioselectivity.

One-Pot, Two-Step Hydroxylation of the Macrocyclic Diterpenoid β-Cembrenediol Catalyzed by P450 BM3 Mutants

Cytochrome P450 monooxygenases (P450s) are involved in the biosynthesis of a wide range of bioactive secondary metabolites. They often introduce several oxy functionalities at different positions of a substrate through multiple steps and produce a range of oxidized derivatives. Herein, we describe a one‐pot two‐step hydroxylation of the diterpenoid β‐cembrenediol isolated from the plant Nicotiana tabacum. This 14‐membered macrocycle shows neuroprotective effects and is, along with its oxidized derivatives, of pharmaceutical interest. Sequential hydroxylations catalyzed by the regioselective P450 BM3 mutants F87A/I263L and V78A/F87G yielded the epimeric (9S,10R/S)‐β‐cembrenetetraols with a diastereomeric ratio of 48:52. The replacement of the mutant V78A/F87G with L75A/V78A/F87G in the second step improves the diastereomeric ratio up to 10:90. Absolute configurations of the newly introduced hydroxy groups were determined by quantum‐mechanical calculations of NMR spectra.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

How accurately do force fields represent protein side chain ensembles

Although the protein backbone is the most fundamental part of the structure, the fine‐tuning of side‐chain conformations is important for protein function, for example, in protein‐protein and protein‐ligand interactions, and also in enzyme catalysis. While several benchmarks testing the performance of protein force fields for side chain properties have already been published, they often considered only a few force fields and were not tested against the same experimental observables; hence, they are not directly comparable. In this work, we explore the ability of twelve force fields, which are different flavors of AMBER, CHARMM, OPLS, or GROMOS, to reproduce average rotamer angles and rotamer populations obtained from extensive NMR studies of the 3J and residual dipolar coupling constants for two small proteins: ubiquitin and GB3. Based on a total of 196u2009μs sampling time, our results reveal that all force fields identify the correct side chain angles, while the AMBER and CHARMM force fields clearly outperform the OPLS and GROMOS force fields in estimating rotamer populations. The three best force fields for representing the protein side chain dynamics are AMBER 14SB, AMBER 99SB*‐ILDN, and CHARMM36. Furthermore, we observe that the side chain ensembles of buried amino acid residues are generally more accurately represented than those of the surface exposed residues.

The evolution of multiple active site configurations in a designed enzyme

Developments in computational chemistry, bioinformatics, and laboratory evolution have facilitated the de novo design and catalytic optimization of enzymes. Besides creating useful catalysts, the generation and iterative improvement of designed enzymes can provide valuable insight into the interplay between the many phenomena that have been suggested to contribute to catalysis. In this work, we follow changes in conformational sampling, electrostatic preorganization, and quantum tunneling along the evolutionary trajectory of a designed Kemp eliminase. We observe that in the Kemp Eliminase KE07, instability of the designed active site leads to the emergence of two additional active site configurations. Evolutionary conformational selection then gradually stabilizes the most efficient configuration, leading to an improved enzyme. This work exemplifies the link between conformational plasticity and evolvability and demonstrates that residues remote from the active sites of enzymes play crucial roles in controlling and shaping the active site for efficient catalysis.Generation and iterative optimization of designed enzymes can provide valuable insights for a more efficient catalysis. Here the authors have followed the iterative improvement of a designed Kemp eliminase and show that remote point mutations could remodel the designed active site via substantial conformational reorganization.

Loop Motion in Triosephosphate Isomerase is not a Simple Open and Shut Case

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 grips the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple open and shut case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Empirical Valence Bond Simulations Suggest a Direct Hydride Transfer Mechanism for Human Diamine Oxidase

Diamine oxidase (DAO) is an enzyme involved in the regulation of cell proliferation and the immune response. This enzyme performs oxidative deamination in the catabolism of biogenic amines, including, among others, histamine, putrescine, spermidine, and spermine. The mechanistic details underlying the reductive half-reaction of the DAO-catalyzed oxidative deamination which leads to the reduced enzyme cofactor and the aldehyde product are, however, still under debate. The catalytic mechanism was proposed to involve a prototropic shift from the substrate–Schiff base to the product–Schiff base, which includes the rate-limiting cleavage of the Cα–H bond by the conserved catalytic aspartate. Our detailed mechanistic study, performed using a combined quantum chemical cluster approach with empirical valence bond simulations, suggests that the rate-limiting cleavage of the Cα–H bond involves direct hydride transfer to the topaquinone cofactor—a mechanism that does not involve the formation of a Schiff base. Additional investigation of the D373E and D373N variants supported the hypothesis that the conserved catalytic aspartate is indeed essential for the reaction; however, it does not appear to serve as the catalytic base, as previously suggested. Rather, the electrostatic contributions of the most significant residues (including D373), together with the proximity of the Cu2+ cation to the reaction site, lower the activation barrier to drive the chemical reaction.