Valentina A. Lakunina

Targeting VDAC-bound hexokinase II: a promising approach for concomitant anti-cancer therapy

Introduction: Enhancement of glucose metabolism and repression of oxidative phosphorylation followed by the Warburg effect is the common hallmark of cancer cells. Hexokinase II (HKII) plays a dual role – first, HKII up-regulation results in increased glycolysis rates. Second, association of VDAC and HKII contributes to inhibition of apoptosis through repression of the formation of mitochondrial permeability transition pores. Areas covered: In this review, the role of HKII in evasion of apoptosis, aspects of HKII expression regulation, novel approaches targeting HKII and VDAC–HKII complexes and their application areas are discussed. Expert opinion: The dual role of HKII in cancer cells makes it an attractive target for anti-cancer therapy. Several agents, either synthetic or plant-derived, that target hexokinase and induce VDAC–HK complex dissociation have been identified to date. Targeting hexokinase, HK–VDAC complexes as well as other glycolytic proteins not only improves the efficacy of commonly used drugs. The most prominent benefit of this approach is the ability to overcome drug resistance, for example, to cisplatin or sorafenib. In some cases, it could create an insurmountable challenge for selection of appropriate therapy. Future studies and trials should address the issue of how to transfer these approaches into clinical practice.

Identification, Expression Analysis, and Target Prediction of Flax Genotroph MicroRNAs Under Normal and Nutrient Stress Conditions.

Cultivated flax (Linum usitatissimum L.) is an important plant valuable for industry. Some flax lines can undergo heritable phenotypic and genotypic changes (LIS-1 insertion being the most common) in response to nutrient stress and are called plastic lines. Offspring of plastic lines, which stably inherit the changes, are called genotrophs. MicroRNAs (miRNAs) are involved in a crucial regulatory mechanism of gene expression. They have previously been assumed to take part in nutrient stress response and can, therefore, participate in genotroph formation. In the present study, we performed high-throughput sequencing of small RNAs (sRNAs) extracted from flax plants grown under normal, phosphate deficient and nutrient excess conditions to identify miRNAs and evaluate their expression. Our analysis revealed expression of 96 conserved miRNAs from 21 families in flax. Moreover, 475 novel potential miRNAs were identified for the first time, and their targets were predicted. However, none of the identified miRNAs were transcribed from LIS-1. Expression of seven miRNAs (miR168, miR169, miR395, miR398, miR399, miR408, and lus-miR-N1) with up- or down-regulation under nutrient stress (on the basis of high-throughput sequencing data) was evaluated on extended sampling using qPCR. Reference gene search identified ETIF3H and ETIF3E genes as most suitable for this purpose. Down-regulation of novel potential lus-miR-N1 and up-regulation of conserved miR399 were revealed under the phosphate deficient conditions. In addition, the negative correlation of expression of lus-miR-N1 and its predicted target, ubiquitin-activating enzyme E1 gene, as well as, miR399 and its predicted target, ubiquitin-conjugating enzyme E2 gene, was observed. Thus, in our study, miRNAs expressed in flax plastic lines and genotrophs were identified and their expression and expression of their targets was evaluated using high-throughput sequencing and qPCR for the first time. These data provide new insights into nutrient stress response regulation in plastic flax cultivars.

Retrotransposon-based molecular markers for analysis of genetic diversity within the Genus Linum.

SSAP method was used to study the genetic diversity of 22 Linum species from sections Linum, Adenolinum, Dasylinum, Stellerolinum, and 46 flax cultivars. All the studied flax varieties were distinguished using SSAP for retrotransposons FL9 and FL11. Thus, the validity of SSAP method was demonstrated for flax marking, identification of accessions in genebank collections, and control during propagation of flax varieties. Polymorphism of Fl1a, Fl1b, and Cassandra insertions were very low in flax varieties, but these retrotransposons were successfully used for the investigation of Linum species. Species clusterization based on SSAP markers was in concordance with their taxonomic division into sections Dasylinum, Stellerolinum, Adenolinum, and Linum. All species of sect. Adenolinum clustered apart from species of sect. Linum. The data confirmed the accuracy of the separation in these sections. Members of section Linum are not as closely related as members of other sections, so taxonomic revision of this section is desirable. L. usitatissimum accessions genetically distant from modern flax cultivars were revealed in our work. These accessions are of utmost interest for flax breeding and introduction of new useful traits into flax cultivars. The chromosome localization of Cassandra retrotransposon in Linum species was determined.

Excess fertilizer responsive miRNAs revealed in Linum usitatissimum L.

