Computational structural genomics and clinical evidence suggest BCKDK gain-of-function may cause a potentially asymptomatic maple syrup urine disease phenotype.

Maple syrup urine disease (MSUD) is a disorder of branched-chain amino acid metabolism caused by a defect in the branched-chain α-ketoacid dehydrogenase (BCKD) complex (OMIM #248600). The hallmark presentation is encephalopathic crisis in neonates, but can also present with metabolic decompensation, developmental delays, and feeding difficulties. Biochemical evidence for MSUD includes elevated branched-chain amino acids (BCAA) and the pathognomonic presence of alloisoleucine. The BCKD complex contains several subunits associated with autosomal recessive MSUD, while its regulatory proteins have less well-defined disease associations. We report on two families with the same BCKDK variant (c.1115C>G (p.Thr372Arg)). Probands were detected on newborn screening and demonstrated biochemical evidence of MSUD. The variant was identified in reportedly asymptomatic parents and additional family members who had elevated BCAA and alloisoleucine, following an autosomal dominant pattern of inheritance. To better define the functional effect of the variant on the kinase, we completed molecular modeling using sequence-based (2D), structural-based (3D), and dynamic-based (4D) analyses. The BCKDK variant modeling indicated a gain-of-function which leads to impaired BCAA catabolism consistent with the biochemical evidence in this cohort. Combining the evidence gained from molecular modeling with the absence of metabolic decompensation in our patients and several adult family members, despite encountering stressors typically problematic in classic MSUD, we suggest that heterozygous gain-of-function variants in BCKDK may represent a novel biochemical phenotype of MSUD with a benign clinical course.

Deep computational phenotyping of genomic variants impacting the SET domain of KMT2C reveal molecular mechanisms for their dysfunction.

Introduction: Kleefstra Syndrome type 2 (KLEFS-2) is a genetic, neurodevelopmental disorder characterized by intellectual disability, infantile hypotonia, severe expressive language delay, and characteristic facial appearance, with a spectrum of other distinct clinical manifestations. Pathogenic mutations in the epigenetic modifier type 2 lysine methyltransferase KMT2C have been identified to be causative in KLEFS-2 individuals. Methods: This work reports a translational genomic study that applies a multidimensional computational approach for deep variant phenotyping, combining conventional genomic analyses, advanced protein bioinformatics, computational biophysics, biochemistry, and biostatistics-based modeling. We use standard variant annotation, paralog annotation analyses, molecular mechanics, and molecular dynamics simulations to evaluate damaging scores and provide potential mechanisms underlying KMT2C variant dysfunction. Results: We integrated data derived from the structure and dynamics of KMT2C to classify variants into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). When compared with controls, these variants show values reflecting alterations in molecular fitness in both structure and dynamics. Discussion: We demonstrate that our 3D models for KMT2C variants suggest distinct mechanisms that lead to their imbalance and are not predictable from sequence alone. Thus, the missense variants studied here cause destabilizing effects on KMT2C function by different biophysical and biochemical mechanisms which we adeptly describe. This new knowledge extends our understanding of how variations in the KMT2C gene cause the dysfunction of its methyltransferase enzyme product, thereby bearing significant biomedical relevance for carriers of KLEFS2-associated genomic mutations.

A multi-layered computational structural genomics approach enhances domain-specific interpretation of Kleefstra syndrome variants in EHMT1.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

A Multi-Layered Computational Structural Genomics Approach Enhances Domain-Specific Interpretation of Kleefstra Syndrome Variants in EHMT1.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretation of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within this domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism of variant dysfunction.

Writers and readers of H3K9me2 form distinct protein networks during the cell cycle that include candidates for H3K9 mimicry.

