Structure-based design of targeted covalent inhibitors.

Covalent inhibition is a rapidly growing discipline within drug discovery. Many historical covalent inhibitors were discovered by serendipity, with such a mechanism of action often regarded as undesirable due to potential toxicity issues. Recent progress has seen a major shift in this outlook, as covalent inhibition shows promise for targets where previous efforts to identify non-covalent small molecule inhibitors have failed. Targeted covalent inhibitors (TCIs) can offer drug discovery scientists the ability to increase the potency and/or selectivity of small molecule inhibitors, by attachment of reactive functional groups designed to covalently bind to specific sites in a target. In this tutorial review we introduce the broader concept of covalent inhibition, discuss the potential benefits and challenges of such an approach, and provide an overview of the current status of the field. We also describe some strategies and computational tools to enable successful covalent drug discovery.

Structural and Computational Insight into the Catalytic Mechanism of Limonene Epoxide Hydrolase Mutants in Stereoselective Transformations.

Directed evolution of limonene epoxide hydrolase (LEH), which catalyzes the hydrolytic desymmetrization reactions of cyclopentene oxide and cyclohexene oxide, results in (R,R)- and (S,S)-selective mutants. Their crystal structures combined with extensive theoretical computations shed light on the mechanistic intricacies of this widely used enzyme. From the computed activation energies of various pathways, we discover the underlying stereochemistry for favorable reactions. Surprisingly, some of the most enantioselective mutants that rapidly convert cyclohexene oxide do not catalyze the analogous transformation of the structurally similar cyclopentene oxide, as shown by additional X-ray structures of the variants harboring this slightly smaller substrate. We explain this puzzling observation on the basis of computational calculations which reveal a disrupted alignment between nucleophilic water and cyclopentene oxide due to the pronounced flexibility of the binding pocket. In contrast, in the stereoselective reactions of cyclohexene oxide, reactive conformations are easily reached. The unique combination of structural and computational data allows insight into mechanistic details of this epoxide hydrolase and provides guidance for future protein engineering in reactions of structurally different substrates.

Methodology Development in Directed Evolution: Exploring Options when Applying Triple-Code Saturation Mutagenesis.

Directed evolution of stereo- or regioselective enzymes as catalysts in asymmetric transformations is of particular interest in organic synthesis. Upon evolving these biocatalysts, screening is the bottleneck. To beat the numbers problem most effectively, methods and strategies for building "small but smart" mutant libraries have been developed. Herein, we compared two different strategies regarding the application of triple-code saturation mutagenesis (TCSM) at multiresidue sites of the Thermoanaerobacter brockii alcohol dehydrogenase by using distinct reduced amino-acid alphabets. By using the synthetically difficult-to-reduce prochiral ketone tetrahydrofuran-3-one as a substrate, highly R- and S-selective variants were obtained (92-99 % ee) with minimal screening. The origin of stereoselectivity was provided by molecular dynamics analyses, which is discussed in terms of the Bürgi-Dunitz trajectory.

Expanding the Armory: Predicting and Tuning Covalent Warhead Reactivity.

Targeted covalent inhibition is an established approach for increasing the potency and selectivity of potential drug candidates, as well as identifying potent and selective tool compounds for target validation studies. It is evident that identification of reversible recognition elements is essential for selective covalent inhibition, but this must also be achieved with the appropriate level of inherent reactivity of the reactive functionality (or "warhead"). Structural changes that increase or decrease warhead reactivity, guided by methods to predict the effect of those changes, have the potential to tune warhead reactivity and negate issues related to potency and/or toxicity. The half-life to adduct formation with glutathione (GSH t1/2) is a useful assay for measuring the reactivity of cysteine-targeting covalent warheads but is limited to synthesized molecules. In this manuscript we assess the ability of several experimental and computational approaches to predict GSH t1/2 for a range of cysteine targeting warheads, including a novel method based on pKa. Furthermore, matched molecular pairs analysis has been performed against our internal compound collection, revealing structure-activity relationships between a selection of different covalent warheads. These observations and methods of prediction will be valuable in the design of new covalent inhibitors with desired levels of reactivity.

Comparing Different Strategies in Directed Evolution of Enzyme Stereoselectivity: Single- versus Double-Code Saturation Mutagenesis.

