Guided Docking as a Data Generation Approach Facilitates Structure-Based Machine Learning on Kinases.

Drug discovery pipelines nowadays rely on machine learning models to explore and evaluate large chemical spaces. While including 3D structural information is considered beneficial, structural models are hindered by the availability of protein-ligand complex structures. Exemplified for kinase drug discovery, we address this issue by generating kinase-ligand complex data using template docking for the kinase compound subset of available ChEMBL assay data. To evaluate the benefit of the created complex data, we use it to train a structure-based E(3)-invariant graph neural network. Our evaluation shows that binding affinities can be predicted with significantly higher precision by models that take synthetic binding poses into account compared to ligand- or drug-target interaction models alone.

High-resolution crystal structure of RNA kinase ArK1 from G. acetivorans.

U47 phosphorylation (Up47) is a novel tRNA modification discovered recently; it can confer thermal stability and nuclease resistance to tRNAs. U47 phosphorylation is catalyzed by Archaeal RNA kinase (Ark1) in an ATP-dependent manner. However, the structural basis for tRNA and/or ATP binding by Ark1 is unclear. Here, we report the expression, purification, and crystallization studies of Ark1 from G. acetivorans (GaArk1). In addition to the Apo-form structure, one GaArk1-ATP complex was also determined in atomic resolution and revealed the detailed basis for ATP binding by GaArk1. The GaArk1-ATP complex represents the only ATP-bound structure of the Ark1 protein. The majority of the ATP-binding residues are conserved, suggesting that GaArk1 and the homologous proteins share similar mechanism in ATP binding. Sequence and structural analysis further indicated that endogenous guanosine will only inhibit the activities of certain Ark1 proteins, such as Ark1 from T. kodakarensis.

The tunicamycin derivative TunR2 exhibits potent antibiotic properties with low toxicity in an in vivo Mycobacterium marinum-zebrafish TB infection model.

Tunicamycins (TUN) are well-defined, Streptomyces-derived natural products that inhibit protein N-glycosylation in eukaryotes, and by a conserved mechanism also block bacterial cell wall biosynthesis. TUN inhibits the polyprenylphosphate-N-acetyl-hexosamine-1-phospho-transferases (PNPT), an essential family of enzymes found in both bacteria and eukaryotes. We have previously published the development of chemically modified TUN, called TunR1 and TunR2, that have considerably reduced activity on eukaryotes but that retain the potent antibacterial properties. A mechanism for this reduced toxicity has also been reported. TunR1 and TunR2 have been tested against mammalian cell lines in culture and against live insect cells but, until now, no in vivo evaluation has been undertaken for vertebrates. In the current work, TUN, TunR1, and TunR2 are investigated for their relative toxicity and antimycobacterial activity in zebrafish using a well-established Mycobacterium marinum (M. marinum) infection system, a model for studying human Mycobacterium tuberculosis infections. We also report the relative ability to activate the unfolded protein response (UPR), the known mechanism for the eukaryotic toxicity observed with TUN treatment. Importantly, TunR1 and TunR2 retained their antimicrobial properties, as evidenced by a reduction in M. marinum bacterial burden, compared to DMSO-treated zebrafish. In summary, findings from this study highlight the characteristics of recently developed TUN derivatives, mainly TunR2, and its potential for use as a novel anti-bacterial agent for veterinary and potential medical purposes.

N-acetylmuramic acid recognition by MurK kinase from the MurNAc auxotrophic oral pathogen Tannerella forsythia.

The bacterial cell wall consists of a three-dimensional peptidoglycan layer, composed of peptides linked to the sugars N-acetylmuramic acid (MurNAc) and GlcNAc. Unlike other bacteria, the pathogenic Tannerella forsythia, a member of the red complex group of bacteria associated with the late stages of periodontitis, lacks biosynthetic pathways for MurNAc production and therefore obtains MurNAc from the environment. Sugar kinases play a crucial role in the MurNAc recycling process, activating the sugar molecules by phosphorylation. In this study, we present the first crystal structures of a MurNAc kinase, called murein sugar kinase (MurK), in its unbound state as well as in complexes with the ATP analog ß-γ-methylene adenosine triphosphate (AMP-PCP) and with MurNAc. We also determined the crystal structures of K1058, a paralogous MurNAc kinase of T. forsythia, in its unbound state and in complex with MurNAc. We identified the active site and residues crucial for MurNAc specificity as the less bulky side chains of S133, P134, and L135, which enlarge the binding cavity for the lactyl ether group, unlike the glutamate or histidine residues present in structural homologs. In establishing the apparent kinetic parameters for both enzymes, we showed a comparable affinity for MurNAc (Km 180 µM and 30 µM for MurK and K1058, respectively), with MurK being over two hundred times faster than K1058 (Vmax 80 and 0.34 µmol min-1 mg-1, respectively). These data might support a structure-guided approach to development of inhibitory MurNAc analogs for pathogen MurK enzymes.

