PKLR mutations in pyruvate kinase deficient Polish patients: Functional characteristics of c.101-1G > A and c.1058delAAG variants.

Pyruvate kinase (PK) deficiency is a rare autosomal recessive disorder characterized by chronic hemolytic anemia of variable severity. Nine Polish patients with severe hemolytic anemia but normal PK activity were found to carry mutations in the PKLR gene encoding PK, five already known ones and one novel (c.178C > T). We characterized two of the known variants by molecular modeling (c.1058delAAG) and minigene splicing analysis (c.101-1G > A). The former gives a partially destabilized PK tetramer, likely of suboptimal activity, and the c.101-1G > A variant gives alternatively spliced mRNA carrying a premature stop codon, encoding a severely truncated PK and likely undergoing nonsense-mediated decay.

Effect of histidine protonation state on ligand binding at the ATP-binding site of human protein kinase CK2.

Histidine residues contribute to numerous molecular interactions, owing to their structure with the ionizable aromatic side chain with pKa close to the physiological pH. Herein, we studied how the two histidine residues, His115 and His160 of the catalytic subunit of human protein kinase CK2, affect the binding of the halogenated heterocyclic ligands at the ATP-binding site. Thermodynamic studies on the interaction between five variants of hCK2α (WT protein and four histidine mutants) and three ionizable bromo-benzotriazoles and their conditionally non-ionizable benzimidazole counterparts were performed with nanoDSF, MST, and ITC. The results allowed us to identify the contribution of interactions involving the particular histidine residues to ligand binding. We showed that despite the well-documented hydrogen bonding/salt bridge formation dragging the anionic ligands towards Lys68, the protonated His160 also contributes to the binding of such ligands by long-range electrostatic interactions. Simultaneously, His 115 indirectly affects ligand binding, placing the hinge region in open/closed conformations.

The relationship between EMG high frequency and low frequency band amplitude changes correlates with tissue inorganic phosphate levels.

Assessing inorganic phosphate levels seems crucial in deciphering the biochemical state of organisms or tissues. The concentration of inorganic phosphate in blood is an order of magnitude smaller than in tissues and, on top of that, it is dynamically used to fill temporary gaps in tissues. This is the reason blood inorganic phosphate level is considered a poor proxy for tissue levels. Therefore, tissue biopsy seems to be the dominant method when assessing inorganic phosphate levels for instance in muscles. In this study, we attempted to derive a non-invasive biomarker for phosphate tissue levels. We analyzed surface electromyography signals taken during 31P spectroscopy of leg muscles in five adult pigs. We induced hypophosphatemia via 20 minutes-long hyperventilation. It turned out that the proportion of the amplitude of the low frequency band and the high frequency band is significantly (p=0.002) correlated with the relative phosphate levels. The electromyographic signal did not correlate significantly with pCO2 levels in the blood, suggesting that the changes in the signal are a result of inorganic phosphate levels, not hyperventilation. The results might lead to the development of a real-time phosphate fluctuations measurement procedure.

Analysis of MT-ATP8 gene variants reported in patients by modeling in silico and in yeast model organism.

Defects in ATP synthase functioning due to the substitutions in its two mitochondrially encoded subunits a and 8 lead to untreatable mitochondrial diseases. Defining the character of variants in genes encoding these subunits is challenging due to their low frequency, heteroplasmy of mitochondrial DNA in patients' cells and polymorphisms of mitochondrial genome. We successfully used yeast S. cerevisiae as a model to study the effects of variants in MT-ATP6 gene and our research led to understand how eight amino acid residues substitutions impact the proton translocation through the channel formed by subunit a and c-ring of ATP synthase at the molecular level. Here we applied this approach to study the effects of the m.8403T>C variant in MT-ATP8 gene. The biochemical data from yeast mitochondria indicate that equivalent mutation is not detrimental for the yeast enzyme functioning. The structural analysis of substitutions in subunit 8 introduced by m.8403T>C and five other variants in MT-ATP8 provides indications about the role of subunit 8 in the membrane domain of ATP synthase and potential structural consequences of substitutions in this subunit.

