The PROSITE database, its status in 1999

The structural insights of L-asparaginase from Pseudomonas aeruginosa CSPS4 at elevated temperatures highlight its thermophilic nature

In the present investigation, a novel thermophilic L-asparaginase (Asn_PA) from Pseudomonas aeruginosa CSPS4 was investigated to explore its structural insights at elevated temperatures. Sequence analysis of Asn_PA depicted three conserved motifs (VVILATGGTIAG, DGIVITHGTDTLEETAYFL, and, LRKQGVQIIRSSHVNAGGF), of them, two motifs exhibit catalytically-important residues i.e., T45 and T125. A homology modelling-based structure model for Asn_PA was generated with 4PGA as the top-matched template. The predicted structure was validated and energy was minimized. Molecular docking was carried out cantered at the active site for asparagine and glutamine as its substrate ligands. The enzyme-substrate interaction analysis showed binding affinities of  - 4.8 and  - 4.1 kcal/mol for asparagine and glutamine respectively. Molecular dynamics (MD) simulation studies showed a better stability of Asn_PA at temperatures of 60 °C, over 40, 50 and, 80 °C, making this enzyme a novel L-asparaginase from other mesophilic P. aeruginosa strain. The trajectory analysis showed that RMSD, Rg, and, SASA values correlate well with each other in the different tested temperatures during the MD analysis. Thus, the present findings encourage extensive characterization of the Asn_PA using laboratory experiments to understand the structural behavior of the active site loop in an open or closed state with and without the substrate molecules.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04072-w.

Detection of Secreted Effector Proteins from Phelipanche aegyptiaca During Invasion of Melon Roots

Parasites can interact with their host plants through the induction and delivery of secreted effector proteins that facilitate plant colonization by decomposing plant cell walls and inhibiting plant immune response to weaken the defense ability of the host. Yet effectors mediating parasitic plant-host interactions are poorly understood. Phelipanche aegyptiaca is an obligate root parasite plant causing severe yield and economic losses in agricultural fields worldwide. Host resistance against P. aegyptiaca occurred during the attachment period of parasitism. Comparative transcriptomics was used to assess resistant and susceptible interactions simultaneously between P. aegyptiaca and two contrasting melon cultivars. In total, 2,740 secreted proteins from P. aegyptiaca were identified here. Combined with transcriptome profiling, 209 candidate secreted effector proteins (CSEPs) were predicted, with functional annotations such as cell wall degrading enzymes, protease inhibitors, transferases, kinases, and elicitor proteins. A heterogeneous expression system in Nicotiana benthamiana was used to investigate the functions of 20 putatively effector genes among the CSEPs. Cluster 15140.0 can suppress BAX-triggered programmed cell death in N. benthamiana. These findings showed that the prediction of P. aegyptiaca effector proteins based on transcriptomic analysis and multiple bioinformatics software is effective and more accurate, providing insights into understanding the essential molecular nature of effectors and laying the foundation of revealing the parasite mechanism of P. aegyptiaca, which is helpful in understanding parasite-host plant interaction.

Protein Myristoylation Plays a Role in the Nuclear Entry of the Parvovirus Minute Virus of Mice

Being nonpathogenic to humans, rodent parvoviruses (PVs) are naturally oncolytic viruses with great potential as anti-cancer agents. As these viruses replicate in the host cell nucleus, they must gain access to the nucleus during infection. The PV minute virus of mice (MVM) and several other PVs transiently disrupt the nuclear envelope (NE) and enter the nucleus through the resulting breaks. However, the molecular basis of this unique nuclear entry pathway remains uncharacterized. In this study, we used MVM as a model to investigate the molecular mechanism by which PVs induce NE disruption during viral nuclear entry. By combining bioinformatics analyses, metabolic labeling assays, mutagenesis, and pharmacological inhibition, we identified a functional myristoylation site at the sequence ⁷⁸GGKVGH⁸³ of the unique portion of the capsid protein VP1 (VP1u) of MVM. Performing proteolytic cleavage studies with a peptide containing this myristoylation site or with purified virions, we found tryptophan at position 77 of MVM VP1u is susceptible to chymotrypsin cleavage, implying this cleavage exposes G (glycine) 78 at the N-terminus of VP1u for myristoylation. Subsequent experiments using inhibitors of myristoylation and cellular proteases with MVM-infected cells, or an imaging-based quantitative NE permeabilization assay, further indicate protein myristoylation and a chymotrypsin-like activity are essential for MVM to locally disrupt the NE during viral nuclear entry. We thus propose a model for the nuclear entry of MVM in which NE disruption is mediated by VP1u myristoylation after the intact capsid undergoes proteolytic processing to expose the required N-terminal G for myristoylation. IMPORTANCE Rodent parvoviruses (PVs), including minute virus of mice (MVM), have the ability to infect and kill cancer cells and thereby possess great potential in anti-cancer therapy. In fact, some of these viruses are currently being investigated in both preclinical studies and clinical trials to treat a wide variety of cancers. However, the detailed mechanism of how PVs enter the cell nucleus remains unknown. In this study, we for the first time demonstrated a chemical modification called "myristoylation" of a MVM protein plays an essential role in the nuclear entry of the virus. We also showed, in addition to protein myristoylation, a chymotrypsin-like activity, which may come from cellular proteasomes, is required for MVM to get myristoylated and enter the nucleus. These findings deepen our understanding on how MVM and other related PVs infect host cells and provide new insights for the development of PV-based anti-cancer therapies.

Production of a Novel Marine Pseudomonas aeruginosa Recombinant L-Asparaginase: Insight on the Structure and Biochemical Characterization

