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.

SilE-R and SilE-S-DABB Proteins Catalying Enantiospecific Hydrolysis of Organosilyl Ethers

Silyl ethers fulfil a fundamental role in synthetic organic chemistry as protecting groups and their selective cleavage is an important factor in their application. We present here for the first time two enzymes, SilE-R and SilE-S, which are able to hydrolyse silyl ethers. They belong to the stress-response dimeric A/B barrel domain (DABB) family and are able to cleave the Si-O bond with opposite enantiopreference. Silyl ethers containing aromatic, cyclic or aliphatic alcohols and, depending on the alcohol moiety, silyl functions as large as TBDMS are accepted. The X-ray crystal structure of SilE-R, determined to a resolution of 1.98 Å, in combination with mutational studies, revealed an active site featuring two histidine residues, H8 and H79, which likely act synergistically as nucleophile and Brønsted base in the hydrolytic mechanism, which has not previously been described for enzymes. Although the natural function of SilE-R and SilE-S is unknown, we propose that these 'silyl etherases' may have significant potential for synthetic applications.

Assessing the effect of a series of mutations on the dynamic behavior of phosphite dehydrogenase using molecular docking, molecular dynamics and quantum mechanics/molecular mechanics simulations

Discovered in Pseudomonas stutzeri, phosphite dehydrogenase (PTDH) is an enzyme that catalyzes the oxidation of phosphite to phosphate while simultaneously reducing NAD⁺ to NADH. Despite several investigations into the mechanism of reaction and cofactor regeneration, only a few studies have focused on improving the activity and stability of PTDH. In this study, we combine molecular docking, molecular dynamics (MD) simulation, and Quantum Mechanics/Molecular Mechanics (QM/MM) to identify the impact of 30 mutations on the activity and stability of PTDH. Molecular docking results suggest that E266Q, K76A, K76M, K76R, K76C, and R237K can act on the NAD⁺ binding site through relatively weak bond development due to their high free binding energy. Moreover, Mulliken population analysis and potential energy barrier indicate that T101A, E175A, E175A/A176R, A176R, and E266Q act on phosphite oxidation. The mutants M53N, M53A, K76R, D79N, D79A, T101A, W134A, W134F Y139F, A146S, E175A, F198I, F198M, E266Q, H292K, S295A, R301K, and R301A were found to act on the structural dynamic of PTDH. The remaining mutants cause the loss of the nitrogen atom of R237 and H292, respectively, inactivating the enzyme. This study provides specific explanations of how mutations affect weak interactions of PTDH. The results should allow researchers to conduct experimental studies to improve PTDH activity and stability.Communicated by Ramaswamy H. Sarma.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

Recent Advances in Biocatalysis with Chemical Modification and Expanded Amino Acid Alphabet

The two main strategies for enzyme engineering, directed evolution and rational design, have found widespread applications in improving the intrinsic activities of proteins. Although numerous advances have been achieved using these ground-breaking methods, the limited chemical diversity of the biopolymers, restricted to the 20 canonical amino acids, hampers creation of novel enzymes that Nature has never made thus far. To address this, much research has been devoted to expanding the protein sequence space via chemical modifications and/or incorporation of noncanonical amino acids (ncAAs). This review provides a balanced discussion and critical evaluation of the applications, recent advances, and technical breakthroughs in biocatalysis for three approaches: (i) chemical modification of cAAs, (ii) incorporation of ncAAs, and (iii) chemical modification of incorporated ncAAs. Furthermore, the applications of these approaches and the result on the functional properties and mechanistic study of the enzymes are extensively reviewed. We also discuss the design of artificial enzymes and directed evolution strategies for enzymes with ncAAs incorporated. Finally, we discuss the current challenges and future perspectives for biocatalysis using the expanded amino acid alphabet.

Participation of S and Se in hydrogen and chalcogen bonds

The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.

Examining the Mechanism of Phosphite Dehydrogenase with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