Divergent contribution of the MVA and MEP pathways to the formation of polyprenols and dolichols in Arabidopsis.

Isoprenoids, including dolichols (Dols) and polyprenols (Prens), are ubiquitous components of eukaryotic cells. In plant cells, there are two pathways that produce precursors utilized for isoprenoid biosynthesis: the mevalonate (MVA) pathway and the methylerythritol phosphate (MEP) pathway. In this work, the contribution of these two pathways to the biosynthesis of Prens and Dols was addressed using an in planta experimental model. Treatment of plants with pathway-specific inhibitors and analysis of the effects of various light conditions indicated distinct biosynthetic origin of Prens and Dols. Feeding with deuteriated, pathway-specific precursors revealed that Dols, present in leaves and roots, were derived from both MEP and MVA pathways and their relative contributions were modulated in response to precursor availability. In contrast, Prens, present in leaves, were almost exclusively synthesized via the MEP pathway. Furthermore, results obtained using a newly introduced here 'competitive' labeling method, designed so as to neutralize the imbalance of metabolic flow resulting from feeding with a single pathway-specific precursor, suggest that under these experimental conditions one fraction of Prens and Dols is synthesized solely from endogenous precursors (deoxyxylulose or mevalonate), while the other fraction is synthesized concomitantly from endogenous and exogenous precursors. Additionally, this report describes a novel methodology for quantitative separation of 2H and 13C distributions observed for isotopologues of metabolically labeled isoprenoids. Collectively, these in planta results show that Dol biosynthesis, which uses both pathways, is significantly modulated depending on pathway productivity, while Prens are consistently derived from the MEP pathway.

Probing the pathogenicity of patient-derived variants of MT-ATP6 in yeast.

The list of mitochondrial DNA (mtDNA) variants detected in individuals with neurodegenerative diseases is constantly growing. Evaluating their functional consequences and pathogenicity is not easy, especially when they are found in only a limited number of patients together with wild-type mtDNA (heteroplasmy). Owing to its amenability to mitochondrial genetic transformation and incapacity to stably maintain heteroplasmy, and the strong evolutionary conservation of the proteins encoded in mitochondria, Saccharomyces cerevisiae provides a convenient model to investigate the functional consequences of human mtDNA variants. We herein report the construction and energy-transducing properties of yeast models of eight MT-ATP6 gene variants identified in patients with various disorders: m.8843T>C, m.8950G>A, m.9016A>G, m.9025G>A, m.9029A>G, m.9058A>G, m.9139G>A and m.9160T>C. Significant defect in growth dependent on respiration and deficits in ATP production were observed in yeast models of m.8950G>A, m.9025G>A and m.9029A>G, providing evidence of pathogenicity for these variants. Yeast models of the five other variants showed very mild, if any, effect on mitochondrial function, suggesting that the variants do not have, at least alone, the potential to compromise human health.

Structures of annexin A2-PS DNA complexes show dominance of hydrophobic interactions in phosphorothioate binding.

The introduction of phosphorothioate (PS) linkages to the backbone of therapeutic nucleic acids substantially increases their stability and potency. It also affects their interactions with cellular proteins, but the molecular mechanisms that underlie this effect are poorly understood. Here, we report structural and biochemical studies of interactions between annexin A2, a protein that does not possess any known canonical DNA binding domains, and phosphorothioate-modified antisense oligonucleotides. We show that a unique mode of hydrophobic interactions between a sulfur atom of the phosphorothioate group and lysine and arginine residues account for the enhanced affinity of modified nucleic acid for the protein. Our results demonstrate that this mechanism of interaction is observed not only for nucleic acid-binding proteins but can also account for the association of PS oligonucleotides with other proteins. Using the anomalous diffraction of sulfur, we showed that preference for phosphorothioate stereoisomers is determined by the hydrophobic environment around the PS linkage that comes not only from protein but also from additional structural features within the ASO such as 5-Me groups on cytosine nucleobases.