The present study focused on the cloning, expression, and characterization of L-asparaginase of marine Pseudomonas aeruginosa HR03 isolated from fish intestine. Thus, a gene fragment containing the L-asparaginase sequence of Pseudomonas aeruginosa HR03 isolated from the fish intestine was cloned in the pET21a vector and then expressed in Escherichia coli BL21 (DE3) cells. Thereafter, the recombinant L-asparaginase (HR03Asnase) was purified by nickel affinity chromatography, and the enzymatic properties of HR03Asnase, including the effects of pH and temperature on HR03Asnase activity and its kinetic parameters, were determined. The recombinant enzyme HR03Asnase showed the highest similarity to type I L-asparaginase from Pseudomonas aeruginosa. The three-dimensional (3D) modeling results indicate that HR03Asnase exists as a homotetramer. Its molecular weight was 35 kDa, and the maximum activity of the purified enzyme was observed at pH8 and at 40 °C. The k_m and V_max of the enzyme obtained with L-asparagine as substrate were 10.904 mM and 3.44 × 10^-2 mM/min, respectively. The maximum activity of HR03Asnase was reduced by 50% at 90 °C after 10-min incubation; however, the enzyme maintained more than 20% of its activity after 30-min incubation. This enzyme also maintained almost 50% of its activity at pH 12 after 40-min incubation. The evaluation of pH and temperature stability of HR03Asnase showed that the enzyme has a wide range of activity, which is a suitable characteristic for its application in different industries. Overall, the results of the present study indicate that marine sources are promising biological reservoirs for enzymes to be used for biotechnological purposes, and marine thermostable HR03Asnase is likely a potential candidate for its future usage in the pharmaceutical and food industries.

High Molecular Weight Kininogen: A Review of the Structural Literature

Kininogens are multidomain glycoproteins found in the blood of most vertebrates. High molecular weight kininogen demonstrate both carrier and co-factor activity as part of the intrinsic pathway of coagulation, leading to thrombin generation. Kininogens are the source of the vasoactive nonapeptide bradykinin. To date, attempts to crystallize kininogen have failed, and very little is known about the shape of kininogen at an atomic level. New advancements in the field of cryo-electron microscopy (cryoEM) have enabled researchers to crack the structure of proteins that has been refractory to traditional crystallography techniques. High molecular weight kininogen is a good candidate for structural investigation by cryoEM. The goal of this review is to summarize the findings of kininogen structural studies.

The Arac-like transcriptional regulator YqhC is involved in pathogenicity of Erwinia amylovora

AIMS

This study aimed to identify virulence-associated genes and functions that affect disease development on pear caused by Erwinia amylovora EaUMG3 isolated from Iran.

METHODS AND RESULTS

A mini-Tn5 transposon library was generated in EaUMG3. An E. amylovora mutant that had lost its ability to cause lesions on immature pear fruits, was selected for further analysis. This mutant was shown to have a transposon insertion in yqhC, a gene belongs to the AraC family of transcriptional regulators. A mutant of the wild-type EaUMG3 carrying an unmarked deletion of the yqhC gene was created using pDMS197. The Ea∆yqhC mutant showed reduced disease progression on immature pear fruits and pear plants, reduced motility and significantly lower levels of the virulence factors siderophore and amylovoran. Complementation with yqhC cloned in pBBR1MCS restored disease progression and the level of virulence factors to near wild type.

CONCLUSION

YqhC transcriptional regulator is necessary for full virulence of E. amylovora. In addition, this regulator affects virulence factors such as siderophore production, amylovoran production, and motility.

SIGNIFICANCE AND IMPACT OF STUDY

The identification of a novel transcriptional regulator with strong impact in the pathogenesis of E. amylovora, an organism causing significant economic losses in fruit production.

SHR4z, a novel decoy effector from the haustorium of the parasitic weed Striga gesnerioides, suppresses host plant immunity

Cowpea (Vigna unguiculata) cultivar B301 is resistant to races SG4 and SG3 of the root parasitic weed Striga gesnerioides, developing a hypersensitive response (HR) at the site of parasite attachment. By contrast, race SG4z overcomes B301 resistance and successfully parasitises the plant. Comparative transcriptomics and in silico analysis identified a small secreted effector protein dubbed Suppressor of Host Resistance 4z (SHR4z) in the SG4z haustorium that upon transfer to the host roots causes a loss of host immunity (i.e. decreased HR and increased parasite growth). SHR4z has significant homology to the short leucine-rich repeat (LRR) domain of SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) family proteins and functions by binding to VuPOB1, a host BTB-BACK domain-containing ubiquitin E3 ligase homologue, leading to its rapid turnover. VuPOB1 is shown to be a positive regulator of HR since silencing of VuPOB1 expression in transgenic B301 roots lowers the frequency of HR and increases the levels of successful SG4 parasitism and overexpression decreases parasitism by SG4z. These findings provide new insights into how parasitic weeds overcome host defences and could potentially contribute to the development of novel strategies for controlling Striga and other parasitic weeds thereby enhancing crop productivity and food security globally.

Genome-Wide Mining of Disease Resistance Gene Analogs Using Conserved Domains

The production of legume crop species is severely affected by disease, imposing a significant yield loss annually worldwide. Plant resistance gene analogs (RGAs) play specific roles in plant resistance responses, and their identification and subsequent application in breeding programs help to reduce this yield loss. RGAs contain conserved domains and motifs, which can be used for their identification and classification. Nucleotide-binding site-leucine-rich repeat (NLR), receptor like kinase (RLK), and receptor like protein (RLP) genes are the main types of RGAs. Computational identification and characterization of RGAs has been performed successfully among different plant species. Here, we explain the computational workflow for genome-wide RGA identification in legumes.

FoldX accurate structural protein-DNA binding prediction using PADA1 (Protein Assisted DNA Assembly 1)