Saturation mutagenesis at sites lining the binding pockets of enzymes constitutes a viable protein engineering technique for enhancing or inverting stereoselectivity. Statistical analysis shows that oversampling in the screening step (the bottleneck) increases astronomically as the number of residues in the randomization site increases, which is the reason why reduced amino acid alphabets have been employed, in addition to splitting large sites into smaller ones. Limonene epoxide hydrolase (LEH) has previously served as the experimental platform in these methodological efforts, enabling comparisons between single-code saturation mutagenesis (SCSM) and triple-code saturation mutagenesis (TCSM); these employ either only one or three amino acids, respectively, as building blocks. In this study the comparative platform is extended by exploring the efficacy of double-code saturation mutagenesis (DCSM), in which the reduced amino acid alphabet consists of two members, chosen according to the principles of rational design on the basis of structural information. The hydrolytic desymmetrization of cyclohexene oxide is used as the model reaction, with formation of either (R,R)- or (S,S)-cyclohexane-1,2-diol. DCSM proves to be clearly superior to the likewise tested SCSM, affording both R,R- and S,S-selective mutants. These variants are also good catalysts in reactions of further substrates. Docking computations reveal the basis of enantioselectivity.

Quantum Mechanics/Molecular Mechanics Modeling of Drug Metabolism: Mexiletine N-Hydroxylation by Cytochrome P450 1A2.

The mechanism of cytochrome P450(CYP)-catalyzed hydroxylation of primary amines is currently unclear and is relevant to drug metabolism; previous small model calculations have suggested two possible mechanisms: direct N-oxidation and H-abstraction/rebound. We have modeled the N-hydroxylation of (R)-mexiletine in CYP1A2 with hybrid quantum mechanics/molecular mechanics (QM/MM) methods, providing a more detailed and realistic model. Multiple reaction barriers have been calculated at the QM(B3LYP-D)/MM(CHARMM27) level for the direct N-oxidation and H-abstraction/rebound mechanisms. Our calculated barriers indicate that the direct N-oxidation mechanism is preferred and proceeds via the doublet spin state of Compound I. Molecular dynamics simulations indicate that the presence of an ordered water molecule in the active site assists in the binding of mexiletine in the active site, but this is not a prerequisite for reaction via either mechanism. Several active site residues play a role in the binding of mexiletine in the active site, including Thr124 and Phe226. This work reveals key details of the N-hydroxylation of mexiletine and further demonstrates that mechanistic studies using QM/MM methods are useful for understanding drug metabolism.

Whole-Cell-Catalyzed Multiple Regio- and Stereoselective Functionalizations in Cascade Reactions Enabled by Directed Evolution.

Biocatalytic cascade reactions using isolated stereoselective enzymes or whole cells in one-pot processes lead to value-added chiral products in a single workup. The concept has been restricted mainly to starting materials and intermediate products that are accepted by the respective wild-type enzymes. In the present study, we exploited directed evolution as a means to create E. coli whole cells for regio- and stereoselective cascade sequences that are not possible using man-made catalysts. The approach is illustrated using P450-BM3 in combination with appropriate alcohol dehydrogenases as catalysts in either two-, three-, or four-step cascade reactions starting from cyclohexane, cyclohexanol, or cyclohexanone, respectively, leading to either (R,R)-, (S,S)-, or meso-cyclohexane-1,2-diol. The one-pot conversion of cyclohexane into (R)- or (S)-2-hydroxycyclohexanone in the absence of ADH is also described.

In pursuit of an accurate spatial and temporal model of biomolecules at the atomistic level: a perspective on computer simulation.

Despite huge advances in the computational techniques available for simulating biomolecules at the quantum-mechanical, atomistic and coarse-grained levels, there is still a widespread perception amongst the experimental community that these calculations are highly specialist and are not generally applicable by researchers outside the theoretical community. In this article, the successes and limitations of biomolecular simulation and the further developments that are likely in the near future are discussed. A brief overview is also provided of the experimental biophysical methods that are commonly used to probe biomolecular structure and dynamics, and the accuracy of the information that can be obtained from each is compared with that from modelling. It is concluded that progress towards an accurate spatial and temporal model of biomacromolecules requires a combination of all of these biophysical techniques, both experimental and computational.

Reduction of α,ß-Unsaturated Ketones by Old Yellow Enzymes: Mechanistic Insights from Quantum Mechanics/Molecular Mechanics Calculations.

Enoate reductases catalyze the reduction of activated C═C bonds with high enantioselectivity. The oxidative half-reaction, which involves the addition of a hydride and a proton to opposite faces of the C═C bond, has been studied for the first time by hybrid quantum mechanics/molecular mechanics (QM/MM). The reduction of 2-cyclohexen-1-one by YqjM from Bacillus subtilis was selected as the model system. Two-dimensional QM/MM (B3LYP-D/OPLS2005) reaction pathways suggest that the hydride and proton are added as distinct steps, with the former step preceding the latter. Furthermore, we present interesting insights into the reactivity of this enzyme, including the weak binding of the substrate in the active site, the role of the two active site histidine residues for polarization of the substrate C═O bond, structural details of the transition states to hydride and proton transfer, and the role of Tyr196 as proton donor. The information presented here will be useful for the future design of enantioselective YqjM mutants for other substrates.