ß-Boronic Acid-Substituted Bodipy Dyes for Fluorescence Anisotropy Analysis of Carbohydrate Binding.

Boronic acids are widely used for labeling catechols and carbohydrates in analytical (bio)chemistry due to their high binding affinities for diols. Here, we present two asymmetrically substituted Bodipy dyes with a boronic acid at the ß-position (BBB). We present a green-emitting BBB, gBBB, and, by expanding the conjugated system of the Bodipy core at α-position, a red-emitting rBBB. Especially, gBBB shows a bathochromic shift of the electronic spectra upon binding to saccharides and polyols, whereas the fluorescence lifetime of rBBB is more sensitive to hydroxy-ligand binding. Moreover, gBBB constantly shows higher binding affinities than rBBB, reaching Kb ≈ 103 M-1 at pH 8.5 for fructose. Finally, time-resolved fluorescence anisotropy allows to distinguish the number of saccharide units of oligosaccharides as the bond along the transition dipole moment ensures that the fluorescence anisotropy only decays due to the rotational diffusion of labeled carbohydrates. ß-substituted BODIPY dyes are, thus, foreseen as fluorescence anisotropy labels for macromolecule sizing, where conventional fluorophores fail to discriminate due to the chemical similarity of recognition sites.

Metadynamics Simulations to Study the Structural Ensembles and Binding Processes of Intrinsically Disordered Proteins.

The structures of intrinsically disordered proteins (IDPs) are highly dynamic. It is hard to characterize the structures of these proteins experimentally. Molecular dynamics (MD) simulation is a powerful tool in the understanding of protein dynamic structures and function. This chapter describes the application of metadynamics-based enhanced sampling methods in the study of phosphorylation regulation on the structure of kinase-inducible domains (KID). The structural properties of free pKID and KID were obtained by parallel tempering metadynamics combined with well-tempered ensemble (PTMetaD WTE) method, and the binding free energy surfaces of pKID/KID and KIX were characterized by bias-exchanged metadynamics (BE-MetaD) simulations.

Installing a Green Engine To Drive an Enzyme Cascade: A Light-Powered In Vitro Biosystem for Poly(3-hydroxybutyrate) Synthesis.

Many existing in vitro biosystems harness power from the chemical energy contained in substrates and co-substrates, and light or electric energy provided from abiotic parts, leading to a compromise in atom economy, incompatibility between biological and abiotic parts, and most importantly, incapability to spatiotemporally co-regenerate ATP and NADPH. In this study, we developed a light-powered in vitro biosystem for poly(3-hydroxybutyrate) (PHB) synthesis using natural thylakoid membranes (TMs) to regenerate ATP and NADPH for a five-enzyme cascade. Through effective coupling of cofactor regeneration and mass conversion, 20 mM PHB was yielded from 50 mM sodium acetate with a molar conversion efficiency of carbon of 80.0 % and a light-energy conversion efficiency of 3.04 %, which are much higher than the efficiencies of similar in vitro PHB synthesis biosystems. This suggests the promise of installing TMs as a green engine to drive more enzyme cascades.

Protocol for determining the regulation of lipid kinases and changes in phospholipids in vitro.

The regulation of lipid kinases has remained elusive given the difficulties of assessing changes in lipid levels. Here, we describe the isolation of protein and lipid kinases to determine the regulation of lipid kinases in vitro. This can be followed by analysis of effects of regulators on lipid kinase-mediated changes in phospholipids without the use of radioactivity, with a specific focus on PI(5)P generation by the enzyme PIKfyve. For complete details on the use and execution of this protocol, please refer to Karabiyik et al. (2021).