Histone H3 lysine 9 methylation (H3K9me), which is written by the Euchromatic Histone Lysine Methyltransferases EHMT1 and EHMT2 and read by the heterochromatin protein 1 (HP1) chromobox (CBX) protein family, is dysregulated in many types of cancers. Approaches to inhibit regulators of this pathway are currently being evaluated for therapeutic purposes. Thus, knowledge of the complexes supporting the function of these writers and readers during the process of cell proliferation is critical for our understanding of their role in carcinogenesis. Here, we immunopurified each of these proteins and used mass spectrometry to define their associated non-histone proteins, individually and at two different phases of the cell cycle, namely G1/S and G2/M. Our findings identify novel binding proteins for these writers and readers, as well as corroborate known interactors, to show the formation of distinct protein complex networks in a cell cycle phase-specific manner. Furthermore, there is an organizational switch between cell cycle phases for interactions among specific writer-reader pairs. Through a multi-tiered bioinformatics-based approach, we reveal that many interacting proteins exhibit histone mimicry, based on an H3K9-like linear motif. Gene ontology analyses, pathway enrichment, and network reconstruction inferred that these comprehensive EHMT and CBX-associated interacting protein networks participate in various functions, including transcription, DNA repair, splicing, and membrane disassembly. Combined, our data reveals novel complexes that provide insight into key functions of cell cycle-associated epigenomic processes that are highly relevant for better understanding these chromatin-modifying proteins during cell cycle and carcinogenesis.

Structural bioinformatics enhances the interpretation of somatic mutations in KDM6A found in human cancers.

The histone demethylase KDM6A has recently elicited significant attention because its mutations are associated with a rare congenital disorder (Kabuki syndrome) and various types of human cancers. However, distinguishing KDM6A mutations that are deleterious to the enzyme and their underlying mechanisms of dysfunction remain to be fully understood. Here, we report the results from a multi-tiered approach evaluating the impact of 197 KDM6A somatic mutations using information derived from combining conventional genomics data with computational biophysics. This comprehensive approach incorporates multiple scores derived from alterations in protein sequence, structure, and molecular dynamics. Using this method, we classify the KDM6A mutations into 136 damaging variants (69.0%), 32 tolerated variants (16.2%), and 29 variants of uncertain significance (VUS, 14.7%), which is a significant improvement from the previous classification based on the conventional tools (over 40% VUS). We further classify the damaging variants into 15 structural variants (SV), 88 dynamic variants (DV), and 33 structural and dynamic variants (SDV). Comparison with variant scoring methods used in current clinical diagnosis guidelines demonstrates that our approach provides a more comprehensive evaluation of damaging potential and reveals mechanisms of dysfunction. Thus, these results should be taken into consideration for clinical assessment of the damaging potential of each mutation, as they provide hypotheses for experimental validation and critical information for the development of mutant-specific drugs to fight diseases caused by KDM6A dysfunctions.

Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis.

Disruptor of telomeric silencing 1-like (DOT1L) is the only non-SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time-dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1-AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L-AF9 complex are less extensive and highly dynamic, resembling a swivel-like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

Molecular mechanics and dynamic simulations of well-known Kabuki syndrome-associated KDM6A variants reveal putative mechanisms of dysfunction.

BACKGROUND: Kabuki syndrome is a genetic disorder that affects several body systems and presents with variations in symptoms and severity. The syndrome is named for a common phenotype of faces resembling stage makeup used in a Japanese traditional theatrical art named kabuki. The most frequent cause of this syndrome is mutations in the H3K4 family of histone methyltransferases while a smaller percentage results from genetic alterations affecting the histone demethylase, KDM6A. Because of the rare presentation of the latter form of the disease, little is known about how missense changes in the KDM6A protein sequence impact protein function. RESULTS: In this study, we use molecular mechanic and molecular dynamic simulations to enhance the annotation and mechanistic interpretation of the potential impact of eleven KDM6A missense variants found in Kabuki syndrome patients. These variants (N910S, D980V, S1025G, C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W, and R1351Q) are predicted to be pathogenic, likely pathogenic or of uncertain significance by sequence-based analysis. Here, we demonstrate, for the first time, that although Kabuki syndrome missense variants are found outside the functionally critical regions, they could affect overall function by significantly disrupting global and local conformation (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), chemical environment (C1153R, C1153Y, P1195L, L1200F, Q1212R, Q1248R, R1255W and R1351Q), and/or molecular dynamics of the catalytic domain (all variants). In addition, our approaches predict that many mutations, in particular C1153R, could allosterically disrupt the key enzymatic interactions of KDM6A. CONCLUSIONS: Our study demonstrates that the KDM6A Kabuki syndrome variants may impair histone demethylase function through various mechanisms that include altered protein integrity, local environment, molecular interactions and protein dynamics. Molecular dynamics simulations of the wild type and the variants are critical to gain a better understanding of molecular dysfunction. This type of comprehensive structure- and MD-based analyses should help develop improved impact scoring systems to interpret the damaging effects of variants in this protein and other related proteins as well as provide detailed mechanistic insight that is not currently predictable from sequence alone.