The projected decline of available phosphorus necessitates alternative methods to derive usable phosphate for fertilizer and other applications. Phosphite dehydrogenase oxidizes phosphite to phosphate with the cofactor NAD⁺ serving as the hydride acceptor. In addition to producing phosphate, this enzyme plays an important role in NADH cofactor regeneration processes. Mixed quantum mechanical/molecular mechanical free energy simulations were performed to elucidate the mechanism of this enzyme and to identify the protonation states of the substrate and product. Specifically, the finite temperature string method with umbrella sampling was used to generate the free energy surfaces and determine the minimum free energy paths for six different initial conditions that varied in the protonation state of the substrate and the position of the nucleophilic water molecule. In contrast to previous studies, the mechanism predicted by all six independent strings is a concerted but asynchronous dissociative mechanism in which hydride transfer from the phosphite substrate to NAD⁺ occurs prior to attack by the nucleophilic water molecule. His292 is identified as the most likely general base that deprotonates the attacking water molecule. However, Arg237 could also serve as this base if it were deprotonated and His292 were protonated prior to the main chemical transformation, although this scenario is less probable. The simulations indicate that the phosphite substrate is monoanionic in its active form and that the most likely product is dihydrogen phosphate. These mechanistic insights may be helpful for designing mutant enzymes or artificial constructs that convert phosphite to phosphate and NAD⁺ to NADH more effectively.

π-Hole Interactions Involving Nitro Aromatic Ligands in Protein Structures

Studying noncanonical intermolecular interactions between a ligand and a protein constitutes an emerging research field. Identifying synthetically accessible molecular fragments that can engage in intermolecular interactions is a key objective in this area. Here, it is shown that so-called "π-hole interactions" are present between the nitro moiety in nitro aromatic ligands and lone pairs within protein structures (water and protein carbonyls and sulfurs). Ample structural evidence was found in a PDB analysis and computations reveal interaction energies of about -5 kcal mol-1 for ligand-protein π-hole interactions. Several examples are highlighted for which a π-hole interaction is implicated in the superior binding affinity or inhibition of a nitro aromatic ligand versus a similar non-nitro analogue. The discovery that π-hole interactions with nitro aromatics are significant within protein structures parallels the finding that halogen bonds are biologically relevant. This has implications for the interpretation of ligand-protein complexation phenomena, for example, involving the more than 50 approved drugs that contain a nitro aromatic moiety.

Temperature-Independent Kinetic Isotope Effects as Evidence for a Marcus-like Model of Hydride Tunneling in Phosphite Dehydrogenase

Phosphite dehydrogenase catalyzes the transfer of a hydride from phosphite to NAD⁺, producing phosphate and NADH. We have evaluated the role of hydride tunneling in a thermostable variant of this enzyme (17X-PTDH) by measuring the temperature dependence of the primary ²H kinetic isotope effects (KIEs) between 5 and 45 °C. Pre-steady-state kinetic measurements were used to demonstrate that the hydride transfer is rate-determining across this temperature range and that the observed KIEs are equal to the intrinsic isotope effect on the chemical step. The KIEs on the pre-exponential factor (A_H/A_D) and the activation energy (ΔE_a) were 1.6 ± 0.1 and 0.21 ± 0.05 kcal/mol, respectively, suggesting that 17X-PTDH facilitates extensive tunneling of both isotopes via a Marcus-like model. Site-directed mutagenesis was used to evaluate the role of an active site threonine (Thr104) found on the back face of the nicotinamide in promoting the close packing of the substrates. In mutants with reduced steric bulk at this position, values of A_H/A_D and ΔE_a fall within the range describing semiclassical "over the barrier" reactivity, suggesting that Thr104 acts as a steric backstop to promote tunneling in 17X-PTDH. Whereas hydrogen tunneling is now a widely appreciated feature of C-H activating enzymes, these observations with a P-H activating system are consistent with the proposal that tunneling is likely to be a common feature on all enzymes that catalyze hydrogen transfers.

Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.

Intramolecular π–hole interactions with nitro aromatics

A thorough CSD evaluation and DFT calculations reveal that intramolecular ElR⋯NO₂ π-hole interactions can stabilize conformers of nitro aromatics by up to 2–3 kcal mol⁻¹ when ElR and N separated by ≥4 bonds.

18O Kinetic Isotope Effects Reveal an Associative Transition State for Phosphite Dehydrogenase Catalyzed Phosphoryl Transfer