Effective fertilizer application is necessary to increase crop yields and reduce risk of plant overdosing. It is known that expression level of microRNAs (miRNAs) alters in plants under different nutrient concentrations in soil. The aim of our study was to identify and characterize miRNAs with expression alterations under excessive fertilizer in agriculturally important crop - flax (Linum usitatissimum L.). We have sequenced small RNAs in flax grown under normal and excessive fertilizer using Illumina GAIIx. Over 14 million raw reads was obtained for two small RNA libraries. 84 conserved miRNAs from 20 families were identified. Differential expression was revealed for several flax miRNAs under excessive fertilizer according to high-throughput sequencing data. For 6 miRNA families (miR395, miR169, miR408, miR399, miR398 and miR168) expression level alterations were evaluated on the extended sampling using qPCR. Statistically significant up-regulation was revealed for miR395 under excessive fertilizer. It is known that target genes of miR395 are involved in sulfate uptake and assimilation. However, according to our data alterations of the expression level of miR395 could be associated not only with excess sulfur application, but also with redundancy of other macro- and micronutrients. Furthermore expression level was evaluated for miRNAs and their predicted targets. The negative correlation between miR399 expression and expression of its predicted target ubiquitin-conjugating enzyme E2 gene was shown in flax for the first time. So we suggested miR399 involvement in phosphate regulation in L. usitatissimum. Revealed in our study expression alterations contribute to miRNA role in flax response to excessive fertilizer.

Basal Glutathionylation of Na,K-ATPase α-Subunit Depends on Redox Status of Cells during the Enzyme Biosynthesis

Many viruses induce oxidative stress and cause S-glutathionylation of Cys residues of the host and viral proteins. Changes in cell functioning during viral infection may be associated with glutathionylation of a number of key proteins including Na,K-ATPase which creates a gradient of sodium and potassium ions. It was found that Na,K-ATPase α-subunit has a basal glutathionylation which is not abrogated by reducing agent. We have shown that acute hypoxia leads to increase of total glutathionylation level of Na,K-ATPase α-subunit; however, basal glutathionylation of α-subunit increases under prolonged hypoxia only. The role of basal glutathionylation in Na,K-ATPase function remains unclear. Understanding significance of basal glutathionylation is complicated by the fact that there are no X-ray structures of Na,K-ATPase with the identified glutathione molecules. We have analyzed all X-ray structures of the Na,K-ATPase α-subunit from pig kidney and found that there are a number of isolated cavities with unresolved electron density close to the relevant cysteine residues. Analysis of the structures showed that this unresolved density in the structure can be occupied by glutathione associated with cysteine residues. Here, we discuss the role of basal glutathionylation of Na,K-ATPase α-subunit and provide evidence supporting the view that this modification is cotranslational.

Direct interaction of beta-amyloid with Na,K-ATPase as a putative regulator of the enzyme function.

By maintaining the Na+ and K+ transmembrane gradient mammalian Na,K-ATPase acts as a key regulator of neuronal electrotonic properties. Na,K-ATPase has an important role in synaptic transmission and memory formation. Accumulation of beta-amyloid (Aβ) at the early stages of Alzheimer’s disease is accompanied by reduction of Na,K-ATPase functional activity. The molecular mechanism behind this phenomenon is not known. Here we show that the monomeric Aβ(1-42) forms a tight (Kd of 3 μM), enthalpy-driven equimolar complex with α1β1 Na,K-ATPase. The complex formation results in dose-dependent inhibition of the enzyme hydrolytic activity. The binding site of Aβ(1-42) is localized in the “gap” between the alpha- and beta-subunits of Na,K-ATPase, disrupting the enzyme functionality by preventing the subunits from shifting towards each other. Interaction of Na,K-ATPase with exogenous Aβ(1-42) leads to a pronounced decrease of the enzyme transport and hydrolytic activity and Src-kinase activation in neuroblastoma cells SH-SY5Y. This interaction allows regulation of Na,K-ATPase activity by short-term increase of the Aβ(1-42) level. However prolonged increase of Aβ(1-42) level under pathological conditions could lead to chronical inhibition of Na,K-ATPase and disruption of neuronal function. Taken together, our data suggest the role of beta-amyloid as a novel physiological regulator of Na,K-ATPase.

Cysteine residues 244 and 458–459 within the catalytic subunit of Na,K-ATPase control the enzyme's hydrolytic and signaling function under hypoxic conditions

