The association of protein-bound methionine sulfoxide with proteomic basis for aging in beech seeds.

BACKGROUND: European beech (Fagus sylvatica L.) trees produce seeds irregularly; therefore, it is necessary to store beech seeds for forestation. Despite the acquisition of desiccation tolerance during development, beech seeds are classified as intermediate because they lose viability during long-term storage faster than typical orthodox seeds. In this study, beech seeds stored for short (3 years) or long (20 years) periods under optimal conditions and displaying 92 and 30% germination capacity, respectively, were compared. RESULTS: Aged seeds displayed increased membrane damage, manifested as electrolyte leakage and lipid peroxidation levels. Analyses have been based on embryonic axes, which contained higher levels of reactive oxygen species (ROS) and higher levels of protein-bound methionine sulfoxide (MetO) in aged seeds. Using label-free quantitative proteomics, 3,949 proteins were identified, of which 2,442 were reliably quantified pointing to 24 more abundant proteins and 35 less abundant proteins in beech seeds under long-term storage conditions. Functional analyses based on gene ontology annotations revealed that nucleic acid binding activity (molecular function), ribosome organization or biogenesis and transmembrane transport (cellular processes), translational proteins (protein class) and membranous anatomical entities (cellular compartment) were affected in aged seeds. To verify whether MetO, the oxidative posttranslational modification of proteins that can be reversed via the action of methionine sulfoxide reductase (Msr) enzymes, is involved in the aging of beech seeds, we identified and quantified 226 MetO-containing proteins, among which 9 and 19 exhibited significantly up- and downregulated MetO levels, respectively, in beech seeds under long-term storage conditions. Several Msr isoforms were identified and recognized as MsrA1-like, MsrA4, MsrB5 and MsrB5-like in beech seeds. Only MsrA1-like displayed decreased abundance in aged seeds. CONCLUSIONS: We demonstrated that the loss of membrane integrity reflected in the elevated abundance of membrane proteins had a higher impact on seed aging progress than the MetO/Msr system. Proteome analyses enabled us to propose protein Sec61 and glyceraldehyde-3-phosphate dehydrogenase as potential longevity modulators in beech seeds.

Characterization and expression of a novel thaumatin-like protein (CcTLP1) from papaveraceous plant Corydalis cava.

Thaumatin-like proteins (TLPs, osmotins) form a protein family which shares a significant sequence homology to the sweet-tasting thaumatin from the plant Thaumatococcus daniellii. TLPs are not sweet-tasting and are involved in response to biotic stresses and developmental processes. Recently it has been shown using a proteomic approach that the tuber extract from Corydalis cava (Papaveraceae) contains a TLP protein. The aim of this work was to characterize the structure and expression of TLP from C. cava tubers. The results obtained using a PCR approach with degenerate primers demonstrated a coding sequence of a novel protein, named CcTLP1. It consists of 225 aa, has a predicted molecular weight of 24.2 kDa (NCBI GenBank accession no. KJ513303) and has 16 strictly conserved cysteine residues, which form 8 disulfide bridges and stabilize the 3D structure. CcTLP1 may be classified into class IX of plant TLPs. The highest CcTLP1 expression levels were shown by qPCR in the stem of the plant compared to other organs and in the medium-size plants compared to other growth phases. The results confirm that CcTLP1 is expressed during plant growth and development until flowering, with a possible defensive function against different stress conditions.

The OB-fold proteins of the Trypanosoma brucei editosome execute RNA-chaperone activity.

Sequence-deficient mitochondrial pre-mRNAs in African trypanosomes are substrates of a U-nucleotide-specific RNA editing reaction to generate translation-competent mRNAs. The reaction is catalyzed by a macromolecular protein complex termed the editosome. Editosomes execute RNA-chaperone activity to overcome the highly folded nature of pre-edited substrate mRNAs. The molecular basis for this activity is unknown. Here we test five of the OB-fold proteins of the Trypanosoma brucei editosome as candidates. We demonstrate that all proteins execute RNA-chaperone activity albeit to different degrees. We further show that the activities correlate to the surface areas of the proteins and we map the protein-induced RNA-structure changes using SHAPE-chemical probing. To provide a structural context for our findings we calculate a coarse-grained model of the editosome. The model has a shell-like structure: Structurally well-defined protein domains are separated from an outer shell of intrinsically disordered protein domains, which suggests a surface-driven mechanism for the chaperone activity.