The speed at which new genomes are being sequenced highlights the need for genome-wide methods capable of predicting protein–DNA interactions. Here, we present PADA1, a generic algorithm that accurately models structural complexes and predicts the DNA-binding regions of resolved protein structures. PADA1 relies on a library of protein and double-stranded DNA fragment pairs obtained from a training set of 2103 DNA–protein complexes. It includes a fast statistical force field computed from atom-atom distances, to evaluate and filter the 3D docking models. Using published benchmark validation sets and 212 DNA–protein structures published after 2016 we predicted the DNA-binding regions with an RMSD of <1.8 Å per residue in >95% of the cases. We show that the quality of the docked templates is compatible with FoldX protein design tool suite to identify the crystallized DNA molecule sequence as the most energetically favorable in 80% of the cases. We highlighted the biological potential of PADA1 by reconstituting DNA and protein conformational changes upon protein mutagenesis of a meganuclease and its variants, and by predicting DNA-binding regions and nucleotide sequences in proteins crystallized without DNA. These results opens up new perspectives for the engineering of DNA–protein interfaces.

Characterization of a highly conserved Antheraea pernyi spermidine synthase gene

In the present study, we isolated a spermidine synthase gene from Antheraea pernyi (ApSpds) using expressed sequence tag method. The obtained cDNA sequence of 1483 bp contains an open-reading frame of 864 bp encoding a polypeptide of 287 amino acids. Sequence analysis revealed that ApSpds belonged to class I of AdoMet-MTase family, and exhibited 30% identity to those from bacteria, 45-48% identity to fungi, 36-47% identity to plants, 52-54% identity to vertebrates and 53-80% identity to invertebrates. Phylogenetic analysis found that the used Spds protein sequences were well divided into five groups corresponding to bacteria, fungi, plants, invertebrates and vertebrates, respectively. These results further confirmed that Spds is highly conserved through evolution of life organisms. The ApSpds mRNA is expressed during all four developmental stages and is present in all examined tissues with the highest abundance in the muscle, in which the relative mRNA expression level was 1.6 times higher than in the fat body. Although not significant, the mRNA level decreased after high-temperature exposure suggesting that the Spds gene may not be involved in temperature stress tolerance in A. pernyi. Taken together, our results suggested that ApSpds play an important role in development of silkworm.

Combining bioinformatics, cheminformatics, functional genomics and whole organism approaches for identifying epigenetic drug targets in Schistosoma mansoni

Schistosomiasis endangers the lives of greater than 200 million people every year and is predominantly controlled by a single class chemotherapy, praziquantel (PZQ). Development of PZQ replacement (to combat the threat of PZQ insensitivity/resistance arising) or combinatorial (to facilitate the killing of PZQ-insensitive juvenile schistosomes) chemotherapies would help sustain this control strategy into the future. Here, we re-categorise two families of druggable epigenetic targets in Schistosoma mansoni, the histone methyltransferases (HMTs) and the histone demethylases (HDMs). Amongst these, a S. mansoni Lysine Specific Demethylase 1 (SmLSD1, Smp_150560) homolog was selected for further analyses. Homology modelling of SmLSD1 and in silico docking of greater than four thousand putative inhibitors identified seven (L1 – L7) showing more favourable binding to the target pocket of SmLSD1 vs Homo sapiens HsLSD1; six of these seven (L1 – L6) plus three structural analogues of L7 (L8 – L10) were subsequently screened against schistosomula using the Roboworm anthelmintic discovery platform. The most active compounds (L10 - pirarubicin > L8 – danunorubicin hydrochloride) were subsequently tested against juvenile (3 wk old) and mature (7 wk old) schistosome stages and found to impede motility, suppress egg production and affect tegumental surfaces. When compared to a surrogate human cell line (HepG2), a moderate window of selectivity was observed for the most active compound L10 (selectivity indices - 11 for schistosomula, 9 for juveniles, 1.5 for adults). Finally, RNA interference of SmLSD1 recapitulated the egg-laying defect of schistosomes co-cultivated in the presence of L10 and L8. These preliminary results suggest that SmLSD1 represents an attractive new target for schistosomiasis; identification of more potent and selective SmLSD1 compounds, however, is essential. Nevertheless, the approaches described herein highlight an interdisciplinary strategy for selecting and screening novel/repositioned anti-schistosomals, which can be applied to any druggable (epigenetic) target encoded by the parasite's genome.

•

Schistosoma mansoni contains 27 histone methyltransferases (HMTs) and 14 histone demethylases (HDMs).

•

S. mansoni lysine specific demethylase 1 (SmLSD1) is a druggable target.

•

Schistosomes treated with the putative SmLSD1 inhibitor pirarubicin or siRNAs targeting SmLSD1 are less fecund.

BrEPS 2.0: Optimization of sequence pattern prediction for enzyme annotation

The prediction of gene functions is crucial for a large number of different life science areas. Faster high throughput sequencing techniques generate more and larger datasets. The manual annotation by classical wet-lab experiments is not suitable for these large amounts of data. We showed earlier that the automatic sequence pattern-based BrEPS protocol, based on manually curated sequences, can be used for the prediction of enzymatic functions of genes. The growing sequence databases provide the opportunity for more reliable patterns, but are also a challenge for the implementation of automatic protocols. We reimplemented and optimized the BrEPS pattern generation to be applicable for larger datasets in an acceptable timescale. Primary improvement of the new BrEPS protocol is the enhanced data selection step. Manually curated annotations from Swiss-Prot are used as reliable source for function prediction of enzymes observed on protein level. The pool of sequences is extended by highly similar sequences from TrEMBL and SwissProt. This allows us to restrict the selection of Swiss-Prot entries, without losing the diversity of sequences needed to generate significant patterns. Additionally, a supporting pattern type was introduced by extending the patterns at semi-conserved positions with highly similar amino acids. Extended patterns have an increased complexity, increasing the chance to match more sequences, without losing the essential structural information of the pattern. To enhance the usability of the database, we introduced enzyme function prediction based on consensus EC numbers and IUBMB enzyme nomenclature. BrEPS is part of the Braunschweig Enzyme Database (BRENDA) and is available on a completely redesigned website and as download. The database can be downloaded and used with the BrEPScmd command line tool for large scale sequence analysis. The BrEPS website and downloads for the database creation tool, command line tool and database are freely accessible at http://breps.tu-bs.de.

Transcriptomic profiling and its implications for the H2 production of a non-methanogen deficient in the frhAGB-encoding hydrogenase