A multiscale approach to modelling drug metabolism by membrane-bound cytochrome P450 enzymes.

Cytochrome P450 enzymes are found in all life forms. P450s play an important role in drug metabolism, and have potential uses as biocatalysts. Human P450s are membrane-bound proteins. However, the interactions between P450s and their membrane environment are not well-understood. To date, all P450 crystal structures have been obtained from engineered proteins, from which the transmembrane helix was absent. A significant number of computational studies have been performed on P450s, but the majority of these have been performed on the solubilised forms of P450s. Here we present a multiscale approach for modelling P450s, spanning from coarse-grained and atomistic molecular dynamics simulations to reaction modelling using hybrid quantum mechanics/molecular mechanics (QM/MM) methods. To our knowledge, this is the first application of such an integrated multiscale approach to modelling of a membrane-bound enzyme. We have applied this protocol to a key human P450 involved in drug metabolism: CYP3A4. A biologically realistic model of CYP3A4, complete with its transmembrane helix and a membrane, has been constructed and characterised. The dynamics of this complex have been studied, and the oxidation of the anticoagulant R-warfarin has been modelled in the active site. Calculations have also been performed on the soluble form of the enzyme in aqueous solution. Important differences are observed between the membrane and solution systems, most notably for the gating residues and channels that control access to the active site. The protocol that we describe here is applicable to other membrane-bound enzymes.

Biocatalytic route to chiral acyloins: P450-catalyzed regio- and enantioselective α-hydroxylation of ketones.

P450-BM3 and mutants of this monooxygenase generated by directed evolution are excellent catalysts for the oxidative α-hydroxylation of ketones with formation of chiral acyloins with high regioselectivity (up to 99%) and enantioselectivity (up to 99% ee). This constitutes a new route to a class of chiral compounds that are useful intermediates in the synthesis of many kinds of biologically active compounds.

Reshaping an enzyme binding pocket for enhanced and inverted stereoselectivity: use of smallest amino acid alphabets in directed evolution.

Directed evolution based on saturation mutagenesis at sites lining the binding pocket is a commonly practiced strategy for enhancing or inverting the stereoselectivity of enzymes for use in organic chemistry or biotechnology. However, as the number of residues in a randomization site increases to five or more, the screening effort for 95 % library coverage increases astronomically until it is no longer feasible. We propose the use of a single amino acid for saturation mutagenesis at superlarge randomization sites comprising 10 or more residues. When used to reshape the binding pocket of limonene epoxide hydrolase, this strategy, which drastically reduces the search space and thus the screening effort, resulted in R,R- and S,S-selective mutants for the hydrolytic desymmetrization of cyclohexene oxide and other epoxides. X-ray crystal structures and docking studies of the mutants unveiled the source of stereoselectivity and shed light on the mechanistic intricacies of this enzyme.

Bioorthogonal Enzymatic Activation of Caged Compounds.

Engineered cytochrome P450 monooxygenase variants are reported as highly active and selective catalysts for the bioorthogonal uncaging of propargylic and benzylic ether protected substrates, including uncaging in living E. coli. observed selectivity is supported by induced-fit docking and molecular dynamics simulations. This proof-of-principle study points towards the utility of bioorthogonal enzyme/protecting group pairs for applications in the life sciences.

Extreme synergistic mutational effects in the directed evolution of a baeyer-villiger monooxygenase as catalyst for asymmetric sulfoxidation.

Structure-based directed evolution utilizing iterative saturation mutagenesis (ISM) has been applied to phenyl acetone monooxygenase (PAMO), a thermally robust Baeyer-Villiger monooxygenase, in the quest to access a mutant which displays reversed enantioselectivity in the asymmetric sulfoxidation of prochiral thioethers. Whereas WT PAMO leads to 90% ee in the sulfoxidation of p-methylbenzyl methyl thioether with preference for the (S)-sulfoxide, the evolved mutant I67Q/P440F/A442N/L443I is 95% (R)-selective in the reaction of this and other thioethers. Partial deconvolution of the (R)-selective mutant with generation of the respective four single mutants shows that all of them are (S)-selective, which points to pronounced synergism (cooperative nonadditivity) when they interact in concert. Complete deconvolution with formation of all combinatorial forms of the respective double and triple mutants allows the designed construction of a fitness landscape featuring all 24 upward pathways leading from WT to the (R)-selective quadruple mutant. In all 24 trajectories strong cooperative mutational effects were found as well, which indicates that such mutational changes in enzymes constitute nonlinear systems. A theoretical analysis based on induced fit docking explains many of the observed effects on a molecular level.