Dimeric peroxiredoxins are druggable targets in human Burkitt lymphoma.

Burkitt lymphoma is a fast-growing tumor derived from germinal center B cells. It is mainly treated with aggressive chemotherapy, therefore novel therapeutic approaches are needed due to treatment toxicity and developing resistance. Disturbance of red-ox homeostasis has recently emerged as an efficient antitumor strategy. Peroxiredoxins (PRDXs) are thioredoxin-family antioxidant enzymes that scavenge cellular peroxides and contribute to red-ox homeostasis. PRDXs are robustly expressed in various malignancies and critically involved in cell proliferation, differentiation and apoptosis. To elucidate potential role of PRDXs in lymphoma, we studied their expression level in B cell-derived primary lymphoma cells as well as in cell lines. We found that PRDX1 and PRDX2 are upregulated in tumor B cells as compared with normal counterparts. Concomitant knockdown of PRDX1 and PRDX2 significantly attenuated the growth rate of lymphoma cells. Furthermore, in human Burkitt lymphoma cell lines, we isolated dimeric 2-cysteine peroxiredoxins as targets for SK053, a novel thiol-specific small-molecule peptidomimetic with antitumor activity. We observed that treatment of lymphoma cells with SK053 triggers formation of covalent PRDX dimers, accumulation of intracellular reactive oxygen species, phosphorylation of ERK1/2 and AKT and leads to cell cycle arrest and apoptosis. Based on site-directed mutagenesis and modeling studies, we propose a mechanism of SK053-mediated PRDX crosslinking, involving double thioalkylation of active site cysteine residues. Altogether, our results suggest that peroxiredoxins are novel therapeutic targets in Burkitt lymphoma and provide the basis for new approaches to the treatment of this disease.

Insight into the Meligethes aeneus voltage-sensitive sodium channel structure and an attempt to select the best pyrethroid ligands.

BACKGROUND: The voltage-sensitive sodium channel (VSSC) is a target for the pharmacological action of pyrethroids which are used in controlling pests, including those of agricultural importance. Among these is the pollen beetle (Meligethes aeneus F.) - the most serious pest of Brassica napus. Owing to the heavy use of pyrethroids, a widespread build-up of resistance has occurred. The main cause of pyrethroid insensitivity in M. aeneus is considered to be an increased oxidative metabolism; however, the additional mechanism of resistance associated with mutations in the VSSC might contribute to this phenomenon. RESULTS: We generated a VSSC 3D model to study the docking affinities of pyrethroids to their target site within the channel. Our goal was to identify the pyrethroids for which docking affinity scores were high and not affected by potential mutations in the VSSC. We found that the docking scores of cypermethrin are hardly influenced by the appearance of point mutations. Additionally, tau-fluvalinate, deltamethrin and bifenthrin are VSSC ligands with high affinity scores. CONCLUSIONS: Our docking models suggest that point mutations in the VSSC binding pocket might affect the stability of ligand interactions and change the pattern of ligand docking locations, which might have a potential effect on VSSC gating properties.

Structure and intrinsic disorder of the proteins of the Trypanosoma brucei editosome.

Mitochondrial pre-mRNAs in trypanosomatids undergo RNA editing to be converted into translatable mRNAs. The reaction is characterized by the insertion and deletion of uridine residues and is catalyzed by a macromolecular protein complex called the editosome. Despite intensive research, structural information for the majority of editosome proteins is still missing and no high resolution structure for the editosome exists. Here we present a comprehensive structural bioinformatics analysis of all proteins of the Trypanosoma brucei editosome. We specifically focus on the interplay between intrinsic order and disorder. According to computational predictions, editosome proteins involved in the basal reaction steps of the processing cycle are mostly ordered. By contrast, thirty percent of the amino acid content of the editosome is intrinsically disordered, which includes most prominently proteins with OB-fold domains. Based on the data we suggest a functional model, in which the structurally disordered domains of the complex are correlated with the RNA binding and RNA unfolding activity of the T. brucei editosome.

Identification of STAT1 and STAT3 specific inhibitors using comparative virtual screening and docking validation.

