A combinatorially complete epistatic fitness landscape in an enzyme active site.

Protein engineering often targets binding pockets or active sites which are enriched in epistasis-nonadditive interactions between amino acid substitutions-and where the combined effects of multiple single substitutions are difficult to predict. Few existing sequence-fitness datasets capture epistasis at large scale, especially for enzyme catalysis, limiting the development and assessment of model-guided enzyme engineering approaches. We present here a combinatorially complete, 160,000-variant fitness landscape across four residues in the active site of an enzyme. Assaying the native reaction of a thermostable ß-subunit of tryptophan synthase (TrpB) in a nonnative environment yielded a landscape characterized by significant epistasis and many local optima. These effects prevent simulated directed evolution approaches from efficiently reaching the global optimum. There is nonetheless wide variability in the effectiveness of different directed evolution approaches, which together provide experimental benchmarks for computational and machine learning workflows. The most-fit TrpB variants contain a substitution that is nearly absent in natural TrpB sequences-a result that conservation-based predictions would not capture. Thus, although fitness prediction using evolutionary data can enrich in more-active variants, these approaches struggle to identify and differentiate among the most-active variants, even for this near-native function. Overall, this work presents a large-scale testing ground for model-guided enzyme engineering and suggests that efficient navigation of epistatic fitness landscapes can be improved by advances in both machine learning and physical modeling.

The ß-subunit of tryptophan synthase is a latent tyrosine synthase.

Aromatic amino acids and their derivatives are diverse primary and secondary metabolites with critical roles in protein synthesis, cell structure and integrity, defense and signaling. All de novo aromatic amino acid production relies on a set of ancient and highly conserved chemistries. Here we introduce a new enzymatic transformation for L-tyrosine synthesis by demonstrating that the ß-subunit of tryptophan synthase-which natively couples indole and L-serine to form L-tryptophan-can act as a latent 'tyrosine synthase'. A single substitution of a near-universally conserved catalytic residue unlocks activity toward simple phenol analogs and yields exclusive para carbon-carbon bond formation to furnish L-tyrosines. Structural and mechanistic studies show how a new active-site water molecule orients phenols for a nonnative mechanism of alkylation, with additional directed evolution resulting in a net >30,000-fold rate enhancement. This new biocatalyst can be used to efficiently prepare valuable L-tyrosine analogs at gram scales and provides the missing chemistry for a conceptually different pathway to L-tyrosine.

Iron Heme Enzyme-Catalyzed Cyclopropanations with Diazirines as Carbene Precursors: Computational Explorations of Diazirine Activation and Cyclopropanation Mechanism.

The mechanism of cyclopropanations with diazirines as air-stable and user-friendly alternatives to commonly employed diazo compounds within iron heme enzyme-catalyzed carbene transfer reactions has been studied by means of density functional theory (DFT) calculations of model systems, quantum mechanics/molecular mechanics (QM/MM) calculations, and molecular dynamics (MD) simulations of the iron carbene and the cyclopropanation transition state in the enzyme active site. The reaction is initiated by a direct diazirine-diazo isomerization occurring in the active site of the enzyme. In contrast, an isomerization mechanism proceeding via the formation of a free carbene intermediate in lieu of a direct, one-step isomerization process was observed for model systems. Subsequent reaction with benzyl acrylate takes place through stepwise C-C bond formation via a diradical intermediate, delivering the cyclopropane product. The origin of the observed diastereo- and enantioselectivity in the enzyme was investigated through MD simulations, which indicate a preferred formation of the cis-cyclopropane by steric control.

MicroED Structure of a Protoglobin Reactive Carbene Intermediate.