Corrigendum: Ni2+-assisted hydrolysis may affect the human proteome; filaggrin degradation ex vivo as an example of possible consequences.

[This corrects the article DOI: 10.3389/fmolb.2022.828674.].

Competition between electrostatic interactions and halogen bonding in the protein-ligand system: structural and thermodynamic studies of 5,6-dibromobenzotriazole-hCK2α complexes.

CK2 is a member of the CMGC group of eukaryotic protein kinases and a cancer drug target. It can be efficiently inhibited by halogenated benzotriazoles and benzimidazoles. Depending on the scaffold, substitution pattern, and pH, these compounds are either neutral or anionic. Their binding poses are dictated by a hydrophobic effect (desolvation) and a tug of war between a salt bridge/hydrogen bond (to K68) and halogen bonding (to E114 and V116 backbone oxygens). Here, we test the idea that binding poses might be controllable by pH for ligands with near-neutral pKa, using the conditionally anionic 5,6-DBBt and constitutively anionic TBBt as our models. We characterize the binding by low-volume Differential Scanning Fluorimetry (nanoDSF), Isothermal Calorimetry (ITC), Hydrogen/Deuterium eXchange (HDX), and X-ray crystallography (MX). The data indicate that the ligand pose away from the hinge dominates for the entire tested pH range (5.5-8.5). The insensitivity of the binding mode to pH is attributed to the perturbation of ligand pKa upon binding that keeps it anionic in the ligand binding pocket at all tested pH values. However, a minor population of the ligand, detectable only by HDX, shifts towards the hinge in acidic conditions. Our findings demonstrate that electrostatic (ionic) interactions predominate over halogen bonding.

Isotope effects observed in diluted D2O/H2O mixtures identify HOD-induced low-density structures in D2O but not H2O.

Normal and heavy water are solvents most commonly used to study the isotope effect. The isotope effect of a solvent significantly influences the behavior of a single molecule in a solution, especially when there are interactions between the solvent and the solute. The influence of the isotope effect becomes more significant in D2O/H2O since the hydrogen bond in H2O is slightly weaker than its counterpart (deuterium bond) in D2O. Herein, we characterize the isotope effect in a mixture of normal and heavy water on the solvation of a HOD molecule. We show that the HOD molecule affects the proximal solvent molecules, and these disturbances are much more significant in heavy water than in normal water. Moreover, in D2O, we observe the formation of low-density structures indicative of an ordering of the solvent around the HOD molecule. The qualitative differences between HOD interaction with D2O and H2O were consistently confirmed with Raman spectroscopy and NMR diffusometry.

Improvement of native structure-based peptides as efficient inhibitors of protein-protein interactions of SARS-CoV-2 spike protein and human ACE2.

New pathogens responsible for novel human disease outbreaks in the last two decades are mainly the respiratory system viruses. Not different was the last pandemic episode, caused by infection of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). One of the extensively explored targets, in the recent scientific literature, as a possible way for rapid development of COVID-19 specific drug(s) is the interaction between the receptor-binding domain of the virus' spike (S) glycoprotein and human receptor angiotensin-converting enzyme 2 (hACE2). This protein-protein recognition process is involved in the early stages of the SARS-CoV-2 life cycle leading to the host cell membrane penetration. Thus, disrupting this interaction may block or significantly reduce the infection caused by the novel pathogen. Previously we have designed (by in silico structure-based analysis) three very short peptides having sequences inspirited by hACE2 native fragments, which effectively bind to the SARS-CoV-2 S protein and block its interaction with the human receptor. In continuation of the above mentioned studies, here we presented an application of molecular modeling approach resulting in improved binding affinity of the previously proposed ligand and its enhanced ability to inhibit meaningful host-virus protein-protein interaction. The new optimized hexapeptide binds to the virus protein with affinity one magnitude higher than the initial ligand and, as a very short peptide, has also great potential for further drug development. The peptide-based strategy is rapid and cost-effective for developing and optimizing efficient protein-protein interactions disruptors and may be successfully applied to discover antiviral candidates against other future emerging human viral infections.