Signal transducers and activators of transcription (STATs) facilitate action of cytokines, growth factors and pathogens. STAT activation is mediated by a highly conserved SH2 domain, which interacts with phosphotyrosine motifs for specific STAT-receptor contacts and STAT dimerization. The active dimers induce gene transcription in the nucleus by binding to a specific DNA-response element in the promoter of target genes. Abnormal activation of STAT signaling pathways is implicated in many human diseases, like cancer, inflammation and auto-immunity. Searches for STAT-targeting compounds, exploring the phosphotyrosine (pTyr)-SH2 interaction site, yielded many small molecules for STAT3 but sparsely for other STATs. However, many of these inhibitors seem not STAT3-specific, thereby questioning the present modeling and selection strategies of SH2 domain-based STAT inhibitors. We generated new 3D structure models for all human (h)STATs and developed a comparative in silico docking strategy to obtain further insight into STAT-SH2 cross-binding specificity of a selection of previously identified STAT3 inhibitors. Indeed, by primarily targeting the highly conserved pTyr-SH2 binding pocket the majority of these compounds exhibited similar binding affinity and tendency scores for all STATs. By comparative screening of a natural product library we provided initial proof for the possibility to identify STAT1 as well as STAT3-specific inhibitors, introducing the 'STAT-comparative binding affinity value' and 'ligand binding pose variation' as selection criteria. In silico screening of a multi-million clean leads (CL) compound library for binding of all STATs, likewise identified potential specific inhibitors for STAT1 and STAT3 after docking validation. Based on comparative virtual screening and docking validation, we developed a novel STAT inhibitor screening tool that allows identification of specific STAT1 and STAT3 inhibitory compounds. This could increase our understanding of the functional role of these STATs in different diseases and benefit the clinical need for more drugable STAT inhibitors with high specificity, potency and excellent bioavailability.

Comparative screening and validation as a novel tool to identify STAT-specific inhibitors.

Signal transducers and activators of transcription (STATs) facilitate action of cytokines, growth factors and pathogens. STAT activation is mediated by a highly conserved SH2 domain, which interacts with phosphotyrosine (pTyr) motifs for specific STAT-receptor contacts and STAT dimerization. The active dimers induce gene transcription in the nucleus by binding to specific DNA-response elements of target genes. Abnormal activation of STAT signaling pathways is implicated in many human diseases, like cancer, inflammation and auto-immunity. STAT inhibitory strategies mostly focus on inhibiting STAT dimerization using small molecules identified by molecular modeling, virtual or library screening, or natural products. Searches for STAT-targeting compounds, exploring the pTyr-SH2 interaction area, yielded many small molecules for STAT3 but sparsely for other STATs. So far, no STAT-targeting drug is approved by the FDA. Moreover, many of these inhibitors do not seem STAT-specific, thereby questioning the present selection strategies of SH2 domain-based STAT inhibitors. This illustrates the need for better models, and screening and validation tools for more druggable STAT inhibitors with high specificity, potency and excellent bioavailability. Based on newly developed 3D structure models for all human (h)STATs, we propose a pipeline approach that combines comparative in silico docking of STAT-SH2 models with an in vitro STAT phosphorylation assay, as a novel tool to screen multi-million compound libraries and identify specific inhibitors for different STATs. Identification of specific and effective STAT inhibitory compounds could provide a tool to increase our understanding of their functional role in different diseases, and serve as therapeutic strategies in cancer, inflammation and auto-immunity.

Structural and degradative aspects of ornithine decarboxylase antizyme inhibitor 2.

Ornithine decarboxylase (ODC) is the key enzyme in the polyamine biosynthetic pathway. ODC levels are controlled by polyamines through the induction of antizymes (AZs), small proteins that inhibit ODC and target it to proteasomal degradation without ubiquitination. Antizyme inhibitors (AZIN1 and AZIN2) are proteins homologous to ODC that bind to AZs and counteract their negative effect on ODC. Whereas ODC and AZIN1 are well-characterized proteins, little is known on the structure and stability of AZIN2, the lastly discovered member of this regulatory circuit. In this work we first analyzed structural aspects of AZIN2 by combining biochemical and computational approaches. We demonstrated that AZIN2, in contrast to ODC, does not form homodimers, although the predicted tertiary structure of the AZIN2 monomer was similar to that of ODC. Furthermore, we identified conserved residues in the antizyme-binding element, whose substitution drastically affected the capacity of AZIN2 to bind AZ1. On the other hand, we also found that AZIN2 is much more labile than ODC, but it is highly stabilized by its binding to AZs. Interestingly, the administration of the proteasome inhibitor MG132 caused differential effects on the three AZ-binding proteins, having no effect on ODC, preventing the degradation of AZIN1, but unexpectedly increasing the degradation of AZIN2. Inhibitors of the lysosomal function partially prevented the effect of MG132 on AZIN2. These results suggest that the degradation of AZIN2 could be also mediated by an alternative route to that of proteasome. These findings provide new relevant information on this unique regulatory mechanism of polyamine metabolism.