Microcrystal electron diffraction (MicroED) is an emerging technique that has shown great potential for describing new chemical and biological molecular structures. Several important structures of small molecules, natural products, and peptides have been determined using ab initio methods. However, only a couple of novel protein structures have thus far been derived by MicroED. Taking advantage of recent technological advances, including higher acceleration voltage and using a low-noise detector in counting mode, we have determined the first structure of an Aeropyrum pernix protoglobin (ApePgb) variant by MicroED using an AlphaFold2 model for phasing. The structure revealed that mutations introduced during directed evolution enhance carbene transfer activity by reorienting an α helix of ApePgb into a dynamic loop, making the catalytic active site more readily accessible. After exposing the tiny crystals to the substrate, we also trapped the reactive iron-carbenoid intermediate involved in this engineered ApePgb's new-to-nature activity, a challenging carbene transfer from a diazirine via a putative metallo-carbene. The bound structure discloses how an enlarged active site pocket stabilizes the carbene bound to the heme iron and, presumably, the transition state for the formation of this key intermediate. This work demonstrates that improved MicroED technology and the advancement in protein structure prediction now enable investigation of structures that was previously beyond reach.

Biocatalytic Carbene Transfer Using Diazirines.

Biocatalytic carbene transfer from diazo compounds is a versatile strategy in asymmetric synthesis. However, the limited pool of stable diazo compounds constrains the variety of accessible products. To overcome this restriction, we have engineered variants of Aeropyrum pernix protoglobin (ApePgb) that use diazirines as carbene precursors. While the enhanced stability of diazirines relative to their diazo isomers enables access to a diverse array of carbenes, they have previously resisted catalytic activation. Our engineered ApePgb variants represent the first example of catalysts for selective carbene transfer from these species at room temperature. The structure of an ApePgb variant, determined by microcrystal electron diffraction (MicroED), reveals that evolution has enhanced access to the heme active site to facilitate this new-to-nature catalysis. Using readily prepared aryl diazirines as model substrates, we demonstrate the application of these highly stable carbene precursors in biocatalytic cyclopropanation, N-H insertion, and Si-H insertion reactions.

Structural analysis of histone deacetylase 8 mutants associated with Cornelia de Lange Syndrome spectrum disorders.

Cornelia de Lange Syndrome (CdLS) and associated spectrum disorders are characterized by one or more congenital anomalies including distinctive facial features, upper limb abnormalities, intellectual disability, and other symptoms. The molecular genetic basis of CdLS is linked to defects in cohesin, a protein complex that functions in sister chromatid cohesion, chromatin organization, and transcriptional regulation. Histone deacetylase 8 (HDAC8) plays an important role in cohesin function by catalyzing the deacetylation of SMC3, which is required for efficient recycling of the cohesin complex. Missense mutations in HDAC8 have been identified in children diagnosed with CdLS spectrum disorders, and here we outline structure-function relationships for four of these mutations. Specifically, we report the 1.50 Å-resolution structure of the I45T HDAC8-suberoylanilide hydroxamic acid complex, the 1.84 Å-resolution structure of E66D/Y306F HDAC8 complexed with a peptide assay substrate, and the 2.40 Å-resolution structure of G320R HDAC8 complexed with the inhibitor M344. Additionally, we present a computationally generated model of D176G HDAC8. These structures illuminate new structure-function relationships for HDAC8 and highlight the importance of long-range interactions in the protein scaffold that can influence catalytic function.

Unusual zinc-binding mode of HDAC6-selective hydroxamate inhibitors.

Histone deacetylases (HDACs) regulate myriad cellular processes by catalyzing the hydrolysis of acetyl-l-lysine residues in histone and nonhistone proteins. The Zn2+-dependent class IIb enzyme HDAC6 regulates microtubule function by deacetylating α-tubulin, which suppresses microtubule dynamics and leads to cell cycle arrest and apoptosis. Accordingly, HDAC6 is a target for the development of selective inhibitors that might be useful in new therapeutic approaches for the treatment of cancer, neurodegenerative diseases, and other disorders. Here, we present high-resolution structures of catalytic domain 2 from Danio rerio HDAC6 (henceforth simply "HDAC6") complexed with compounds that selectively inhibit HDAC6 while maintaining nanomolar inhibitory potency: N-hydroxy-4-[(N(2-hydroxyethyl)-2-phenylacetamido)methyl)-benzamide)] (HPB), ACY-1215 (Ricolinostat), and ACY-1083. These structures reveal that an unusual monodentate Zn2+ coordination mode is exploited by sterically bulky HDAC6-selective phenylhydroxamate inhibitors. We additionally report the ultrahigh-resolution structure of the HDAC6-trichostatin A complex, which reveals two Zn2+-binding conformers for the inhibitor: a major conformer (70%) with canonical bidentate hydroxamate-Zn2+ coordination geometry and a minor conformer (30%) with monodentate hydroxamate-Zn2+ coordination geometry, reflecting a free energy difference of only 0.5 kcal/mol. The minor conformer is not visible in lower resolution structure determinations. Structural comparisons of HDAC6-inhibitor complexes with class I HDACs suggest active site features that contribute to the isozyme selectivity observed in biochemical assays.