Protein lysine 43 methylation by EZH1 promotes AML1-ETO transcriptional repression in leukemia.

The oncogenic fusion protein AML1-ETO retains the ability of AML1 to interact with the enhancer core DNA sequences, but blocks AML1-dependent transcription. Previous studies have shown that post-translational modification of AML1-ETO may play a role in its regulation. Here we report that AML1-ETO-positive patients, with high histone lysine methyltransferase Enhancer of zeste homolog 1 (EZH1) expression, show a worse overall survival than those with lower EZH1 expression. EZH1 knockdown impairs survival and proliferation of AML1-ETO-expressing cells in vitro and in vivo. We find that EZH1 WD domain binds to the AML1-ETO NHR1 domain and methylates AML1-ETO at lysine 43 (Lys43). This requires the EZH1 SET domain, which augments AML1-ETO-dependent repression of tumor suppressor genes. Loss of Lys43 methylation by point mutation or domain deletion impairs AML1-ETO-repressive activity. These findings highlight the role of EZH1 in non-histone lysine methylation, indicating that cooperation between AML1-ETO and EZH1 and AML1-ETO site-specific lysine methylation promote AML1-ETO transcriptional repression in leukemia.

ErbB3-binding protein 1 (EBP1) represses HNF4α-mediated transcription and insulin secretion in pancreatic ß-cells.

HNF4α (hepatocyte nuclear factor 4α) is one of the master regulators of pancreatic ß-cell development and function, and mutations in the HNF4α gene are well-known monogenic causes of diabetes. As a member of the nuclear receptor family, HNF4α exerts its gene regulatory function through various molecular interactions; however, there is a paucity of knowledge of the different functional complexes in which HNF4α participates. Here, to find HNF4α-binding proteins in pancreatic ß-cells, we used yeast two-hybrid screening, a mammalian two-hybrid assay, and glutathione S-transferase pulldown approaches, which identified EBP1 (ErbB3-binding protein 1) as a factor that binds HNF4α in a LXXLL motif-mediated manner. In the ß-cells, EBP1 suppressed the expression of HNF4α target genes that are implicated in insulin secretion, which is impaired in HNF4α mutation-driven diabetes. The crystal structure of the HNF4α ligand-binding domain in complex with a peptide harboring the EBP1 LXXLL motif at 3.15Å resolution hinted at the molecular basis of the repression. The details of the structure suggested that EBP1's LXXLL motif competes with HNF4α coactivators for the same binding pocket and thereby prevents recruitment of additional transcriptional coactivators. These findings provide further evidence that EBP1 plays multiple cellular roles and is involved in nuclear receptor-mediated gene regulation. Selective disruption of the HNF4α-EBP1 interaction or tissue-specific EBP1 inactivation can enhance HNF4α activities and thereby improve insulin secretion in ß-cells, potentially representing a new strategy for managing diabetes and related metabolic disorders.

Dimerization defective MODY mutations of hepatocyte nuclear factor 4α.