In silico simulations of STAT1 and STAT3 inhibitors predict SH2 domain cross-binding specificity.

Signal transducers and activators of transcription (STATs) comprise a family of transcription factors that are structurally related and which participate in signaling pathways activated by cytokines, growth factors and pathogens. Activation of STAT proteins is mediated by the highly conserved Src homology 2 (SH2) domain, which interacts with phosphotyrosine motifs for specific contacts between STATs and receptors and for STAT dimerization. By generating new models for human (h)STAT1, hSTAT2 and hSTAT3 we applied comparative in silico docking to determine SH2-binding specificity of the STAT3 inhibitor stattic, and of fludarabine (STAT1 inhibitor). Thus, we provide evidence that by primarily targeting the highly conserved phosphotyrosine (pY+0) SH2 binding pocket stattic is not a specific hSTAT3 inhibitor, but is equally effective towards hSTAT1 and hSTAT2. This was confirmed in Human Micro-vascular Endothelial Cells (HMECs) in vitro, in which stattic inhibited interferon-α-induced phosphorylation of all three STATs. Likewise, fludarabine inhibits both hSTAT1 and hSTAT3 phosphorylation, but not hSTAT2, by competing with the highly conserved pY+0 and pY-X binding sites, which are less well-preserved in hSTAT2. Moreover we observed that in HMECs in vitro fludarabine inhibits cytokine and lipopolysaccharide-induced phosphorylation of hSTAT1 and hSTAT3 but does not affect hSTAT2. Finally, multiple sequence alignment of STAT-SH2 domain sequences confirmed high conservation between hSTAT1 and hSTAT3, but not hSTAT2, with respect to stattic and fludarabine binding sites. Together our data offer a molecular basis that explains STAT cross-binding specificity of stattic and fludarabine, thereby questioning the present selection strategies of SH2 domain-based competitive small inhibitors.

Ratio of mutated versus wild-type coat protein sequences in Pepino mosaic virus determines the nature and severity of yellowing symptoms on tomato plants.

Recently, Pepino mosaic virus (PepMV) infections causing severe yellowing symptoms in tomato plants have been reported in glasshouse tomato crops. When studying this phenomenon in commercial glasshouses, two different types of yellowing symptoms, occurring in adjacent plants, were distinguished: interveinal leaf yellowing and yellow mosaics. After several weeks, the interveinal leaf yellowing symptoms gradually disappeared and the plant heads became green again, with yellow mosaic patterns on the leaves as an intermediate stage. The sequencing of multiple isolates causing interveinal leaf yellowing identified two point mutations, occurring in positions 155 and 166 of the coat protein (CP), as unique to the yellowing pathotype. Site-directed mutagenesis of infectious clones confirmed that both CP mutations are determinants of the interveinal leaf yellowing symptoms. Sequencing of CP clones from plants or plant parts with the yellow mosaic symptoms resulted in a mixture of wild-type and mutated sequences, whereas sequencing of CP clones from the green heads of recovered plants resulted in only wild-type sequences. Yellow mosaic symptoms could be reproduced by inoculation of an artificial 1:1 mixture of RNA transcripts from the wild-type and mutated infectious clones. These results show that the ratio of mutated versus wild-type sequences can determine the nature and severity of symptom development. The gradual recovery of the plants, which coincides with the disappearance of the yellowing mutations, suggests that selection pressure acts to the advantage of the wild-type virus. Experiments with wild-type and mutated infectious clones showed that reverse mutation events from mutant to wild-type occur and that the wild-type virus does not have a replicative advantage over the mutant. These results suggest that reverse mutation events occur, with subsequent selection pressure acting in favour of the wild-type virus in the growing plant parts, possibly related to a lower long-distance movement efficiency of the mutant.

Identification of a Coding Sequence and Structure Modeling of a Glycine-Rich RNA-Binding Protein (CmGRP1) from Chelidonium majus L.