Quantum mechanics/molecular mechanics modeling of regioselectivity of drug metabolism in cytochrome P450 2C9.

Cytochrome P450 enzymes (P450s) are important in drug metabolism and have been linked to adverse drug reactions. P450s display broad substrate reactivity, and prediction of metabolites is complex. QM/MM studies of P450 reactivity have provided insight into important details of the reaction mechanisms and have the potential to make predictions of metabolite formation. Here we present a comprehensive study of the oxidation of three widely used pharmaceutical compounds (S-ibuprofen, diclofenac, and S-warfarin) by one of the major drug-metabolizing P450 isoforms, CYP2C9. The reaction barriers to substrate oxidation by the iron-oxo species (Compound I) have been calculated at the B3LYP-D/CHARMM27 level for different possible metabolism sites for each drug, on multiple pathways. In the cases of ibuprofen and warfarin, the process with the lowest activation energy is consistent with the experimentally preferred metabolite. For diclofenac, the pathway leading to the experimentally observed metabolite is not the one with the lowest activation energy. This apparent inconsistency with experiment might be explained by the two very different binding modes involved in oxidation at the two competing positions. The carboxylate of diclofenac interacts strongly with the CYP2C9 Arg108 side chain in the transition state for formation of the observed metabolite-but not in that for the competing pathway. We compare reaction barriers calculated both in the presence and in the absence of the protein and observe a marked improvement in selectivity prediction ability upon inclusion of the protein for all of the substrates studied. The barriers calculated with the protein are generally higher than those calculated in the gas phase. This suggests that active-site residues surrounding the substrate play an important role in controlling selectivity in CYP2C9. The results show that inclusion of sampling (particularly) and dispersion effects is important in making accurate predictions of drug metabolism selectivity of P450s using QM/MM methods.

A practical guide to modelling enzyme-catalysed reactions.

Molecular modelling and simulation methods are increasingly at the forefront of elucidating mechanisms of enzyme-catalysed reactions, and shedding light on the determinants of specificity and efficiency of catalysis. These methods have the potential to assist in drug discovery and the design of novel protein catalysts. This Tutorial Review highlights some of the most widely used modelling methods and some successful applications. Modelling protocols commonly applied in studying enzyme-catalysed reactions are outlined here, and some practical implications are considered, with cytochrome P450 enzymes used as a specific example.

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling.

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.

Does compound I vary significantly between isoforms of cytochrome P450?

The cytochrome P450 (CYP) enzymes are important in many areas, including pharmaceutical development. Subtle changes in the electronic structure of the active species, Compound I, have been postulated previously to account partly for the experimentally observed differences in reactivity between isoforms. Current predictive models of CYP metabolism typically assume an identical Compound I in all isoforms. Here we present a method to calculate the electronic structure and to estimate the Fe-O bond enthalpy of Compound I, and apply it to several human and bacterial CYP isoforms. Conformational flexibility is accounted for by sampling large numbers of structures from molecular dynamics simulations, which are subsequently optimized with density functional theory (B3LYP) based quantum mechanics/molecular mechanics. The observed differences in Compound I between human isoforms are small: They are generally smaller than the spread of values obtained for the same isoform starting from different initial structures. Hence, it is unlikely that the variation in activity between human isoforms is due to differences in the electronic structure of Compound I. A larger difference in electronic structure is observed between the human isoforms and P450(cam) and may be explained by the slightly different hydrogen-bonding environment surrounding the cysteinyl sulfur. The presence of substrate in the active site of all isoforms studied appears to cause a slight decrease in the Fe-O bond enthalpy, apparently due to displacement of water out of the active site, suggesting that Compound I is less stable in the presence of substrate.

Comment on "A stationary-wave model of enzyme catalysis" by Carlo Canepa.

Energy barriers for enzyme-catalyzed reactions calculated with quantum mechanics/molecular mechanics (QM/MM) and empirical valence bond (EVB) methods can be in excellent agreement with activation energies derived from experiments, supporting the applicability of transition state theory for enzymic reactions.

Correction: Simultaneous engineering of an enzyme's entrance tunnel and active site: the case of monoamine oxidase MAO-N.

[This corrects the article DOI: 10.1039/C6SC05381E.].