The F₄₂₀-reducing hydrogenase of methanogens functions in methanogenesis by providing reduced coenzyme F₄₂₀ (F₄₂₀H₂) as an electron donor. In non-methanogens, however, their physiological function has not been identified yet. In this study, we constructed an ΔfrhA mutant, whose frhA gene encoding the hydrogenase α subunit was deleted, in the non-methanogenic Thermococcus onnurineus NA1 as a model organism. There was no significant difference in the formate-dependent growth between the mutant and the wild-type strains. Interestingly, the mutation in the frhA gene affected the expression of genes involved in various cellular functions such as H₂ oxidation, chemotactic signal transduction, and carbon monoxide (CO) metabolism. Among these genes, the CO oxidation gene cluster, enabling CO-dependent growth and H₂ production, showed a 2.8- to 7.0-fold upregulation by microarray-based whole transcriptome expression profiling. The levels of proteins produced by this gene cluster were also significantly increased not only under the formate condition but also under the CO condition. In a controlled bioreactor, where 100% CO was continuously fed, the ΔfrhA mutant exhibited significant increases in cell growth (2.8-fold) and H₂ production (3.4-fold). These findings strongly imply that this hydrogenase is functional in non-methanogens and is related to various cellular metabolic processes through an unidentified mechanism. An understanding of the mechanism by which the frhA gene deletion affected the expression of other genes will provide insights that can be applied to the development of strategies for the enhancement of H₂ production using CO as a substrate.

A novel Meloidogyne enterolobii effector MeTCTP promotes parasitism by suppressing programmed cell death in host plants

Meloidogyne enterolobii is one of the most important plant-parasitic nematodes that can overcome the Mi-1 resistance gene and damage many economically important crops. Translationally controlled tumour protein (TCTP) is a multifunctional protein that exists in various eukaryotes and plays an important role in parasitism. In this study, a novel M. enterolobii TCTP effector, named MeTCTP, was identified and functionally characterized. MeTCTP was specifically expressed within the dorsal gland and was up-regulated during M. enterolobii parasitism. Transient expression of MeTCTP in protoplasts from tomato roots showed that MeTCTP was localized in the cytoplasm of the host cells. Transgenic Arabidopsis thaliana plants overexpressing MeTCTP were more susceptible to M. enterolobii infection than wild-type plants in a dose-dependent manner. By contrast, in planta RNA interference (RNAi) targeting MeTCTP suppressed the expression of MeTCTP in infecting nematodes and attenuated their parasitism. Furthermore, MeTCTP could suppress programmed cell death triggered by the pro-apoptotic protein BAX. These results demonstrate that MeTCTP is a novel plant-parasitic nematode effector that promotes parasitism, probably by suppressing programmed cell death in host plants.

Two fibrinogen-like proteins, FGL1 and FGL2 are disulfide-linked subunits of oligomers that specifically bind nonviable spermatozoa

Nevertheless, a nonviable sperm population is present in the cauda epididymidis of many species. Degenerating spermatozoa release enzymes that could have detrimental effects on the viability of neighboring cells, and they are source of autoantigens that induce an autoimmune response if they escape the blood-epididymis barrier. Does the epididymis have specialized protective mechanism(s) to segregate the viable sperm population from defective spermatozoa? Previously, we identified a fibrinogen-like protein-2 (fgl2) that specifically binds to and polymerizes into a cocoon-like complex coating defective spermatozoa and sperm fragments. The objective of the present study is to identify the subunit composition of the fgl2-containing oligomers both in the soluble and cocoon-like complex. Our proteomic studies indicate that the 260/280kDa oligomers (termed eFGL) contain two distinct disulfide-linked subunits; 64kDa fgl2 and 33kDa fgl1. Utilizing a PCR-based cloning strategy, the 33kDa polypeptide has been identified as fibrinogen-like protein-1 (fgl1). Immunocytochemical studies revealed that fgl1 selectively binds to defective spermatozoa in the cauda epididymidis. Northern blot analysis and in situ hybridization demonstrated the high expression of fgl1 in the principal cells of the proximal cauda epididymidis. Co-immunoprecipitation analyses of cauda epididymal fluid, using anti-fgl2, demonstrate that both fgl1 and fgl2 are present in the soluble eFGL. Our study is the first to show an association of fgl1 and fgl2 both in the soluble and in the sperm-associated eFGL. We conclude that our results provide new insights into the mechanisms by which the potentially unique epididymal protein functions in the recognition and elimination of defective spermatozoa.

Exhaustive identification of human class II basic helix-loop-helix proteins by virtual library screening

Cellular proliferation, specification and differentiation in developing tissues are tightly coordinated by groups of transcription factors in response to extrinsic and intrinsic signals. Furthermore, renewable pools of stem cells in adult tissues are subject to similar regulation. Basic helix-loop-helix (bHLH) proteins are a group of transcription factors that exert such a determinative influence on a variety of developmental pathways from C. elegans to humans, and we wished to exclusively identify novel members from within the whole human bHLH family. We have, therefore, developed an 'empirical custom fingerprint', to define the class II bHLH domain and exclusively identify these proteins in silico. We have identified nine previously uncharacterised human class II proteins, four of which were novel, by interrogating conceptual translations of the GenBank HTGS database. RT-PCR and mammalian 2-hybrid analysis of a subset of the factors demonstrated that they were indeed expressed, and were able to interact with an appropriate binding partner in vitro. Thus, we are now approaching an almost complete listing of human class II bHLH factors.

Mammalian Target of Rapamycin Mediates Kidney Injury Molecule 1-Dependent Tubule Injury in a Surrogate Model