Entropy as a Driver of Selectivity for Inhibitor Binding to Histone Deacetylase 6.

Among the metal-dependent histone deacetylases, the class IIb isozyme HDAC6 is remarkable because of its role in the regulation of microtubule dynamics in the cytosol. Selective inhibition of HDAC6 results in microtubule hyperacetylation, leading to cell cycle arrest and apoptosis, which is a validated strategy for cancer chemotherapy and the treatment of other disorders. HDAC6 inhibitors generally consist of a Zn2+-binding group such as a hydroxamate, a linker, and a capping group; the capping group is a critical determinant of isozyme selectivity. Surprisingly, however, even "capless" inhibitors exhibit appreciable HDAC6 selectivity. To probe the chemical basis for this selectivity, we now report high-resolution crystal structures of HDAC6 complexed with capless cycloalkyl hydroxamate inhibitors 1-4. Each inhibitor hydroxamate group coordinates to the catalytic Zn2+ ion with canonical bidentate geometry. Additionally, the olefin moieties of compounds 2 and 4 bind in an aromatic crevice between the side chains of F583 and F643. Reasoning that similar binding could be achieved in the representative class I isozyme HDAC8, we employed isothermal titration calorimetry to study the thermodynamics of inhibitor binding. These measurements indicate that the entropy of inhibitor binding is generally positive for binding to HDAC6 and negative for binding to HDAC8, resulting in ≤313-fold selectivity for binding to HDAC6 relative to HDAC8. Thus, favorable binding entropy contributes to HDAC6 selectivity. Notably, cyclohexenyl hydroxamate 2 represents a promising lead for derivatization with capping groups that may further enhance its impressive 313-fold thermodynamic selectivity for HDAC6 inhibition.

Structural and Functional Influence of the Glycine-Rich Loop G302GGGY on the Catalytic Tyrosine of Histone Deacetylase 8.

Histone deacetylase 8 (HDAC8) catalyzes the hydrolysis of acetyl-l-lysine to yield products l-lysine and acetate through a mechanism in which a nucleophilic water molecule is activated by a histidine general base and a catalytic metal ion (Zn2+ or Fe2+). Acetyl-l-lysine also requires activation by metal coordination and a hydrogen bond with catalytic tyrosine Y306, which also functions in transition state stabilization. Interestingly, Y306 is located in the conserved glycine-rich loop G302GGGY. The potential flexibility afforded by the tetraglycine segment may facilitate induced-fit conformational changes in Y306 between "in" and "out" positions, as observed in related deacetylases. To probe the catalytic importance of the glycine-rich loop in HDAC8, we rigidified this loop by preparing the G302A, G303A, G304A, and G305A mutants and measured their steady state kinetics and determined their X-ray crystal structures. Substantial losses of catalytic efficiency are observed (10-500-fold based on kcat/KM), particularly for G304A HDAC8 and G305A HDAC8. These mutants also exhibit the greatest structural changes for catalytic tyrosine Y306 (1.3-1.7 Å shifts of the phenolic hydroxyl group). Molecular dynamics simulations further indicate that G304 and G305 undergo pronounced structural changes as residue 306 undergoes a transition between "in" and "out" conformations. Thus, the G304A and G305A substitutions likely compromise the position and conformational changes of Y306 required for substrate activation and transition state stabilization. The G302A and G303A substitutions have less severe catalytic consequences, and these substitutions may influence an internal channel through which product acetate is believed to exit.