The family of glycine-rich plant proteins (GRPs) is a large and complex group of proteins that share, as a common feature, the presence of glycine-rich domains arranged in (Gly)n-X repeats that are suggested to be involved in protein-protein interactions, RNA binding, and nucleolar targeting. These proteins are implicated in several independent physiological processes. Some are components of cell walls of many higher plants, while others are involved in molecular responses to environmental stress, and mediated by post-transcriptional regulatory mechanisms. The goals of this study are to identify the coding sequence of a novel glycine-rich RNA-binding protein from Chelidonium majus and to propose its structural model. DNA fragments obtained using degenerate PCR primers showed high sequence identities with glycine-rich RNA-binding protein coding sequences from different plant species. A 439-bp nucleotide sequence is identified coding for a novel polypeptide composed of 146 amino acids, designated as CmGRP1 (C. majus glycine-rich protein 1), with a calculated MW of 14,931 Da (NCBI GenBank accession no. HM173636). Using NCBI CDD and GeneSilico MetaServer, a single conserved domain, the RNA recognition motif (RRM), was detected in CmGRP1. The C-terminal region of CmGRP1 is a glycine-rich motif (GGGGxxGxGGGxxG), and it is predicted to be disordered. Based on a 1fxl crystal structure, a 3D model of CmGRP1 is proposed. CmGRP1 can be classified as a class IVa plant GRP, implicated to play a role in plant defense.

Molecular evolution of dihydrouridine synthases.

BACKGROUND: Dihydrouridine (D) is a modified base found in conserved positions in the D-loop of tRNA in Bacteria, Eukaryota, and some Archaea. Despite the abundant occurrence of D, little is known about its biochemical roles in mediating tRNA function. It is assumed that D may destabilize the structure of tRNA and thus enhance its conformational flexibility. D is generated post-transcriptionally by the reduction of the 5,6-double bond of a uridine residue in RNA transcripts. The reaction is carried out by dihydrouridine synthases (DUS). DUS constitute a conserved family of enzymes encoded by the orthologous gene family COG0042. In protein sequence databases, members of COG0042 are typically annotated as "predicted TIM-barrel enzymes, possibly dehydrogenases, nifR3 family". RESULTS: To elucidate sequence-structure-function relationships in the DUS family, a comprehensive bioinformatic analysis was carried out. We performed extensive database searches to identify all members of the currently known DUS family, followed by clustering analysis to subdivide it into subfamilies of closely related sequences. We analyzed phylogenetic distributions of all members of the DUS family and inferred the evolutionary tree, which suggested a scenario for the evolutionary origin of dihydrouridine-forming enzymes. For a human representative of the DUS family, the hDus2 protein suggested as a potential drug target in cancer, we generated a homology model. While this article was under review, a crystal structure of a DUS representative has been published, giving us an opportunity to validate the model. CONCLUSIONS: We compared sequences and phylogenetic distributions of all members of the DUS family and inferred the phylogenetic tree, which provides a framework to study the functional differences among these proteins and suggests a scenario for the evolutionary origin of dihydrouridine formation. Our evolutionary and structural classification of the DUS family provides a background to study functional differences among these proteins that will guide experimental analyses.

Analysis of two strains of Peanut stunt virus: satRNA-associated and satRNA free.

Peanut stunt virus (PSV) is a pathogen of legumes, vegetables, trees, and weeds occurring worldwide. The species is characterized by significant genetic variability. PSV strains are classified into four subgroups on the basis of their nucleotide sequence homology. Here, we are presenting two further, fully sequenced PSV strains-PSV-Ag and PSV-G, that could be considered as I subgroup representatives. However, their sequence homology with other typical I subgroups members, similarly as another strain-PSV-P, characterized by our group previously, is lower than 90%. This lead us to propose further subdivision of the I subgroup into IA, IB, and IC units, and to classify PSV-Ag and PSV-G strains to the last one. In this article, we are showing that identity level of PSV-Ag and PSV-G is very high and apart from the presence of satRNA in the first one, they differ only by a few nucleotides in their genomic RNAs. Nevertheless, symptoms they cause on host plants might differ significantly, just as the levels in infected plants. Effect of single amino acid changes between strains on the three-dimensional structure of viral proteins was analyzed. Differences occur mainly on the protein surfaces which can possibly affect protein-protein interaction in infected cells, which is discussed.

Identification and modeling of a phosphatase-like domain in a tRNA 2'-O-ribosyl phosphate transferase Rit1p.