Kidney injury molecule 1 (KIM-1), an epithelial phagocytic receptor, is markedly upregulated in the proximal tubule in various forms of acute and chronic kidney injury in humans and many other species. Whereas acute expression of KIM-1 has adaptive anti-inflammatory effects, chronic expression may be maladaptive in mice. Here, we characterized the zebrafish Kim family, consisting of Kim-1, Kim-3, and Kim-4. Kim-1 was markedly upregulated in kidney after gentamicin-induced injury and had conserved phagocytic activity in zebrafish. Both constitutive and tamoxifen-induced expression of Kim-1 in zebrafish kidney tubules resulted in loss of the tubule brush border, reduced GFR, pericardial edema, and increased mortality. Kim-1-induced kidney injury was associated with reduction of growth of adult fish. Kim-1 expression led to activation of the mammalian target of rapamycin (mTOR) pathway, and inhibition of this pathway with rapamycin increased survival. mTOR pathway inhibition in KIM-1-overexpressing transgenic mice also significantly ameliorated serum creatinine level, proteinuria, tubular injury, and kidney inflammation. In conclusion, persistent Kim-1 expression results in chronic kidney damage in zebrafish through a mechanism involving mTOR. This observation predicted the role of the mTOR pathway and the therapeutic efficacy of mTOR-targeted agents in KIM-1-mediated kidney injury and fibrosis in mice, demonstrating the utility of the Kim-1 renal tubule zebrafish models.

Transcriptome sequencing and annotation of the polychaete Hermodice carunculata (Annelida, Amphinomidae)

Background

The amphinomid polychaete Hermodice carunculata is a cosmopolitan and ecologically important omnivore in coral reef ecosystems, preying on a diverse suite of reef organisms and potentially acting as a vector for coral disease. While amphinomids are a key group for determining the root of the Annelida, their phylogenetic position has been difficult to resolve, and their publically available genomic data was scarce.

Results

We performed deep transcriptome sequencing (Illumina HiSeq) and profiling on Hermodice carunculata collected in the Western Atlantic Ocean. We focused this study on 58,454 predicted Open Reading Frames (ORFs) of genes longer than 200 amino acids for our homology search, and Gene Ontology (GO) terms and InterPro IDs were assigned to 32,500 of these ORFs. We used this de novo assembled transcriptome to recover major signaling pathways and housekeeping genes. We also identify a suite of H. carunculata genes related to reproduction and immune response.

Conclusions

We provide a comprehensive catalogue of annotated genes for Hermodice carunculata and expand the knowledge of reproduction and immune response genes in annelids, in general. Overall, this study vastly expands the available genomic data for H. carunculata, of which previously consisted of only 279 nucleotide sequences in NCBI. This underscores the utility of Illumina sequencing for de novo transcriptome assembly in non-model organisms as a cost-effective and efficient tool for gene discovery and downstream applications, such as phylogenetic analysis and gene expression profiling.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1565-6) contains supplementary material, which is available to authorized users.

Prediction of multi-drug resistance transporters using a novel sequence analysis method

There are many examples of groups of proteins that have similar function, but the determinants of functional specificity may be hidden by lack of sequence similarity, or by large groups of similar sequences with different functions. Transporters are one such protein group in that the general function, transport, can be easily inferred from the sequence, but the substrate specificity can be impossible to predict from sequence with current methods. In this paper we describe a linguistic-based approach to identify functional patterns from groups of unaligned protein sequences and its application to predict multi-drug resistance transporters (MDRs) from bacteria. We first show that our method can recreate known patterns from PROSITE for several motifs from unaligned sequences. We then show that the method, MDRpred, can predict MDRs with greater accuracy and positive predictive value than a collection of currently available family-based models from the Pfam database. Finally, we apply MDRpred to a large collection of protein sequences from an environmental microbiome study to make novel predictions about drug resistance in a potential environmental reservoir.

Prediction of multi-drug resistance transporters using a novel sequence analysis method

There are many examples of groups of proteins that have similar function, but the determinants of functional specificity may be hidden by lack of sequence similarity, or by large groups of similar sequences with different functions. Transporters are one such protein group in that the general function, transport, can be easily inferred from the sequence, but the substrate specificity can be impossible to predict from sequence with current methods. In this paper we describe a linguistic-based approach to identify functional patterns from groups of unaligned protein sequences and its application to predict multi-drug resistance transporters (MDRs) from bacteria. We first show that our method can recreate known patterns from PROSITE for several motifs from unaligned sequences. We then show that the method, MDRpred, can predict MDRs with greater accuracy and positive predictive value than a collection of currently available family-based models from the Pfam database. Finally, we apply MDRpred to a large collection of protein sequences from an environmental microbiome study to make novel predictions about drug resistance in a potential environmental reservoir.

Impact of next-generation technologies on exploring socioeconomically important parasites and developing new interventions

High-throughput molecular and computer technologies have become instrumental for systems biological explorations of pathogens, including parasites. For instance, investigations of the transcriptomes of different developmental stages of parasitic nematodes give insights into gene expression, regulation and function in a parasite, which is a significant step to understanding their biology, as well as interactions with their host(s) and disease. This chapter (1) gives a background on some key parasitic nematodes of socioeconomic importance, (2) describes sequencing and bioinformatic technologies for large-scale studies of the transcriptomes and genomes of these parasites, (3) provides some recent examples of applications and (4) emphasizes the prospects of fundamental biological explorations of parasites using these technologies for the development of new interventions to combat parasitic diseases.

Structural and functional analysis of various globulin proteins from soy seed

Storage proteins of soybean mostly consist of globulins, which are classified according to their sedimentation coefficient. Among 4 major types: 2S, 7S, 11S, and 15S of globulins, 7S and 11S constitute major fraction. The 11S fraction consists only of glycinin and 7S fraction majorly consists of β-conglycinin, small amounts of γ-conglycinin and basic 7S globulin (Bg7S). Glycinin exist as a hexamer while β-conglycinin as a trimer and Bg7S as a tetramer. Glycinin subunits are coded by 5 genes of a family, whereas about 15 genes are present for β-conglycinin subunits. Bg7S gene is present in four copies in soybean genome. Synthesis of all proteins takes place as a single polypeptide chain, which is cleaved after folding to yield different chains or subunits. Glycinin and β-Conglycinin are made for storage purpose. However, Bg7S has potential xylanase inhibition activity and protein kinase activity. Primary structure of Bg7S reveals 12 conserved cysteine residues involved in forming 6 disulfide bonds, which provides appreciable stability to protein. Secondary structure is predominately rich in β-sheets with few alpha helices. Bg7S shares structural similarity with various aspartic-proteases. In this review, our aim is to discuss sequence, structure, and function of various globulins present in Glycine max.