Structural basis of transposon end recognition explains central features of Tn7 transposition systems.

Tn7 is a bacterial transposon with relatives containing element-encoded CRISPR-Cas systems mediating RNA-guided transposon insertion. Here, we present the 2.7 Å cryoelectron microscopy structure of prototypic Tn7 transposase TnsB interacting with the transposon end DNA. When TnsB interacts across repeating binding sites, it adopts a beads-on-a-string architecture, where the DNA-binding and catalytic domains are arranged in a tiled and intertwined fashion. The DNA-binding domains form few base-specific contacts leading to a binding preference that requires multiple weakly conserved sites at the appropriate spacing to achieve DNA sequence specificity. TnsB binding imparts differences in the global structure of the protein-bound DNA ends dictated by the spacing or overlap of binding sites explaining functional differences in the left and right ends of the element. We propose a model of the strand-transfer complex in which the terminal TnsB molecule is rearranged so that its catalytic domain is in a position conducive to transposition.

Variants in the pancreatic CUB and zona pellucida-like domains 1 (CUZD1) gene in early-onset chronic pancreatitis - A possible new susceptibility gene.

OBJECTIVE: Non-alcoholic chronic pancreatitis (NACP) frequently develops in the setting of genetic susceptibility associated with alterations in genes that are highly expressed in the pancreas. However, the genetic basis of NACP remains unresolved in a significant number of patients warranting a search for further risk genes. DESIGN: We analyzed CUZD1, which encodes the CUB and zona pellucida-like domains 1 protein that is found in high levels in pancreatic acinar cells. We sequenced the coding region in 1163 European patients and 2018 European controls. In addition, we analyzed 297 patients and 1070 controls from Japan. We analyzed secretion of wild-type and mutant CUZD1 from transfected cells using Western blotting. RESULTS: In the European cohort, we detected 30 non-synonymous variants. Using different prediction tools (SIFT, CADD, PROVEAN, PredictSNP) or the combination of these tools, we found accumulation of predicted deleterious variants in patients (p-value range 0.002-0.013; OR range 3.1-5.2). No association was found in the Japanese cohort, in which 13 non-synonymous variants were detected. Functional studies revealed >50% reduced secretion of 7 variants, however, these variants were not significantly enriched in European CP patients. CONCLUSION: Our data indicate that CUZD1 might be a novel susceptibility gene for NACP. How these variants predispose to pancreatitis remains to be elucidated.

The Rysto immune receptor recognises a broadly conserved feature of potyviral coat proteins.

Knowledge of the immune mechanisms responsible for viral recognition is critical for understanding durable disease resistance and successful crop protection. We determined how potato virus Y (PVY) coat protein (CP) is recognised by Rysto , a TNL immune receptor. We applied structural modelling, site-directed mutagenesis, transient overexpression, co-immunoprecipitation, infection assays and physiological cell death marker measurements to investigate the mechanism of Rysto -CP interaction. Rysto associates directly with PVY CP in planta that is conditioned by the presence of a CP central 149 amino acids domain. Each deletion that affects the CP core region impairs the ability of Rysto to trigger defence. Point mutations in the amino acid residues Ser125 , Arg157 , and Asp201 of the conserved RNA-binding pocket of potyviral CP reduce or abolish Rysto binding and Rysto -dependent responses, demonstrating that appropriate folding of the CP core is crucial for Rysto -mediated recognition. Rysto recognises the CPs of at least 10 crop-damaging viruses that share a similar core region. It confers immunity to plum pox virus and turnip mosaic virus in both Solanaceae and Brassicaceae systems, demonstrating potential utility in engineering virus resistance in various crops. Our findings shed new light on how R proteins detect different viruses by sensing conserved structural patterns.