The Escherichia coli two-component signal sensor BarA binds protonated acetate via a conserved hydrophobic-binding pocket.

The BarA/UvrY two-component signal transduction system is widely conserved in γ-proteobacteria and provides a link between the metabolic state of the cells and the Csr posttranscriptional regulatory system. In Escherichia coli, the BarA/UvrY system responds to the presence of acetate and other short-chain carboxylic acids by activating transcription of the noncoding RNAs, CsrB and CsrC, which sequester the RNA-binding protein CsrA, a global regulator of gene expression. However, the state of the carboxyl group in the acetate molecule, which serves as the BarA stimulus, and the signal reception site of BarA remain unknown. In this study, we show that the deletion or replacement of the periplasmic domain of BarA and also the substitution of certain hydroxylated and hydrophobic amino acid residues in this region, result in a sensor kinase that remains unresponsive to its physiological stimulus, demonstrating that the periplasmic region of BarA constitutes a functional detector domain. Moreover, we provide evidence that the protonated state of acetate or formate serves as the physiological stimulus of BarA. In addition, modeling of the BarA sensor domain and prediction of the signal-binding site, by blind molecular docking, revealed a calcium channels and chemotaxis receptors domain with a conserved binding pocket, which comprised uncharged polar and hydrophobic amino acid residues. Based on the comparative sequence and phylogenetic analyses, we propose that, at least, two types of BarA orthologues diverged and evolved separately to acquire distinct signal-binding properties, illustrating the wide adaptability of the bacterial sensor kinase proteins.

Reconstitution of an intact clock reveals mechanisms of circadian timekeeping.

Circadian clocks control gene expression to provide an internal representation of local time. We report reconstitution of a complete cyanobacterial circadian clock in vitro, including the central oscillator, signal transduction pathways, downstream transcription factor, and promoter DNA. The entire system oscillates autonomously and remains phase coherent for many days with a fluorescence-based readout that enables real-time observation of each component simultaneously without user intervention. We identified the molecular basis for loss of cycling in an arrhythmic mutant and explored fundamental mechanisms of timekeeping in the cyanobacterial clock. We find that SasA, a circadian sensor histidine kinase associated with clock output, engages directly with KaiB on the KaiC hexamer to regulate period and amplitude of the central oscillator. SasA uses structural mimicry to cooperatively recruit the rare, fold-switched conformation of KaiB to the KaiC hexamer to form the nighttime repressive complex and enhance rhythmicity of the oscillator, particularly under limiting concentrations of KaiB. Thus, the expanded in vitro clock reveals previously unknown mechanisms by which the circadian system of cyanobacteria maintains the pace and rhythmicity under variable protein concentrations.

A fully reversible 25-hydroxy steroid kinase involved in oxygen-independent cholesterol side-chain oxidation.

The degradation of cholesterol and related steroids by microbes follows fundamentally different strategies in aerobic and anaerobic environments. In anaerobic bacteria, the primary C26 of the isoprenoid side chain is hydroxylated without oxygen via a three-step cascade: (i) water-dependent hydroxylation at the tertiary C25, (ii) ATP-dependent dehydration to form a subterminal alkene, and (iii) water-dependent hydroxylation at the primary C26 to form an allylic alcohol. However, the enzymes involved in the ATP-dependent dehydration have remained unknown. Here, we isolated an ATP-dependent 25-hydroxy-steroid kinase (25-HSK) from the anaerobic bacterium Sterolibacterium denitrificans. This highly active enzyme preferentially phosphorylated the tertiary C25 of steroid alcohols, including metabolites of cholesterol and sitosterol degradation or 25-OH-vitamin D3. Kinetic data were in agreement with a sequential mechanism via a ternary complex. Remarkably, 25-HSK readily catalyzed the formation of γ-(18O)2-ATP from ADP and the C25-(18O)2-phosphoester. The observed full reversibility of 25-HSK with an equilibrium constant below one can be rationalized by an unusual high phosphoryl transfer potential of tertiary steroid C25-phosphoesters, which is ≈20 kJ mol-1 higher than that of standard sugar phosphoesters and even slightly greater than the ß,γ-phosphoanhydride of ATP. In summary, 25-HSK plays an essential role in anaerobic bacterial degradation of zoo- and phytosterols and shows only little similarity to known phosphotransferases.