Overview of software options for processing, analysis and interpretation of mass spectrometric proteomic data

Recently, the interests in proteomics have been intensively increased, and the proteomic methods have been widely applied to many problems in cell biology. If the age of 1990s is considered to be a decade of genomics, we can claim that the following years of the new century is a decade of proteomics. The rapid evolution of proteomics has continued through these years, with a series of innovations in separation techniques and the core technologies of two-dimensional gel electrophoresis and MS. Both technologies are fueled by automation and high throughput computation for profiling of proteins from biological systems. As Patterson ever mentioned, 'data analysis is the Achilles heel of proteomics and our ability to generate data now outstrips our ability to analyze it'. The development of automatic and high throughput technologies for rapid identification of proteins is essential for large-scale proteome projects and automatic protein identification and characterization is essential for high throughput proteomics. This review provides a snap shot of the tools and applications that are available for mass spectrometric high throughput biocomputation. The review starts with a brief introduction of proteomics and MS. Computational tools that can be employed at various stages of analysis are presented, including that for data processing, identification, quantification, and the understanding of the biological functions of individual proteins and their dynamic interactions. The challenges of computation software development and its future trends in MS-based proteomics have also been speculated.

Effectiveness and durability of the rice pi-ta gene in Yunnan province of China

Rice blast is one of the most damaging diseases of rice worldwide. In the present study, we analyzed DNA sequence variation of avirulence (AVR) genes of AVR-Pita1 in field isolates of Magnaporthe oryzae in order to understand the effectiveness of the resistance gene Pi-ta in China. Genomic DNA of 366 isolates of M. oryzae collected from Yunnan province of China were used for polymerase chain reaction (PCR) amplification to examine the existence of AVR-Pita1 using gene-specific PCR markers. Results of PCR products revealed that 218 isolates of M. oryzae carry AVR-Pita1. Among of them, 62.5, 56.3, 58.5, 46.7, 72.4, and 57.4% of M. oryzae carry AVR-Pita1 from northeastern, southeast, western, northwest, southwestern, and central Yunnan province, respectively. The detection rate of AVR-Pita1 was, in order: southwestern > northeastern > western > central > southeastern > northwestern Yunnan province. Moreover, in total, 18 AVR-Pita1 haplotypes encoding 13 novel AVR-Pita1 variants were identified among 60 isolates. Most DNA sequence variation was found to occur in the exon region, resulting in amino acid substitution. Six virulent haplotypes of AVR-Pita1 to Pita were identified among 60 field isolates. The AVR-Pita1 has evolved to virulence from avirulent origins via base substitution. These findings demonstrate that AVR-Pita1 is under positive selection and mutations of AVR-Pita1 are responsible for defeating race-specific resistance in nature.

Phytophthora viruses

Phytophthora sp. is a genus in the oomycetes, which are similar to filamentous fungi in morphology and habitat, but phylogenetically more closely related to brown algae and diatoms and fall in the kingdom Stramenopila. In the past few years, several viruses have been characterized in Phytophthora species, including four viruses from Phytophthora infestans, the late blight pathogen, and an endornavirus from an unnamed Phytophthora species from Douglas fir. Studies on Phytophthora viruses have revealed several interesting systems. Phytophthora infestans RNA virus 1 (PiRV-1) and PiRV-2 are likely the first members of two new virus families; studies on PiRV-3 support the establishment of a new virus genus that is not affiliated with established virus families; PiRV-4 is a member of Narnaviridae, most likely in the genus Narnavirus; and Phytophthora endornavirus 1 (PEV1) was the first nonplant endornavirus at the time of reporting. Viral capsids have not been found in any of the above-mentioned viruses. PiRV-1 demonstrated a unique genome organization that requires further examination, and PiRV-2 may have played a role in late blight resurgence in 1980s-1990s.

Protein databases on the internet

Protein databases have become a crucial part of modern biology. Huge amounts of data for protein structures, functions, and particularly sequences are being generated. Searching databases is often the first step in the study of a new protein. Comparison between proteins or between protein families provides information about the relationship between proteins within a genome or across different species, and hence offers much more information than can be obtained by studying only an isolated protein. In addition, secondary databases derived from experimental databases are also widely available. These databases reorganize and annotate the data or provide predictions. The use of multiple databases often helps researchers understand the structure and function of a protein. Although some protein databases are widely known, they are far from being fully utilized in the protein science community. This unit provides a starting point for readers to explore the potential of protein databases on the Internet.

Genome-wide computational approach for the prediction of duplications generating protein localization signals

Investigating the possible generation of motifs accountable for aberrant protein dislocation subsequent to the rise of short tandem duplications is interesting, given the pathogenic potential of this mechanism, as demonstrated in diseases such adult myeloid leukemia (AML). In this paper we introduce a new computational method for predicting genomic points which, after hypothetical mutation events such as micro-duplications, might encode molecular patterns such as localization or export signals. The proposed framework allows to study motifs of unconstrained length defined as regular expressions at a genome-wide level, providing an in silico platform capable of analyzing the potential effect of duplications on abnormal cellular localization.

Molecular characterization of a Cyrtochilum loxense Somatic Embryogenesis Receptor-like Kinase (SERK) gene expressed during somatic embryogenesis