Chemodivergent C(sp3)-H and C(sp2)-H Cyanomethylation Using Engineered Carbene Transferases.

The ubiquity of C-H bonds presents an attractive opportunity to elaborate and build complexity in organic molecules. Methods for selective functionalization, however, often must differentiate among multiple chemically similar and, in some cases indistinguishable, C-H bonds. An advantage of enzymes is that they can be finely tuned using directed evolution to achieve control over divergent C-H functionalization pathways. Here, we demonstrate engineered enzymes that effect a new-to-nature C-H alkylation with unparalleled selectivity: two complementary carbene C-H transferases derived from a cytochrome P450 from Bacillus megaterium deliver an α-cyanocarbene into the α-amino C(sp3)-H bonds or the ortho-arene C(sp2)-H bonds of N-substituted arenes. These two transformations proceed via different mechanisms, yet only minimal changes to the protein scaffold (nine mutations, less than 2% of the sequence) were needed to adjust the enzyme's control over the site-selectivity of cyanomethylation. The X-ray crystal structure of the selective C(sp3)-H alkylase, P411-PFA, reveals an unprecedented helical disruption which alters the shape and electrostatics in the enzyme active site. Overall, this work demonstrates the advantages of enzymes as C-H functionalization catalysts for divergent molecular derivatization.

Exploring Structural Determinants of Inhibitor Affinity and Selectivity in Complexes with Histone Deacetylase 6.

Inhibition of histone deacetylase 6 (HDAC6) has emerged as a promising therapeutic strategy for the treatment of cancer, chemotherapy-induced peripheral neuropathy, and neurodegenerative disease. The recent X-ray crystal structure determination of HDAC6 enables an understanding of structural features directing affinity and selectivity in the active site. Here, we present the X-ray crystal structures of five HDAC6-inhibitor complexes that illuminate key molecular features of the inhibitor linker and capping groups that facilitate and differentiate binding to HDAC6. In particular, aromatic and heteroaromatic linkers nestle within an aromatic cleft defined by F583 and F643, and different aromatic linkers direct the capping group toward shallow pockets defined by the L1 loop, the L2 loop, or somewhere in between these pockets. These results expand our understanding of factors contributing to the selective inhibition of HDAC6, particularly regarding interactions that can be targeted in the region of the L2 pocket.

Design and Synthesis of Dihydroxamic Acids as HDAC6/8/10 Inhibitors.

We report the synthesis and evaluation of a class of selective multitarget agents for the inhibition of HDAC6, HDAC8, and HDAC10. The concept for this study grew out of a structural analysis of the two selective inhibitors Tubastatin A (HDAC6/10) and PCI-34051 (HDAC8), which we recognized share the same N-benzylindole core. Hybridization of the two inhibitor structures resulted in dihydroxamic acids with benzyl-indole and -indazole core motifs. These substances exhibit potent activity against HDAC6, HDAC8, and HDAC10, while retaining selectivity over HDAC1, HDAC2, and HDAC3. The best substance inhibited the viability of the SK-N-BE(2)C neuroblastoma cell line with an IC50 value similar to a combination treatment with Tubastatin A and PCI-34051. This compound class establishes a proof of concept for such hybrid molecules and could serve as a starting point for the further development of enhanced HDAC6/8/10 inhibitors.

Preparation of a new construct of human histone deacetylase 8 for the crystallization of enzyme-inhibitor complexes.