Our previous findings suggested that reversible thiol modifications of cysteine residues within the actuator (AD) and nucleotide binding domain (NBD) of the Na,K-ATPase may represent a powerful regulatory mechanism conveying redox- and oxygen-sensitivity of this multifunctional enzyme. S-glutathionylation of Cys244 in the AD and Cys 454-458-459 in the NBD inhibited the enzyme and protected cysteines’ thiol groups from irreversible oxidation under hypoxic conditions. In this study mutagenesis approach was used to assess the role these cysteines play in regulation of the Na,K-ATPase hydrolytic and signaling functions. Several constructs of mouse α1 subunit of the Na,K-ATPase were produced in which Cys244, Cys 454-458-459 or Cys 244-454-458-459 were replaced by alanine. These constructs were expressed in human HEK293 cells. Non-transfected cells and those expressing murine α1 subunit were exposed to hypoxia or treated with oxidized glutathione (GSSG). Both conditions induced inhibition of the wild type Na,K-ATPase. Enzymes containing mutated mouse α1 lacking Cys244 or all four cysteines (Cys 244-454-458-459) were insensitive to hypoxia. Inhibitory effect of GSSG was observed for wild type murine Na,K-ATPase, but was less pronounced in Cys454-458-459Ala mutant and completely absent in the Cys244Ala and Cys 244-454-458-459Ala mutants. In cells, expressing wild type enzyme, ouabain induced activation of Src and Erk kinases under normoxic conditions, whereas under hypoxic conditions this effect was inversed. Cys454-458-459Ala substitution abolished Src kinase activation in response to ouabain treatment, uncoupled Src from Erk signaling, and interfered with O2-sensitivity of Na,K-ATPase signaling function. Moreover, modeling predicted that S-glutathionylation of Cys 458 and 459 should prevent inhibitory binding of Src to NBD. Our data indicate for the first time that cysteine residues within the AD and NBD influence hydrolytic as well as receptor function of the Na,K-ATPase and alter responses of the enzyme to hypoxia or upon treatment with cardiotonic steroids.

Changes in the receptor function of Na,K-ATPase during hypoxia and ischemia

Na,K-ATPase maintains sodium and potassium homeostasis. It is the only known receptor for cardiotonic steroids such as ouabain. Binding of ouabain to Na,K-ATPase leads to the activation of Src kinase and the subsequent initiation of intracellular signaling pathways, including the induction of apoptosis. Changes in Na,K-ATPase activity is one of the earliest responses to hypoxia and is most critical for cell survival. However, it is not known how the hypoxia affects the functioning of Na,K-ATPase as a receptor. We have shown that, under the conditions of hypoxia and ischemia, ouabain is less toxic for murine fibroblast cells (SC-1 cell line) and ouabain does not cause an increase in the level of reactive oxygen species, which is typically observed at 20% pO2. Under hypoxia, the treatment of cells with ouabain also does not lead to the activation of Na,K-ATPase-associated Src kinase. Thus, at low oxygen content, the receptor function of Na,K-ATPase is altered, and cells become less sensitive to cardiotonic steroids. The decrease in sensitivity to cardiotonic steroids, which is evident at hypoxic conditions, should be taken into account in clinical practice. At the same time, in the presence of ouabain the cells are less sensitive to hypoxia, which indicates that cardiotonic steroids can be protective in acute ischemia.

MethyMer: Design of combinations of specific primers for bisulfite sequencing of complete CpG islands

We present MethyMer, a Python-based tool aimed at selecting primers for amplification of complete CpG islands. These regions are difficult in terms of selecting appropriate primers because of their low-complexity, high GC content. Moreover, bisulfite treatment, in fact, leads to the reduction of the 4-letter alphabet (ATGC) to 3-letter one (ATG, except for methylated cytosines), and this also reduces region complexity and increases mispriming potential. MethyMer has a flexible scoring system, which optimizes the balance between various characteristics such as nucleotide composition, thermodynamic features (melting temperature, dimers [Formula: see text]G, etc.), the presence of CpG sites and polyN tracts, and primer specificity, which is assessed with aligning primers to the bisulfite-treated genome using bowtie (up to three mismatches are allowed). Users are able to customize desired or limit ranges of various parameters as well as penalties for non-desired values. Moreover, MethyMer allows picking up the optimal combination of PCR primer pairs to perform the amplification of a large genomic locus, e.g. CpG island or other hard-to-study region, with minimal overlap of the individual amplicons. MethyMer incorporates ENCODE genome annotation records (promoter/enhancer/insulator), The Cancer Genome Atlas (TCGA) CpG methylation data derived with Illumina Infinium 450K microarrays, and records on correlations between TCGA RNA-Seq and CpG methylation data for 20 cancer types. These databases are included in the MethyMer release. Our tool is available at https://sourceforge.net/projects/methymer/ .

Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes

BackgroundThe species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes.ResultsHigh-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH.Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14).ConclusionsHigh-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.