Cytoplasmic initiator tRNAs from plants and fungi are excluded from participating in translational elongation by the presence of a unique 2'-phosphoribosyl modification of purine 64, introduced posttranscriptionally by the enzyme Rit1p. Members of the Rit1p family show no obvious similarity to other proteins or domains, there is no structural information available to guide experimental analyses, and the mechanism of action of this enzyme remains a mystery. Using protein fold recognition, we identified a phosphatase-like domain in the C-terminal part of Rit1p. A comparative model of the C-terminal domain was constructed and used to predict the function of conserved residues and to propose the mechanism of action of Rit1p. The model will facilitate experimental analyses of Rit1p and its interactions with the initiator tRNA substrate.

STAT1 as a novel therapeutical target in pro-atherogenic signal integration of IFNγ, TLR4 and IL-6 in vascular disease.

Inflammation participates importantly in host defenses against infectious agents and injury, but it also contributes to the pathophysiology of atherosclerosis. Recruitment of blood leukocytes to the injured vascular endothelium characterizes the initiation and progression of atherosclerosis and involves many inflammatory mediators, modulated by cells of both innate and adaptive immunity. The pro-inflammatory cytokine, interferon (IFN)-γ derived from T cells, is vital for both innate and adaptive immunity and is also expressed at high levels in atherosclerotic lesions. As such IFN-γ plays a crucial role in the pathology of atherosclerosis through activation of signal transducer and activator of transcription (STAT) 1. Toll-like receptors (TLRs) are innate immune pattern recognition receptors (PRRs) expressed on a variety of cells, and thus initiate and sustain the inflammatory response in atherosclerosis. More recent studies have revealed that STAT1 is involved in the signaling events mediated by TLR4, leading to increased expression of several pro-inflammatory and pro-atherogenic mediators. By upregulating members of the Suppressors Of Cytokine Signaling (SOCS) family that regulate cellular responsiveness to immune signals, IFNγ and TLR4-activated pathways have also shown to inhibit IL-6 STAT3-dependent anti-inflammatory signaling and potentially shift IL-6 to a STAT1 activating pro-inflammatory cytokine. Consequently, STAT1 has been identified as a point of convergence for the cross-talk between the pro-atherogenic IFN-γ, TLR4 and IL-6 activated pathways in immune as well as vascular cells, as such amplifying pro-inflammatory signals. This results in augmented smooth muscle cell (SMC) and leukocyte migration, leukocyte to endothelial cell (EC) adhesion and foam cell formation, and could encompass a novel mechanism involved in the initiation and progression of atherosclerosis. Therefore, application of small inhibitory compounds that specifically interact with the SH2-phosphotyrosine pocket of STAT1, proposed here as a novel working mechanism for the known STAT1 inhibitor fludarabine, could be a promising tool in the development of a therapeutical strategy for atherosclerosis.

Tridimensional model structure and patterns of molecular evolution of Pepino mosaic virus TGBp3 protein.

BACKGROUND: Pepino mosaic virus (PepMV) is considered one of the most dangerous pathogens infecting tomatoes worldwide. The virus is highly diverse and four distinct genotypes, as well as inter-strain recombinants, have already been described. The isolates display a wide range on symptoms on infected plant species, ranging from mild mosaic to severe necrosis. However, little is known about the mechanisms and pattern of PepMV molecular evolution and about the role of individual proteins in host-pathogen interactions. METHODS: The nucleotide sequences of the triple gene block 3 (TGB3) from PepMV isolates varying in symptomatology and geographic origin have been analyzed. The modes and patterns of molecular evolution of the TGBp3 protein were investigated by evaluating the selective constraints to which particular amino acid residues have been subjected during the course of diversification. The tridimensional structure of TGBp3 protein has been modeled de novo using the Rosetta algorithm. The correlation between symptoms development and location of specific amino acids residues was analyzed. RESULTS: The results have shown that TGBp3 has been evolving mainly under the action of purifying selection operating on several amino acid sites, thus highlighting its functional role during PepMV infection. Interestingly, amino acid 67, which has been previously shown to be a necrosis determinant, was found to be under positive selection. CONCLUSIONS: Identification of diverse selection events in TGB3p3 will help unraveling its biological functions and is essential to an understanding of the evolutionary constraints exerted on the Potexvirus genome. The estimated tridimensional structure of TGBp3 will serve as a platform for further sequence, structural and function analysis and will stimulate new experimental advances.

MODOMICS: a database of RNA modification pathways. 2008 update.