The metal-dependent histone deacetylases (HDACs) are critical regulatory enzymes that modulate myriad cellular processes. Implicated in cancer, neurodegenerative diseases, and other clinical disorders, various HDAC isozymes serve as validated drug targets. However, structural similarities among the HDAC isozymes challenge efforts in targeting a single isozyme for therapeutic intervention with an inhibitor. X-ray crystallography remains the premiere technique for studying the chemistry of isozyme-selective inhibition. While crystal structures of many HDAC-inhibitor complexes have been determined, especially with the class I isozyme HDAC8, the study of complexes with large inhibitors is complicated by flexible regions of the protein structure that can hinder crystallization. Here, we outline an approach for the identification of regions in HDAC8 that may hinder crystallization. We also describe protocols for the design and preparation of a truncated HDAC8 construct, HDAC8374, that enabled the successful crystallization and structure determination of the HDAC8-Trapoxin A complex at 1.24Å resolution.

Aza-proline effectively mimics l-proline stereochemistry in triple helical collagen.

The prevalence of l-amino acids in biomolecules has been shown to have teleological importance in biomolecular structure and self-assembly. Recently, biophysical studies have demonstrated that natural l-amino acids can be replaced with non-natural achiral aza-amino acids in folded protein structures such as triple helical collagen. However, the structural consequences of achiral aza-amino acid incorporation has not been elucidated in the context of any relevant folded biomolecule. Herein, we use X-ray crystallography to provide the first atomic resolution crystal structure of an achiral aza-amino acid residue embedded within a folded protein structure, definitively illustrating that achiral aza-proline has the capacity to effectively mimic the stereochemistry of natural amino acids within the context of triple helical collagen. We further corroborate this finding with density functional theory computational analysis showing that the natural l-amino acid stereochemistry for aza-proline is energetically favored when arranged in the aza-proline-hydroxyproline-glycine motif. In addition to providing fundamental insight into peptide and protein structure, the incorporation of achiral stereochemical mimics such as aza-amino acids could have far reaching impacts in areas ranging from synthetic materials to drug design.

Synthesis and Biological Investigation of Phenothiazine-Based Benzhydroxamic Acids as Selective Histone Deacetylase 6 Inhibitors.

The phenothiazine system was identified as a favorable cap group for potent and selective histone deacetylase 6 (HDAC6) inhibitors. Here, we report the preparation and systematic variation of phenothiazines and their analogues containing a benzhydroxamic acid moiety as the zinc-binding group. We evaluated their ability to selectively inhibit HDAC6 by a recombinant HDAC enzyme assay, by determining the protein acetylation levels in cells by western blotting (tubulin vs histone acetylation), and by assessing their effects on various cancer cell lines. Structure-activity relationship studies revealed that incorporation of a nitrogen atom into the phenothiazine framework results in increased potency and selectivity for HDAC6 (more than 500-fold selectivity relative to the inhibition of HDAC1, HDAC4, and HDAC8), as rationalized by molecular modeling and docking studies. The binding mode was confirmed by co-crystallization of the potent azaphenothiazine inhibitor with catalytic domain 2 from Danio rerio HDAC6.

Molecular Basis for the Selective Inhibition of Histone Deacetylase 6 by a Mercaptoacetamide Inhibitor.

Mercaptoacetamide histone deacetylase inhibitors are neuroprotective agents that do not exhibit the genotoxicity associated with more commonly used hydroxamate inhibitors. Here, we present the crystal structure of a selective mercaptoacetamide complexed with the C-terminal catalytic domain of HDAC6. When compared with the structure of a mercaptoacetamide bound to the class I isozyme HDAC8, different interactions are observed with the conserved tandem histidine pair in the active site. These differences likely contribute to the selectivity for inhibition of HDAC6, an important target for cancer chemotherapy and the treatment of neurodegenerative disease.

Histone Deacetylase 6-Selective Inhibitors and the Influence of Capping Groups on Hydroxamate-Zinc Denticity.

Four crystal structures are presented of histone deacetylase 6 (HDAC6) complexes with para-substituted phenylhydromaxamate inhibitors, including bulky peptoids. These structures provide insight regarding the design of capping groups that confer selectivity for binding to HDAC6, specifically with regard to interactions in a pocket formed by the L1 loop. Capping group interactions may also influence hydroxamate-Zn2+ coordination with monodentate or bidentate geometry.