Modeling a human CLP1 mutation in mouse identifies an accumulation of tyrosine pre-tRNA fragments causing pontocerebellar hypoplasia type 10.

Cleavage factor polyribonucleotide kinase subunit 1 (CLP1), an RNA kinase, plays essential roles in protein complexes involved in the 3'-end formation and polyadenylation of mRNA and the tRNA splicing endonuclease complex, which is involved in precursor tRNA splicing. The mutation R140H in human CLP1 causes pontocerebellar hypoplasia type 10 (PCH10), which is characterized by microcephaly and axonal peripheral neuropathy. Previously, we reported that RNA fragments derived from isoleucine pre-tRNA introns (Ile-introns) accumulate in fibroblasts of patients with PCH10. Therefore, it has been suggested that this intronic RNA fragment accumulation may trigger PCH10 onset. However, the molecular mechanism underlying PCH10 pathogenesis remains elusive. Thus, we generated knock-in mutant mice that harbored a CLP1 mutation consistent with R140H. As expected, these mice showed progressive loss of the upper motor neurons, resulting in impaired locomotor activity, although the phenotype was milder than that of the human variant. Mechanistically, we found that the R140H mutation causes intracellular accumulation of Ile-introns derived from isoleucine pre-tRNAs and 5' tRNA fragments derived from tyrosine pre-tRNAs, suggesting that these two types of RNA fragments were cooperatively or independently involved in the onset and progression of the disease. Taken together, the CLP1-R140H mouse model provided new insights into the pathogenesis of neurodegenerative diseases, such as PCH10, caused by genetic mutations in tRNA metabolism-related molecules.

Identification of kinase inhibitors that rule out the CYP27B1-mediated activation of vitamin D: an integrated machine learning and structure-based drug designing approach.

CYP27B1, a cytochrome P450-containing hydroxylase enzyme, converts vitamin D precursor calcidiol (25-hydroxycholecalciferol) to its active form calcitriol (1α,25(OH)2D3). Tyrosine kinase inhibitor such as imatinib is reported to interfere with the activation of vitamin D3 by inhibiting CYP27B1 enzyme. Consequently, there is a decrease in the serum levels of active vitamin D that in turn may increase the relapse risk among the cancer patients treated with imatinib. Within this framework, the current study focuses on identifying other possible kinase inhibitors that may affect the calcitriol level in the body by inhibiting CYP27B1. To achieve this, we explored multiple machine learning approaches including support vector machine (SVM), random forest (RF), and artificial neural network (ANN) to identify possible CYP27B1 inhibitors from a pool of kinase inhibitors database. The most reliable classification model was obtained from the SVM approach with Matthews correlation coefficient of 0.82 for the external test set. This model was further employed for the virtual screening of kinase inhibitors from the binding database (DB), which tend to interfere with the CYP27B1-mediated activation of vitamin D. This screening yielded around 4646 kinase inhibitors that were further subjected to structure-based analyses using the homology model of CYP27B1, as the 3D structure of CYP27B1 complexed with heme was not available. Overall, five kinase inhibitors including two well-known drugs, i.e., AT7867 (Compound-2) and amitriptyline N-oxide (Compound-3), were found to interact with CYP27B1 in such a way that may preclude the conversion of vitamin D to its active form and hence testify the impairment of vitamin D activation pathway.

Transcriptional regulation of the mannitol phosphotransferase system operon by the ferric uptake regulator (Fur) in Vibrio cholerae El Tor serogroup O1.

The phosphoenolpyruvate (PEP): carbohydrate phosphotransferase system (PTS) allows bacteria to use various carbohydrates as energy resources including mannitol. The mannitol-specific PTS transporter in Vibrio cholerae is encoded by the mtlADR operon. Expression of the mtl operon has been shown to be strictly regulated by CRP, MtlS, and MtlR. In the present study, we investigated the regulation of mtlADR by the ferric uptake regulator (Fur). The results showed that Fur binds to the promoter-proximal DNA region of mtlADR to repress its transcription independent of iron, in mannitol-containing growth medium. The capacity for mannitol fermentation was significantly increased in Δfur relative to that of WT for normal and iron-replete growth media. The level of organic acids produced by Δfur was significantly enhanced relative to that produced by the WT strain in the normal and iron-replete media but not in an iron-starved medium. The results provided for a deeper understanding of the regulation of mtlADR in V. cholerae.