The novel P330L pathogenic variant of aromatic amino acid decarboxylase maps on the catalytic flexible loop underlying its crucial role.

Aromatic amino acid decarboxylase (AADC) deficiency is a rare monogenic disease, often fatal in the first decade, causing severe intellectual disability, movement disorders and autonomic dysfunction. It is due to mutations in the gene coding for the AADC enzyme responsible for the synthesis of dopamine and serotonin. Using whole exome sequencing, we have identified a novel homozygous c.989C > T (p.Pro330Leu) variant of AADC causing AADC deficiency. Pro330 is part of an essential structural and functional element: the flexible catalytic loop suggested to cover the active site as a lid and properly position the catalytic residues. Our investigations provide evidence that Pro330 concurs in the achievement of an optimal catalytic competence. Through a combination of bioinformatic approaches, dynamic light scattering measurements, limited proteolysis experiments, spectroscopic and in solution analyses, we demonstrate that the substitution of Pro330 with Leu, although not determining gross conformational changes, results in an enzymatic species that is highly affected in catalysis with a decarboxylase catalytic efficiency decreased by 674- and 194-fold for the two aromatic substrates. This defect does not lead to active site structural disassembling, nor to the inability to bind the pyridoxal 5'-phosphate (PLP) cofactor. The molecular basis for the pathogenic effect of this variant is rather due to a mispositioning of the catalytically competent external aldimine intermediate, as corroborated by spectroscopic analyses and pH dependence of the kinetic parameters. Altogether, we determined the structural basis for the severity of the manifestation of AADC deficiency in this patient and discussed the rationale for a precision therapy.

Ni2+-Assisted Hydrolysis May Affect the Human Proteome; Filaggrin Degradation Ex Vivo as an Example of Possible Consequences.

Deficiency in a principal epidermal barrier protein, filaggrin (FLG), is associated with multiple allergic manifestations, including atopic dermatitis and contact allergy to nickel. Toxicity caused by dermal and respiratory exposures of the general population to nickel-containing objects and particles is a deleterious side effect of modern technologies. Its molecular mechanism may include the peptide bond hydrolysis in X1-S/T-c/p-H-c-X2 motifs by released Ni2+ ions. The goal of the study was to analyse the distribution of such cleavable motifs in the human proteome and examine FLG vulnerability of nickel hydrolysis. We performed a general bioinformatic study followed by biochemical and biological analysis of a single case, the FLG protein. FLG model peptides, the recombinant monomer domain human keratinocytes in vitro and human epidermis ex vivo were used. We also investigated if the products of filaggrin Ni2+-hydrolysis affect the activation profile of Langerhans cells. We found X1-S/T-c/p-H-c-X2 motifs in 40% of human proteins, with the highest abundance in those involved in the epidermal barrier function, including FLG. We confirmed the hydrolytic vulnerability and pH-dependent Ni2+-assisted cleavage of FLG-derived peptides and FLG monomer, using in vitro cell culture and ex-vivo epidermal sheets; the hydrolysis contributed to the pronounced reduction in FLG in all of the models studied. We also postulated that Ni-hydrolysis might dysregulate important immune responses. Ni2+-assisted cleavage of barrier proteins, including FLG, may contribute to clinical disease associated with nickel exposure.

cis-prenyltransferase 3 and α/ß-hydrolase are new determinants of dolichol accumulation in Arabidopsis.