Binding of the Microbial Cyclic Tetrapeptide Trapoxin A to the Class I Histone Deacetylase HDAC8.

Trapoxin A is a microbial cyclic tetrapeptide that is an essentially irreversible inhibitor of class I histone deacetylases (HDACs). The inhibitory warhead is the α,ß-epoxyketone side-chain of (2S,9S)-2-amino-8-oxo-9,10-epoxydecanoic acid (l-Aoe), which mimics the side-chain of the HDAC substrate acetyl-l-lysine. We now report the crystal structure of the HDAC8-trapoxin A complex at 1.24 Å resolution, revealing that the ketone moiety of l-Aoe undergoes nucleophilic attack to form a zinc-bound tetrahedral gem-diolate that mimics the tetrahedral intermediate and its flanking transition states in catalysis. Mass spectrometry, activity measurements, and isothermal titration calorimetry confirm that trapoxin A binds tightly (Kd = 3 ± 1 nM) and does not covalently modify the enzyme, so the epoxide moiety of l-Aoe remains intact. Comparison of the HDAC8-trapoxin A complex with the HDAC6-HC toxin complex provides new insight regarding the inhibitory potency of l-Aoe-containing natural products against class I and class II HDACs.

Histone deacetylase 10 structure and molecular function as a polyamine deacetylase.

Cationic polyamines such as spermidine and spermine are critical in all forms of life, as they regulate the function of biological macromolecules. Intracellular polyamine metabolism is regulated by reversible acetylation and dysregulated polyamine metabolism is associated with neoplastic diseases such as colon cancer, prostate cancer and neuroblastoma. Here we report that histone deacetylase 10 (HDAC10) is a robust polyamine deacetylase, using recombinant enzymes from Homo sapiens (human) and Danio rerio (zebrafish). The 2.85 Å-resolution crystal structure of zebrafish HDAC10 complexed with a transition-state analogue inhibitor reveals that a glutamate gatekeeper and a sterically constricted active site confer specificity for N8-acetylspermidine hydrolysis and disfavour acetyllysine hydrolysis. Both HDAC10 and spermidine are known to promote cellular survival through autophagy. Accordingly, this work sets a foundation for studying the chemical biology of autophagy through the structure-based design of inhibitors that may also serve as new leads for cancer chemotherapy.

Persistent negative effects of alcohol drinking on aspects of novelty-directed behavior in male rhesus macaques.

Humans with histories of prolonged heavy alcohol use exhibit poorer performance on cognitive tasks associated with problem solving, short-term memory, and visuospatial reasoning, even following the cessation of drinking, when compared with healthy controls. It is unclear, however, whether the cognitive problems are a consequence of alcohol exposure or a contributing factor to alcohol-use disorders. Here, we examined the relationship between performance on a novel object recognition (NOR) task and total alcohol consumption (TAC) in adult male rhesus macaques (n = 12; ETH group; trained to self-administer alcohol). NOR performance in this group was assessed prior to induction of alcohol drinking ("pre") and, again, after a 1-year abstinence period ("post") and was compared to the performance of a second group (n = 6; Control group), which was alcohol-naïve. In the NOR task, difficulty was manipulated across three phases by varying specific object features and/or by varying duration of access to objects. For each monkey, we measured aspects of novelty-related behavior including novelty detection, novelty reactivity, and perseverative behavior. TAC during induction and a "free" access period in which the monkey could choose between water and a 4% w/v ethanol solution also was determined. We found that performance deficits in the NOR task were a consequence of high total alcohol intake instead of a predictor of subsequent high intake. Poor NOR performance in drinkers with the highest intakes was characterized by increased perseverative behavior rather than an inability to detect or react to novelty. Finally, the observed deficits are long-lasting - persisting even after a year of abstinence. Given the prevalent and persistent nature of alcohol-induced cognitive deficits in patients in treatment settings, understanding the nature of the deficit and its neural basis could ultimately offer novel treatment approaches based on the reversal of alcohol-induced impairment.