NMR solution structures of Runella slithyformis RNA 2'-phosphotransferase Tpt1 provide insights into NAD+ binding and specificity.

Tpt1, an essential component of the fungal and plant tRNA splicing machinery, catalyzes transfer of an internal RNA 2'-PO4 to NAD+ yielding RNA 2'-OH and ADP-ribose-1',2'-cyclic phosphate products. Here, we report NMR structures of the Tpt1 ortholog from the bacterium Runella slithyformis (RslTpt1), as apoenzyme and bound to NAD+. RslTpt1 consists of N- and C-terminal lobes with substantial inter-lobe dynamics in the free and NAD+-bound states. ITC measurements of RslTpt1 binding to NAD+ (KD ∼31 µM), ADP-ribose (∼96 µM) and ADP (∼123 µM) indicate that substrate affinity is determined primarily by the ADP moiety; no binding of NMN or nicotinamide is observed by ITC. NAD+-induced chemical shift perturbations (CSPs) localize exclusively to the RslTpt1 C-lobe. NADP+, which contains an adenylate 2'-PO4 (mimicking the substrate RNA 2'-PO4), binds with lower affinity (KD ∼1 mM) and elicits only N-lobe CSPs. The RslTpt1·NAD+ binary complex reveals C-lobe contacts to adenosine ribose hydroxyls (His99, Thr101), the adenine nucleobase (Asn105, Asp112, Gly113, Met117) and the nicotinamide riboside (Ser125, Gln126, Asn163, Val165), several of which are essential for RslTpt1 activity in vivo. Proximity of the NAD+ ß-phosphate to ribose-C1″ suggests that it may stabilize an oxocarbenium transition-state during the first step of the Tpt1-catalyzed reaction.

1,6-Hexanediol, commonly used to dissolve liquid-liquid phase separated condensates, directly impairs kinase and phosphatase activities.

The concept of liquid-liquid phase separation (LLPS) has emerged as an intriguing mechanism for the organization of membraneless compartments in cells. The alcohol 1,6-hexanediol is widely used as a control to dissolve LLPS assemblies in phase separation studies in diverse fields. However, little is known about potential side effects of 1,6-hexanediol, which could compromise data interpretation and mislead the scientific debate. To examine this issue, we analyzed the effect of 1,6-hexanediol on the activities of various enzymes in vitro. Already at 1% volume concentration, 1,6-hexanediol strongly impaired kinases and phosphatases and partly blocked DNA polymerases, while it had no effect on DNase activity. At concentrations that are usually used to dissolve LLPS droplets (5-10%), both kinases and phosphatases were virtually inactive. Given the widespread function of protein phosphorylation in cells, our data argue for a careful review of 1,6-hexanediol in phase separation studies.

A BLUS1 kinase signal and a decrease in intercellular CO2 concentration are necessary for stomatal opening in response to blue light.

Light-induced stomatal opening stimulates CO2 uptake and transpiration in plants. Weak blue light under strong red light effectively induces stomatal opening. Blue light-dependent stomatal opening initiates light perception by phototropins, and the signal is transmitted to a plasma membrane H+-ATPase in guard cells via BLUE LIGHT SIGNALING 1 (BLUS1) kinase. However, it is unclear how BLUS1 transmits the signal to H+-ATPase. Here, we characterized BLUS1 signaling in Arabidopsis thaliana, and showed that the BLUS1 C-terminus acts as an auto-inhibitory domain and that phototropin-mediated Ser-348 phosphorylation within the domain removes auto-inhibition. C-Terminal truncation and phospho-mimic Ser-348 mutation caused H+-ATPase activation in the dark, but did not elicit stomatal opening. Unexpectedly, the plants exhibited stomatal opening under strong red light and stomatal closure under weak blue light. A decrease in intercellular CO2 concentration via red light-driven photosynthesis together with H+-ATPase activation caused stomatal opening. Furthermore, phototropins caused H+-ATPase dephosphorylation in guard cells expressing constitutive signaling variants of BLUS1 in response to blue light, possibly for fine-tuning stomatal opening. Overall, our findings provide mechanistic insights into the blue light regulation of stomatal opening.