HNF4α is a culprit gene product for a monogenic and dominantly-inherited form of diabetes, referred to as MODY1 (Maturity Onset Diabetes of the Young type 1). Reduced HNF4α activities have been linked to impaired insulin secretion and ß-cell function. Numerous mutations have been identified from the patients and they have been instructive as to the individual residue's role in protein structure-function and dysfunction. As a member of the nuclear receptor (NR) superfamily, HNF4α is made of characteristic modular domains and it functions exclusively as a homodimer despite its sequence homology to RXR, a common heterodimer partner of non-steroidal NRs. Transcription factors commonly dimerize to enhance their molecular functions mainly by facilitating the recognition of double helix target DNAs that display an intrinsic pseudo-2-fold symmetry and the recruitment of the remainder of the main transcriptional machinery. HNF4α is no exception and its dimerization is maintained by the ligand binding domain (LBD) mainly through the leucine-zipper-like interactions at the stalk of two interacting helices. Although many MODY1 mutations have been previously characterized, including DNA binding disruptors, ligand binding disruptors, coactivator binding disruptors, and protein stability disruptors, protein dimerization disruptors have not been formally reported. In this report, we present a set of data for the two MODY1 mutations found right at the dimerization interface (L332 P and L328del mutations) which clearly exhibit the disruptive effects of directly affecting dimerization, protein stability, and transcriptional activities. These data reinforced the fact that MODY mutations are loss-of-function mutations and HNF4α dimerization is essential for its optimal function and normal physiology.

Discovery of Novel Pyruvate Dehydrogenase Kinase 4 Inhibitors for Potential Oral Treatment of Metabolic Diseases.

Pyruvate dehydrogenase kinase 4 (PDK4) activation is associated with metabolic diseases including hyperglycemia, insulin resistance, allergies, and cancer. Structural modifications of hit anthraquinone led to the identification of a new series of allosteric PDK4 inhibitors. Among this series, compound 8c showed promising in vitro activity with an IC50 value of 84 nM. Good metabolic stability, pharmacokinetic profiles, and possible metabolites were suggested. Compound 8c improved glucose tolerance in diet-induced obese mice and ameliorated allergic reactions in a passive cutaneous anaphylaxis mouse model. Additionally, compound 8c exhibited anticancer activity by controlling cell proliferation, transformation, and apoptosis. From the molecular docking studies, compound 8c displayed optimal fitting in the lipoamide binding site (allosteric) with a full fitness, providing a new scaffold for drug development toward PDK4 inhibitors.

Crystal structures reveal a new and novel FoxO1 binding site within the human glucose-6-phosphatase catalytic subunit 1 gene promoter.

Human glucose-6-phosphatase plays a vital role in blood glucose homeostasis and holds promise as a therapeutic target for diabetes. Expression of its catalytic subunit gene 1 (G6PC1) is tightly regulated by metabolic-response transcription factors such as FoxO1 and CREB. Although at least three potential FoxO1 binding sites (insulin response elements, IREs) and one CREB binding site (cAMP response element, CRE) within the proximal region of the G6PC1 promoter have been identified, the interplay between FoxO1 and CREB and between FoxO1 bound at multiple IREs has not been well characterized. Here we present the crystal structures of the FoxO1 DNA binding domain in complex with the G6PC1 promoter. These complexes reveal the presence of a new non-consensus FoxO1 binding site that overlaps the CRE, suggesting a mutual exclusion mechanism for FoxO1 and CREB binding at the G6PC1 promoter. Additional findings include (i) non-canonical FoxO1 recognition sites, (ii) incomplete FoxO1 occupancies at the available IRE sites, and (iii) FoxO1 dimeric interactions that may play a role in stabilizing DNA looping. These findings provide insight into the regulation of G6PC1 gene transcription by FoxO1, and demonstrate a high versatility of target gene recognition by FoxO1 that correlates with its diverse roles in biology.

Dub3 inhibition suppresses breast cancer invasion and metastasis by promoting Snail1 degradation.