MODOMICS, a database devoted to the systems biology of RNA modification, has been subjected to substantial improvements. It provides comprehensive information on the chemical structure of modified nucleosides, pathways of their biosynthesis, sequences of RNAs containing these modifications and RNA-modifying enzymes. MODOMICS also provides cross-references to other databases and to literature. In addition to the previously available manually curated tRNA sequences from a few model organisms, we have now included additional tRNAs and rRNAs, and all RNAs with 3D structures in the Nucleic Acid Database, in which modified nucleosides are present. In total, 3460 modified bases in RNA sequences of different organisms have been annotated. New RNA-modifying enzymes have been also added. The current collection of enzymes includes mainly proteins for the model organisms Escherichia coli and Saccharomyces cerevisiae, and is currently being expanded to include proteins from other organisms, in particular Archaea and Homo sapiens. For enzymes with known structures, links are provided to the corresponding Protein Data Bank entries, while for many others homology models have been created. Many new options for database searching and querying have been included. MODOMICS can be accessed at http://genesilico.pl/modomics.

The sequence and model structure analysis of three Polish peanut stunt virus strains.

Peanut stunt virus (PSV) belongs to the Cucumovirus genus of the family Bromoviridae and is widely distributed worldwide, also in Poland. PSV is a common pathogen of a wide range of economically important plants. Its coat protein (CP), similarly as in other viruses, plays an important role in many processes during viral life cycle and has great impact on the infectivity. In this study, we present the results of sequence-structure analysis of CP derived from three Polish strains of PSV: PSV Ag, G, and P. Sequences were determined using RT-PCR amplification followed by sequencing and compared with each other and also with CP from other known PSV viruses. We analyzed their phylogenetic relationship, based on CP sequence, using bioinformatic tools as well as their spatial model using homology-modeling approach with combination of ROSETTA algorithm for de novo modeling. We compared our model with those recently obtained for other cucumoviruses including PSV-Er. Our results have shown that all Polish strains probably belong to the first subgroup of PSV viruses. Homology level between strains Ag and G proved very high. Using theoretical modeling approach we obtained a model very similar to the one resolved previously with the differences caused by slightly different amino acid sequence. We have also undertaken an attempt to analyze its distant regions; however, results are not unequivocal. Analysis of symptoms and their correlation with specific amino acid position was also performed on the basis of results published elsewhere. The definite interpretation is impeded by the presence of satellite RNAs in Ag and P strains modulating symptoms' severity, though.

Structural bioinformatics analysis of enzymes involved in the biosynthesis pathway of the hypermodified nucleoside ms(2)io(6)A37 in tRNA.

TRNAs from all organisms contain posttranscriptionally modified nucleosides, which are derived from the four canonical nucleosides. In most tRNAs that read codons beginning with U, adenosine in the position 37 adjacent to the 3' position of the anticodon is modified to N(6)-(Delta(2)-isopentenyl) adenosine (i(6)A). In many bacteria, such as Escherichia coli, this residue is typically hypermodified to N(6)-isopentenyl-2-thiomethyladenosine (ms(2)i(6)A). In a few bacteria, such as Salmonella typhimurium, ms(2)i(6)A can be further hydroxylated to N(6)-(cis-4-hydroxyisopentenyl)-2-thiomethyladenosine (ms(2)io(6)A). Although the enzymes that introduce the respective modifications (prenyltransferase MiaA, methylthiotransferase MiaB, and hydroxylase MiaE) have been identified, their structures remain unknown and sequence-function relationships remain obscure. We carried out sequence analysis and structure prediction of MiaA, MiaB, and MiaE, using the protein fold-recognition approach. Three-dimensional models of all three proteins were then built using a new modeling protocol designed to overcome uncertainties in the alignments and divergence between the templates. For MiaA and MiaB, the catalytic core was built based on the templates from the P-loop NTPase and Radical-SAM superfamilies, respectively. For MiaB, we have also modeled the C-terminal TRAM domain and the newly predicted N-terminal flavodoxin-fold domain. For MiaE, we confidently predict that it shares the three-dimensional fold with the ferritin-like four-helix bundle proteins and that it has a similar active site and mechanism of action to diiron carboxylate enzymes, in particular, methane monooxygenase (E.C.1.14.13.25) that catalyses the biological hydroxylation of alkanes. Our models provide the first structural platform for enzymes involved in the biosynthesis of i(6)A, ms(2)i(6)A, and ms(2)io(6)A, explain the data available from the literature and will help to design further experiments and interpret their results.