Integrating Multimeric Threading With High-throughput Experiments for Structural Interactome of Escherichia coli.

Genome-wide protein-protein interaction (PPI) determination remains a significant unsolved problem in structural biology. The difficulty is twofold since high-throughput experiments (HTEs) have often a relatively high false-positive rate in assigning PPIs, and PPI quaternary structures are more difficult to solve than tertiary structures using traditional structural biology techniques. We proposed a uniform pipeline, Threpp, to address both problems. Starting from a pair of monomer sequences, Threpp first threads both sequences through a complex structure library, where the alignment score is combined with HTE data using a naïve Bayesian classifier model to predict the likelihood of two chains to interact with each other. Next, quaternary complex structures of the identified PPIs are constructed by reassembling monomeric alignments with dimeric threading frameworks through interface-specific structural alignments. The pipeline was applied to the Escherichia coli genome and created 35,125 confident PPIs which is 4.5-fold higher than HTE alone. Graphic analyses of the PPI networks show a scale-free cluster size distribution, consistent with previous studies, which was found critical to the robustness of genome evolution and the centrality of functionally important proteins that are essential to E. coli survival. Furthermore, complex structure models were constructed for all predicted E. coli PPIs based on the quaternary threading alignments, where 6771 of them were found to have a high confidence score that corresponds to the correct fold of the complexes with a TM-score >0.5, and 39 showed a close consistency with the later released experimental structures with an average TM-score = 0.73. These results demonstrated the significant usefulness of threading-based homologous modeling in both genome-wide PPI network detection and complex structural construction.

Advancements in chemical biology targeting the kinases and phosphatases of RNA polymerase II-mediated transcription.

Phosphorylation of RNA polymerase II (RNAP II) coordinates the temporal progression of eukaryotic transcription. The development and application of chemical genetic methods have enhanced our ability to investigate the intricate and intertwined pathways regulated by the kinases and phosphatases targeting RNAP II to ensure transcription accuracy and efficiency. Although identifying small molecules that modulate these enzymes has been challenging due to their highly conserved structures, powerful new chemical biology strategies such as targeted covalent inhibitors and small molecule degraders have significantly improved chemical probe specificity. The recent success in discovering phosphatase holoenzyme activators and inhibitors, which demonstrates the feasibility of selective targeting of individual phosphatase complexes, opens up new avenues into the study of transcription. Herein, we summarize how chemical biology is used to delineate kinases' identities involved in RNAP II regulation and new concepts in inhibitor/activator design implemented for kinases/phosphatases involved in modulating RNAP II-mediated transcription.

Crystal structure and molecular dynamics simulations of a promiscuous ancestor reveal residues and an epistatic interaction involved in substrate binding and catalysis in the ATP-dependent vitamin kinase family members.

Enzymes with hydroxymethylpyrimidine/phosphomethylpyrimidine kinase activity (HMPPK) are essential in the vitamin B1 (thiamine pyrophosphate) biosynthesis and recycling pathways. In contrast, enzymes with pyridoxal kinase activity (PLK) produce pyridoxal phosphate (vitamin B6), an essential cofactor for various biochemical reactions. In the ATP-dependent vitamin kinases family, the members of PLK/HMPPK-like subfamily have both enzymatic activities. It has been proposed that the promiscuous PLK activity of ancestral HMPPK enzymes could have been the starting point for this activity. In earlier work, we reconstructed the ancestral sequences of this family and characterized the substrate specificity of the common ancestor between PLK/HMPPK-like and HMPPK enzymes (AncC). From these studies, the Gln45Met mutation was proposed as a critical event for the PLK activity emergence. Here, we crystallize and determine the AncC structure by X-ray crystallography and assess the role of the Gln45Met mutation by site-directed mutagenesis. Kinetic characterization of this mutant shows a significant increase in the PL affinity. Through molecular dynamics simulation and MM/PBSA calculations some residues, important for substrate interactions and catalysis, were identified in the wild type and in the mutated ancestor. Interestingly, a strong epistatic interaction responsible for the evolutionary pathway of the PLK activity in PLK/HMPPK-like enzymes was revealed. Also, other putative mutations relevant to PLK activity in modern PLK/HMPPK-like enzymes were identified.