Snail1, a key transcription factor of epithelial-mesenchymal transition (EMT), is subjected to ubiquitination and degradation, but the mechanism by which Snail1 is stabilized in tumours remains unclear. We identify Dub3 as a bona fide Snail1 deubiquitinase, which interacts with and stabilizes Snail1. Dub3 is overexpressed in breast cancer; knockdown of Dub3 resulted in Snail1 destabilization, suppressed EMT and decreased tumour cell migration, invasion, and metastasis. These effects are rescued by ectopic Snail1 expression. IL-6 also stabilizes Snail1 by inducing Dub3 expression, the specific inhibitor WP1130 binds to Dub3 and inhibits the Dub3-mediating Snail1 stabilization in vitro and in vivo. Our study reveals a critical Dub3-Snail1 signalling axis in EMT and metastasis, and provides an effective therapeutic approach against breast cancer.

Repression of HNF1α-mediated transcription by amino-terminal enhancer of split (AES).

HNF1α (Hepatocyte Nuclear Factor 1α) is one of the master regulators in pancreatic beta-cell development and function, and the mutations in Hnf1α are the most common monogenic causes of diabetes mellitus. As a member of the POU transcription factor family, HNF1α exerts its gene regulatory function through various molecular interactions; however, there is a paucity of knowledge in their functional complex formation. In this study, we identified the Groucho protein AES (Amino-terminal Enhancer of Split) as a HNF1α-specific physical binding partner and functional repressor of HNF1α-mediated transcription, which has a direct link to glucose-stimulated insulin secretion in beta-cells that is impaired in the HNF1α mutation-driven diabetes.

Targeting epidermal fatty acid binding protein for treatment of experimental autoimmune encephalomyelitis.

BACKGROUND: Multiple sclerosis (MS) is an autoimmune disease in which dysregulated immune cells attack myelin in the central nervous system (CNS), leading to irreversible neuronal degeneration. Our previous studies have demonstrated that epidermal fatty acid binding protein (E-FABP), widely expressed in immune cells, in particular in dendritic cells (DCs) and T lymphocytes, fuels the overactive immune responses in the mouse model of experimental autoimmune encephalomyelitis (EAE). METHODS: In the present study, we conducted an intensive computational docking analysis to identify novel E-FABP inhibitors for regulation of immune cell functions and for treatment of EAE. RESULTS: We demonstrate that compound [2-(4-acetylphenoxy)-9,10-dimethoxy-6,7-dihydropyrimido[6,1-a]isoquinolin-4-one; designated as EI-03] bound to the lipid binding pocket of E-FABP and enhanced the expression of peroxisome proliferator-activating receptor (PPAR) γ. Further in vitro experiments showed that EI-03 regulated DC functions by inhibition of TNFα production while promoting IL-10 secretion. Moreover, EI-03 treatment counterregulated T cell balance by decreasing effector T cell differentiation (e.g. Th17, Th1) while increasing regulatory T cell development. Most importantly, mice treated with this newly identified compound exhibited reduced clinical symptoms of EAE in mouse models. CONCLUSIONS: Taken together, we have identified a new compound which displays a potential therapeutic benefit for treatment of MS by targeting E-FABP.

Inhibition of tumor growth by a newly-identified activator for epidermal fatty acid binding protein.

Our previous studies have demonstrated that expression of epidermal fatty acid binding protein (E-FABP) in tumor associated macrophages (TAMs) promotes macrophage anti-tumor activity by enhancing IFNß responses in tumor models. Thus, E-FABP represents a new protective factor in enhancing tumor immune surveillance against tumor development. Herein, we report the compound 5-(benzylamino)-2-(3-methylphenyl)-1,3-oxazole-4-carbonitrile (designated EI-05) as a novel E-FABP activator for inhibition of mammary tumor growth. EI-05 was selected from the ZINC compound library using molecular docking analysis based on the crystal structure of E-FABP. Although EI-05 is unable to bind E-FABP directly, it significantly increases E-FABP expression in macrophages during inflammation. Stimulation of macrophages with EI-05 remarkably enhances lipid droplet formation and IFNß production, which further promotes the anti-tumor activity of macrophages. Importantly, administering EI-05 in vivo significantly inhibits mammary tumor growth in a syngeneic mouse model. Altogether, these results suggest that EI-05 may represent a promising drug candidate for anti-tumor treatment through enhancing E-FABP activity and IFNß responses in macrophages.