Somatic embryogenesis is crucial for the propagation of endangered Ecuadorian orchid species, among them Cyrtochilum loxense, in view of the fact that their number in nature or in collections is quite reduced. One of the genes expressed during somatic and zygotic embryogenesis is Somatic Embryogenesis Receptor-like Kinase (SERK). Despite the development of somatic embryogenesis protocols for orchids, no SERK genes have been isolated from this family. This is the first report on the isolation of a full-length orchid SERK sequence, namely that of Cyrtochilum loxense (ClSERK). The identity of ClSERK was inferred by the presence of all domains typical of SERK proteins: a signal peptide, a leucine zipper domain, five LRRs, a serine proline-rich domain, a transmembrane domain, a kinase domain, and the C-terminal region. We have observed that the ClSERK gene is highly expressed in embryogenic calluses generated from protocorms at the time of appearance of embryonic morphological features. At later stages when embryos become well visible on calluses, ClSERK gene expression decreases. Compared to early stages of embryo formation on calluses, the expression detected in leaf tissue is far lower, thus suggesting a role of this gene during development.

ALTERED ABUNDANCE OF MESSENGER RNA TRANSCRIPTS WHOSE PRODUCTS HAVE MEMBRANE ATTACHMENT SITES IN PREGNANCY-INDUCED HYPERTENSION

Pregnancy-induced hypertensive disorders (PIH) are leading causes of maternal mortality. Although the mechanism responsible for initiating and maintaining the disorder is unproven, physiologic molecular attachments in kidney and placenta play a role. The SHHF/Mcc-facp (SHHF) rat has features of the disorder, including abnormal placenta gene expression. To gain a molecular understanding of the gene expression profile associated with PIH, kidneys and placentas of SHHF rats at gestation day 20 were compared to WKY controls using microarray technology. We report that SHHF rats have spontaneous PIH, elevated total placenta weights, and reduced total pup weights than WKY controls and that they also have greater total number of mRNA transcripts expressed in placenta. Kidneys of SHHF rats, on the other hand, not only expressed disproportionately more predicted gene products with attachment sites such as RGD motifs, N-glycosylation sites, and N-myristoylation sites they also responded more profoundly to oral administration of L-arginine. We conclude that the increased abundance of transcripts whose products engage in posttranslational attachments using RGD motifs, N-glycosylation sites, and N-myristoylation sites and the reversal of these increases by oral administration of L-arginine suggests that NO may be of importance in PIH at the level of molecular attachments.

Acidic α-galactosidase is the most abundant nectarin in floral nectar of common tobacco (Nicotiana tabacum)

BACKGROUND AND AIMS

To date, most floral nectarins (nectar proteins) are reported to function in nectar defence, particularly for insect-pollinated outcrossing species. We compared nectarin composition and abundance in selfing common tobacco (Nicotiana tobaccum) with outcrossing ornamental tobacco plants to elucidate the functional difference of nectarins in different reproductive systems.

METHODS

Common tobacco (CT) nectarins were separated by SDS-PAGE and the N terminus of the most abundant nectarin was sequenced via Edman degradation. The full-length nectarin gene was amplified and cloned from genomic DNA and mRNA with hiTail-PCR and RACE (rapid amplification of cDNA ends), and expression patterns were then investigated in different tissues using semi-quantitative reverse transcriptase PCR. Additionally, high-performance liquid chromatography and enzymatic analyses of nectar sugar composition, and other biochemical traits and functions of the novel nectarin were studied.

KEY RESULTS

The most abundant nectarin in CT nectar is an acidic α-galactosidase, here designated NTα-Gal. This compound has a molecular mass of 40 013 Da and a theoretical pI of 5·33. NTα-Gal has a conserved α-Gal characteristic signature, encodes a mature protein of 364 amino acids and is expressed in different organs. Compared with 27 other melliferous plant species from different families, CT floral nectar demonstrated the highest α-Gal activity, which is inhibited by d-galactose. Raffinose family oligosaccharides were not detected in CT nectar, indicating that NTα-Gal does not function in post-secretory hydrolysis. Moreover, tobacco plant fruits did not develop intact skin with galactose inhibition of NTα-Gal activity in nectar, suggesting that NTα-Gal induces cell-wall surface restructuring during the initial stages of fruit development.

CONCLUSIONS

α-Gal was the most abundant nectarin in selfing CT plants, but was not detected in the nectar of strictly outcrossing sister tobacco species. No function was demonstrated in antimicrobial defence. Therefore, floral nectarins in selfing species maintain their functional significance in reproductive organ development.

Heartworm genomics: unprecedented opportunities for fundamental molecular insights and new intervention strategies

Vector-borne diseases, including canine heartworm disease (CHWD), are of major socioeconomic and canine health importance worldwide. Although many studies have provided insights into CHWD, to date there has been limited study of fundamental molecular aspects of Dirofilaria immitis itself, its relationship with the canine host, its vectors, as well as the potential of drug resistance to emerge, using advanced -omic technologies. This article takes a prospective view of the benefits that advanced -omics technologies will have toward understanding D. immitis and CHWD. Tackling key biological questions using these technologies will provide a "systems biology" context and could lead to radically new intervention and management strategies against heartworm.

Lipases: an overview

Lipases are ubiquitous enzymes, widespread in nature. They were first isolated from bacteria in the early nineteenth century and the associated research continuously increased due to the particular characteristics of these enzymes. This chapter reviews the main sources, structural properties, and industrial applications of these highly studied enzymes.

Co-expressed peJK genes of lobster (Jasus edwardsii)

Previous studies have shown that the two novel genes of southern rock lobster (Jasus edwardsii) named peJK2 and peJK3 are implicated in eyestalk hormonal regulation of the lobster moult cycle. Northern blot, in situ hybridisation studies and sequence analyses showed that their putative products might be transmembrane proteins associated with cell signal transduction of hormonal signals in the eyestalk during the intermoult phase of the moult cycle. The aim of this study was to analyse coexpression of peJK genes in different J. edwardsii tissues. Using reverse transcriptase–polymerase chain reaction (RT-PCR), the expression of peJK genes was analysed in seven different tissues (eyestalk, brain, epidermis, hepatopancreas, gill, muscle and heart) of an intermoult lobster. During RT-PCR analysis, a novel sequence was isolated, and was named peJK4. It shares 88% and 86% sequence identity with peJK2 and peJK3 respectively. The peJK2 and peJK4 genes are expressed in all tested tissues. Sequence analyses of the predicted peJK2 and peJK4 proteins revealed two common signal transduction motifs, transmembrane helices and protein kinase C. These results showed that the peJK genes of J. edwardsii are a complex group of genes and possibly involved in different signal transduction pathways.