Dolichols (Dols), ubiquitous components of living organisms, are indispensable for cell survival. In plants, as well as other eukaryotes, Dols are crucial for post-translational protein glycosylation, aberration of which leads to fatal metabolic disorders in humans and male sterility in plants. Until now, the mechanisms underlying Dol accumulation remain elusive. In this study, we have analysed the natural variation of the accumulation of Dols and six other isoprenoids among more than 120 Arabidopsis thaliana accessions. Subsequently, by combining QTL and GWAS approaches, we have identified several candidate genes involved in the accumulation of Dols, polyprenols, plastoquinone and phytosterols. The role of two genes implicated in the accumulation of major Dols in Arabidopsis-the AT2G17570 gene encoding a long searched for cis-prenyltransferase (CPT3) and the AT1G52460 gene encoding an α/ß-hydrolase-is experimentally confirmed. These data will help to generate Dol-enriched plants which might serve as a remedy for Dol-deficiency in humans.

5,6-diiodo-1H-benzotriazole: new TBBt analogue that minutely affects mitochondrial activity.

4,5,6,7-Tetrabromo-1H-benzotriazole is widely used as the reference ATP-competitive inhibitor of protein kinase CK2. Herein, we study its new analogs: 5,6-diiodo- and 5,6-diiodo-4,7-dibromo-1H-benzotriazole. We used biophysical (MST, ITC) and biochemical (enzymatic assay) methods to describe the interactions of halogenated benzotriazoles with the catalytic subunit of human protein kinase CK2 (hCK2α). To trace the biological activity, we measured their cytotoxicity against four reference cancer cell lines and the effect on the mitochondrial inner membrane potential. The results obtained lead to the conclusion that iodinated compounds are an attractive alternative to brominated ones. One of them retains the cytotoxicity against selected cancer cell lines of the reference TBBt with a smaller side effect on mitochondrial activity. Both iodinated compounds are candidate leaders in the further development of CK2 inhibitors.

Calmodulin as Ca2+-Dependent Interactor of FTO Dioxygenase.

FTO is an N6-methyladenosine demethylase removing methyl groups from nucleic acids. Several studies indicate the creation of FTO complexes with other proteins. Here, we looked for regulatory proteins recognizing parts of the FTO dioxygenase region. In the Calmodulin (CaM) Target Database, we found the FTO C-domain potentially binding CaM, and we proved this finding experimentally. The interaction was Ca2+-dependent but independent on FTO phosphorylation. We found that FTO-CaM interaction essentially influences calcium-binding loops in CaM, indicating the presence of two peptide populations-exchanging as CaM alone and differently, suggesting that only one part of CaM interacts with FTO, and the other one reminds free. The modeling of FTO-CaM interaction showed its stable structure when the half of the CaM molecule saturated with Ca2+ interacts with the FTO C-domain, whereas the other part is disconnected. The presented data indicate calmodulin as a new FTO interactor and support engagement of the FTO protein in calcium signaling pathways.

Synthesis of Novel Halogenated Heterocycles Based on o-Phenylenediamine and Their Interactions with the Catalytic Subunit of Protein Kinase CK2.

Protein kinase CK2 is a highly pleiotropic protein kinase capable of phosphorylating hundreds of protein substrates. It is involved in numerous cellular functions, including cell viability, apoptosis, cell proliferation and survival, angiogenesis, or ER-stress response. As CK2 activity is found perturbed in many pathological states, including cancers, it becomes an attractive target for the pharma. A large number of low-mass ATP-competitive inhibitors have already been developed, the majority of them halogenated. We tested the binding of six series of halogenated heterocyclic ligands derived from the commercially available 4,5-dihalo-benzene-1,2-diamines. These ligand series were selected to enable the separation of the scaffold effect from the hydrophobic interactions attributed directly to the presence of halogen atoms. In silico molecular docking was initially applied to test the capability of each ligand for binding at the ATP-binding site of CK2. HPLC-derived ligand hydrophobicity data are compared with the binding affinity assessed by low-volume differential scanning fluorimetry (nanoDSF). We identified three promising ligand scaffolds, two of which have not yet been described as CK2 inhibitors but may lead to potent CK2 kinase inhibitors. The inhibitory activity against CK2α and toxicity against four reference cell lines have been determined for eight compounds identified as the most promising in nanoDSF assay.