Phosphite dehydrogenase (PTDH) catalyzes an unusual phosphoryl transfer reaction in which water displaces a hydride leaving group. Despite extensive effort, it remains unclear whether PTDH catalysis proceeds via an associative or dissociative mechanism. Here, primary ²H and secondary ¹⁸O kinetic isotope effects (KIEs) were determined and used together with computation to characterize the transition state (TS) catalyzed by a thermostable PTDH (17X-PTDH). The large, normal ¹⁸O KIEs suggest an associative mechanism. Various transition state structures were computed within a model of the enzyme active site and ²H and ¹⁸O KIEs were predicted to evaluate the accuracy of each TS. This analysis suggests that 17X-PTDH catalyzes an associative process with little leaving group displacement and extensive nucleophilic participation. This tight TS is likely a consequence of the extremely poor leaving group requiring significant P-O bond formation to expel the hydride. This finding contrasts with the dissociative TSs in most phosphoryl transfer reactions from phosphate mono- and diesters.

A Multiscale Simulation Approach to Modeling Drug-Protein Binding Kinetics

Drug-target binding kinetics has recently emerged as a sometimes critical determinant of in vivo efficacy and toxicity. Its rational optimization to improve potency or reduce side effects of drugs is, however, extremely difficult. Molecular simulations can play a crucial role in identifying features and properties of small ligands and their protein targets affecting the binding kinetics, but significant challenges include the long time scales involved in (un)binding events and the limited accuracy of empirical atomistic force fields (lacking, e.g., changes in electronic polarization). In an effort to overcome these hurdles, we propose a method that combines state-of-the-art enhanced sampling simulations and quantum mechanics/molecular mechanics (QM/MM) calculations at the BLYP/VDZ level to compute association free energy profiles and characterize the binding kinetics in terms of structure and dynamics of the transition state ensemble. We test our combined approach on the binding of the anticancer drug Imatinib to Src kinase, a well-characterized target for cancer therapy with a complex binding mechanism involving significant conformational changes. The results indicate significant changes in polarization along the binding pathways, which affect the predicted binding kinetics. This is likely to be of widespread importance in binding of ligands to protein targets.

Intermolecular π-hole/n→π* interactions with carbon monoxide ligands in crystal structures

A thorough analysis of the Cambridge Structure Database reveals that intermolecular π-hole/n→π* interactions with carbon monoxide ligands are abundant in the solid state and somewhat directional, particularly with fac-like M(CO)3 fragments (P < 4.0). High level DFT calculations suggest interacting energies up to about -10 kcal mol-1 for adducts of charge neutral complexes.

π-Hole/n→π* interactions with acetonitrile in crystal structures

A thorough analysis of the Cambridge Structure Database shows that π-hole/n→π* interactions with acetonitrile are abundant in the solid state, particularly when acetonitrile is coordinated to a metal. The interaction is weakly directional (P ≤ 1.5) and high level computations suggest a complexation energy of about -5 kcal mol-1.

QM/MM simulations identify the determinants of catalytic activity differences between type II dehydroquinase enzymes

Type II dehydroquinase enzymes (DHQ2), recognized targets for antibiotic drug discovery, show significantly different activities dependent on the species: DHQ2 from Mycobacterium tuberculosis (MtDHQ2) and Helicobacter pylori (HpDHQ2) show a 50-fold difference in catalytic efficiency. Revealing the determinants of this activity difference is important for our understanding of biological catalysis and further offers the potential to contribute to tailoring specificity in drug design. Molecular dynamics simulations using a quantum mechanics/molecular mechanics potential, with correlated ab initio single point corrections, identify and quantify the subtle determinants of the experimentally observed difference in efficiency. The rate-determining step involves the formation of an enolate intermediate: more efficient stabilization of the enolate and transition state of the key step in MtDHQ2, mainly by the essential residues Tyr24 and Arg19, makes it more efficient than HpDHQ2. Further, a water molecule, which is absent in MtDHQ2 but involved in generation of the catalytic Tyr22 tyrosinate in HpDHQ2, was found to destabilize both the transition state and the enolate intermediate. The quantification of the contribution of key residues and water molecules in the rate-determining step of the mechanism also leads to improved understanding of higher potencies and specificity of known inhibitors, which should aid ongoing inhibitor design.

Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures

Drug action is inherently multiscale: it connects molecular interactions to emergent properties at cellular and larger scales. Simulation techniques at each of these different scales are already central to drug design and development, but methods capable of connecting across these scales will extend understanding of complex mechanisms and the ability to predict biological effects. Improved algorithms, ever-more-powerful computing architectures and the accelerating growth of rich datasets are driving advances in multiscale modeling methods capable of bridging chemical and biological complexity from the atom to the cell. Particularly exciting is the development of highly detailed, structure-based, physical simulations of biochemical systems, which are now able to access experimentally relevant timescales for large systems and, at the same time, achieve unprecedented accuracy. In this Perspective, we discuss how emerging data-rich, physics-based multiscale approaches are of the cusp of realizing long-promised impact in the discovery, design and development of novel therapeutics. We highlight emerging methods and applications in this growing field, and outline how different scales can be combined in practical modelling and simulation strategies.

Chemical Control in the Battle against Fidelity in Promiscuous Natural Product Biosynthesis: The Case of Trichodiene Synthase

Terpene cyclases catalyze the highly stereospecific molding of polyisoprenes into terpenes, which are precursors to most known natural compounds. The isoprenoids are formed via intricate chemical cascades employing rich, yet highly erratic, carbocation chemistry. It is currently not well understood how these biocatalysts achieve chemical control. Here, we illustrate the catalytic control exerted by trichodiene synthase, and in particular, we discover two features that could be general catalytic tools adopted by other terpenoid cyclases. First, to avoid formation of byproducts, the enzyme raises the energy of bisabolyl carbocation, which is a general mechanistic branching point in many sesquiterpene cyclases, resulting in an essentially concerted cyclization cascade. Second, we identify a sulfur-carbocation dative bonding interaction that anchors the bisabolyl cation in a reactive conformation, avoiding tumbling and premature deprotonation. Specifically, Met73 acts as a chameleon, shifting from an initial sulfur-π interaction in the Michaelis complex to a sulfur-carbocation complex during catalysis.

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.

N-H···N Hydrogen Bonds Involving Histidine Imidazole Nitrogen Atoms: A New Structural Role for Histidine Residues in Proteins

The amino acid histidine can play a significant role in the structure and function of proteins. Its various functions include enzyme catalysis, metal binding activity, and involvement in cation-π, π-π, salt-bridge, and other types of noncovalent interactions. Although histidine's imidazole nitrogens (Nδ and Nε) are known to participate in hydrogen bond (HB) interactions as an acceptor or a donor, a systematic study of N-H···N HBs with the Nδ/Nε atom as the acceptor has not been conducted. In this study, we have examined two data sets of ultra-high-resolution (data set I) and very high-resolution (data set II) protein structures and identified 28 and 4017 examples of HBs of the N-H···Nδ/Nε type from both data sets involving histidine imidazole nitrogen as the acceptor. In nearly 70% of them, the main-chain N-H bond is the HB donor, and a majority of the examples are from the N-H group separated by two residues (Ni+2-Hi+2) from histidine. Quantum chemical calculations using model compounds were performed with imidazole and N-methylacetamide, and they assumed conformations from 19 examples from data set I with N-H···Nδ/Nε HBs. Basis set superposition error-corrected interaction energies varied from -5.0 to -6.78 kcal/mol. We also found that the imidazole nitrogen of 9% of histidine residues forming N-H···Nδ/Nε interactions in data set II participate in bifurcated HBs. Natural bond orbital analyses of model compounds indicate that the strength of each HB is mutually influenced by the other. Histidine residues involved in Ni+2-Hi+2···Nδi/Nεi HBs are frequently observed in a specific N-terminal capping position giving rise to a novel helix-capping motif. Along with their predominant occurrence in loop segments, we propose a new structural role for histidines in protein structures.

Towards design strategies for anion–π interactions in crystal engineering

This highlight article summarizes some of the fundamental aspects of the anion–π interaction leading to several design strategies for generating it in solids. In the main body we highlight some relevant examples that illustrate the viability of these strategies and the importance of anion–π interactions in crystal engineering.

Directionality of π-holes in nitro compounds

A statistical survey of the Cambridge Structural Database reveals that the interaction between the π-holes of nitro groups and electron-rich atoms is somewhat directional. High-level ab initio computations indicate energies up to -6.6 kcal mol(-1).