Crystallization and preliminary crystallographic analysis of a putative glucokinase/hexokinase from Thermus thermophilus

Glucokinase/hexokinase catalyzes the phosphorylation of glucose to glucose 6-phosphate, which is the first step of glycolysis. The open reading frame TTHA0299 of the extreme thermophile Thermus thermophilus encodes a putative glucokinase/hexokinase which contains the consensus sequence for proteins from the repressors, open reading frames and sugar kinases family. In this study, the glucokinase/hexokinase from T. thermophilus was purified and crystallized using polyethylene glycol 8000 as a precipitant. Diffraction data were collected and processed to 2.02 Å resolution. The crystal belonged to space group P2(1), with unit-cell parameters a = 70.93, b = 138.14, c = 75.16 Å, β = 95.41°.

RULEMINER: A KNOWLEDGE SYSTEM FOR SUPPORTING HIGH-THROUGHPUT PROTEIN FUNCTION ANNOTATIONS

In this paper, we present RuleMiner, a knowledge system to facilitate a seamless integration of multi-sequence analysis tools and define profile-based rules for supporting high-throughput protein function annotations. This system consists of three essential components, Protein Function Groups (PFGs), PFG profiles and rules. The PFGs, established from an integrated analysis of current knowledge of protein functions from Swiss-Prot database and protein family-based sequence classifications, cover all possible cellular functions available in the database. The PFG profiles illustrate detailed protein features in the PFGs as in sequence conservations, the occurrences of sequence-based motifs, domains and species distributions. The rules, extracted from the PFG profiles, describe the clear relationships between these PFGs and all possible features. As a result, the RuleMiner is able to provide an enhanced capability for protein function analysis, such as results from the integrated sequence analysis tools for given proteins can be comparatively analyzed due to the clear feature-PFG relationships. Also, much needed guidance is readily available for such analysis. If the rules describe one-to-one (unique) relationships between the protein features and the PFGs, then these features can be utilized as unique functional identifiers and cellular functions of unknown proteins can be reliably determined. Otherwise, additional information has to be provided.

Tandem repeats in proteins: from sequence to structure

The bioinformatics analysis of proteins containing tandem repeats requires special computer programs and databases, since the conventional approaches predominantly developed for globular domains have limited success. Here, I survey bioinformatics tools which have been developed recently for identification and proteome-wide analysis of protein repeats. The last few years have also been marked by an emergence of new 3D structures of these proteins. Appraisal of the known structures and their classification uncovers a straightforward relationship between their architecture and the length of the repetitive units. This relationship and the repetitive character of structural folds suggest rules for better prediction of the 3D structures of such proteins. Furthermore, bioinformatics approaches combined with low resolution structural data, from biophysical techniques, especially, the recently emerged cryo-electron microscopy, lead to reliable prediction of the protein repeat structures and their mode of binding with partners within molecular complexes. This hybrid approach can actively be used for structural and functional annotations of proteomes.

Distinct properties of human HMGN5 reveal a rapidly evolving but functionally conserved nucleosome binding protein

The HMGN family is a family of nucleosome-binding architectural proteins that affect the structure and function of chromatin in vertebrates. We report that the HMGN5 variant, encoded by a gene located on chromosome X, is a rapidly evolving protein with an acidic C-terminal domain that differs among vertebrate species. We found that the intranuclear organization and nucleosome interactions of human HMGN5 are distinct from those of mouse HMGN5 and that the C-terminal region of the protein is the main determinant of the chromatin interaction properties. Despite their apparent differences, both mouse and human HMGN5 proteins interact with histone H1, reduce its chromatin residence time, and can induce large-scale chromatin decompaction in living cells. Analysis of HMGN5 mutants suggests that distinct domains in HMGN5 affect specific steps in the interaction of H1 with chromatin. Elevated levels of either human or mouse HMGN5 affect the transcription of numerous genes, most in a variant-specific manner. Our study identifies HMGN5 as a rapidly evolving vertebrate nuclear protein with species-specific properties. HMGN5 has a highly disordered structure, binds dynamically to nucleosome core particles, modulates the binding of H1 to chromatin, reduces the compaction of the chromatin fiber, and affects transcription.

Major prospects for exploring canine vector borne diseases and novel intervention methods using 'omic technologies

Canine vector-borne diseases (CVBDs) are of major socioeconomic importance worldwide. Although many studies have provided insights into CVBDs, there has been limited exploration of fundamental molecular aspects of most pathogens, their vectors, pathogen-host relationships and disease and drug resistance using advanced, 'omic technologies. The aim of the present article is to take a prospective view of the impact that next-generation, 'omics technologies could have, with an emphasis on describing the principles of transcriptomic/genomic sequencing as well as bioinformatic technologies and their implications in both fundamental and applied areas of CVBD research. Tackling key biological questions employing these technologies will provide a 'systems biology' context and could lead to radically new intervention and management strategies against CVBDs.

Molecular cloning and characterization of a calreticulin cDNA from the pinewood nematode Bursaphelenchus xylophilus

The cloning and characterization of a cDNA encoding a calreticulin from the pinewood nematode Bursaphelenchus xylophilus is described herein. The full-length cDNA (Bx-crt-1) contained a 1200 bp open reading frame that could be translated to a 399 amino acid polypeptide. The deduced protein contained highly conserved regions of a calreticulin gene and had 66.2-70.1% amino acid sequence identity to other calreticulin sequences from nematodes. RNAi, RT-PCR amplification, and southern blot suggest that Bx-crt-1 may be important for the development of B. xylophilus.