FSP4HSP: Frequent sequential patterns for the improved classification of heat shock proteins, their families, and sub-types.

Heat shock proteins (HSPs) from different families and sub-types play a vital role in the folding and unfolding of proteins, in maintaining cellular health, and in preventing serious disorders. Previous computational methods for HSP classification have yielded promising performance. However, most of the existing methods rely heavily on amino acid composition features and still face challenges related to interpretability and accuracy. To overcome these issues, we introduce a novel frequent sequential pattern (FSP)-based analysis and classification method for the classification of HSPs, their families, and sub-types. The proposed method is called FSP4HSP, which stands for "FSP for HSP". It identifies FSPs of amino acids (FSPAAs) and utilizes them for analysis and classification. Besides FSPAAs, sequential rules among amino acids are also discovered. Both binary and multi-class classification scenarios are considered, with the utilization of eight integer-based and four string-based classifiers. The incorporation of FSPAAs in the classification/prediction task enhances the interpretability of FSP4HSP and a comprehensive performance comparison using various evaluation measures demonstrates that it surpasses existing methods for the classification/recognition of HSPs.

Identifying Heat Shock Protein Families from Imbalanced Data by Using Combined Features.

Heat shock proteins (HSPs) are ubiquitous in living organisms. HSPs are an essential component for cell growth and survival; the main function of HSPs is controlling the folding and unfolding process of proteins. According to molecular function and mass, HSPs are categorized into six different families: HSP20 (small HSPS), HSP40 (J-proteins), HSP60, HSP70, HSP90, and HSP100. In this paper, improved methods for HSP prediction are proposed-the split amino acid composition (SAAC), the dipeptide composition (DC), the conjoint triad feature (CTF), and the pseudoaverage chemical shift (PseACS) were selected to predict the HSPs with a support vector machine (SVM). In order to overcome the imbalance data classification problems, the syntactic minority oversampling technique (SMOTE) was used to balance the dataset. The overall accuracy was 99.72% with a balanced dataset in the jackknife test by using the optimized combination feature SAAC+DC+CTF+PseACS, which was 4.81% higher than the imbalanced dataset with the same combination feature. The Sn, Sp, Acc, and MCC of HSP families in our predictive model were higher than those in existing methods. This improved method may be helpful for protein function prediction.

Genome-wide Identification and Characterization of Heat Shock Protein Family Reveals Role in Development and Stress Conditions in Triticum aestivum L.

Heat shock proteins (HSPs) have a significant role in protein folding and are considered as prominent candidates for development of heat-tolerant crops. Understanding of wheat HSPs has great importance since wheat is severely affected by heat stress, particularly during the grain filling stage. In the present study, efforts were made to identify HSPs in wheat and to understand their role during plant development and under different stress conditions. HSPs in wheat genome were first identified by using Position-Specific Scoring Matrix (PSSMs) of known HSP domains and then also confirmed by sequence homology with already known HSPs. Collectively, 753 TaHSPs including 169 TaSHSP, 273 TaHSP40, 95 TaHSP60, 114 TaHSP70, 18 TaHSP90 and 84 TaHSP100 were identified in the wheat genome. Compared with other grass species, number of HSPs in wheat was relatively high probably due to the higher ploidy level. Large number of tandem duplication was identified in TaHSPs, especially TaSHSPs. The TaHSP genes showed random distribution on chromosomes, however, there were more TaHSPs in B and D sub-genomes as compared to the A sub-genome. Extensive computational analysis was performed using the available genomic resources to understand gene structure, gene expression and phylogentic relationship of TaHSPs. Interestingly, apart from high expression under heat stress, high expression of TaSHSP was also observed during seed development. The study provided a list of candidate HSP genes for improving thermo tolerance during developmental stages and also for understanding the seed development process in bread wheat.

Recent Advances in Machine Learning Methods for Predicting Heat Shock Proteins.

BACKGROUND: As molecular chaperones, Heat Shock Proteins (HSPs) not only play key roles in protein folding and maintaining protein stabilities, but are also linked with multiple kinds of diseases. Therefore, HSPs have been regarded as the focus of drug design. Since HSPs from different families play distinct functions, accurately classifying the families of HSPs is the key step to clearly understand their biological functions. In contrast to laborintensive and cost-ineffective experimental methods, computational classification of HSP families has emerged to be an alternative approach. METHODS: We reviewed the paper that described the existing datasets of HSPs and the representative computational approaches developed for the identification and classification of HSPs. RESULTS: The two benchmark datasets of HSPs, namely HSPIR and sHSPdb were introduced, which provided invaluable resources for computationally identifying HSPs. The gold standard dataset and sequence encoding schemes for building computational methods of classifying HSPs were also introduced. The three representative web-servers for identifying HSPs and their families were described. CONCLUSION: The existing machine learning methods for identifying the different families of HSPs indeed yielded quite encouraging results and did play a role in promoting the research on HSPs. However, the number of HSPs with known structures is very limited. Therefore, determining the structure of the HSPs is also urgent, which will be helpful in revealing their functions.

Identification and Functional Characterization of a Cold-Related Protein, BcHHP5, in Pak-Choi (Brassica rapa ssp. chinensis).

In plants, heptahelical proteins (HHPs) have been shown to respond to a variety of abiotic stresses, including cold stress. Up to the present, the regulation mechanism of HHP5 under low temperature stress remains unclear. In this study, BcHHP5 was isolated from Pak-choi (Brassica rapa ssp. chinensis cv. Suzhouqing). Sequence analysis and phylogenetic analysis indicated that BcHHP5 in Pak-choi is similar to AtHHP5 in Arabidopsis thaliana. Structure analysis showed that the structure of the BcHHP5 protein is relatively stable and highly conservative. Subcellular localization indicated that BcHHP5 was localized on the cell membrane and nuclear membrane. Furthermore, real-time quantitative polymerase chain reaction (RT-qPCR) analysis showed that BcHHP5 was induced to express by cold and other abiotic stresses. In Pak-choi, BcHHP5-silenced assay, inhibiting the action of endogenous BcHHP5, indicated that BcHHP5-silenced might have a negative effect on cold tolerance, which was further confirmed. All of these results indicate that BcHHP5 might play a role in abiotic response. This work can serve as a reference for the functional analysis of other cold-related proteins from Pak-choi in the future.

Photosynthetic and Stress Responsive Proteins Are Altered More Effectively in Nicotiana benthamiana Infected with Plum pox virus Aggressive PPV-CR versus Mild PPV-C Cherry-Adapted Isolates.

Plum pox virus (PPV, family Potyviridae) is one of the most important viral pathogens of Prunus spp. causing considerable damage to stone-fruit industry worldwide. Among the PPV strains identified so far, only PPV-C, PPV-CR, and PPV-CV are able to infect cherries under natural conditions. Herein, we evaluated the pathogenic potential of two viral isolates in herbaceous host Nicotiana benthamiana. Significantly higher accumulation of PPV capsid protein in tobacco leaves infected with PPV-CR (RU-30sc isolate) was detected in contrast to PPV-C (BY-101 isolate). This result correlated well with the symptoms observed in the infected plants. To further explore the host response upon viral infection at the molecular level, a comprehensive proteomic profiling was performed. Using reverse-phase ultra-high-performance liquid chromatography followed by label-free mass spectrometry quantification, we identified 38 unique plant proteins as significantly altered due to the infection. Notably, the abundances of photosynthesis-related proteins, mainly from the Calvin-Benson cycle, were found more aggressively affected in plants infected with PPV-CR isolate than those of PPV-C. This observation was accompanied by a significant reduction in the amount of photosynthetic pigments extracted from the leaves of PPV-CR infected plants. Shifts in the abundance of proteins that are involved in stimulation of photosynthetic capacity, modification of amino acid, and carbohydrate metabolism may affect plant growth and initiate energy formation via gluconeogenesis in PPV infected N. benthamiana. Furthermore, we suggest that the higher accumulation of H2O2 in PPV-CR infected leaves plays a crucial role in plant defense and development by activating the glutathione synthesis.

Characterization of six small HSP genes from Chironomus riparius (Diptera, Chironomidae): Differential expression under conditions of normal growth and heat-induced stress.

Small heat shock proteins (sHSPs) comprise the most numerous, structurally diverse, and functionally uncharacterized family of heat shock proteins. Several Hsp genes (Hsp 90, 70, 40, and 27) from the insect Chironomus riparius are widely used in aquatic toxicology as biomarkers for environmental toxins. Here, we conducted a comparative study and characterized secondary structure of the six newly identified sHsp genes Hsp17, Hsp21, Hsp22, Hsp23, Hsp24, and Hsp34. A characteristic α-crystallin domain is predicted in all the new proteins. Phylogenetic analysis suggests a strong relation to other sHSPs from insects and interesting evidence regarding evolutionary origin and duplication events. Comparative analysis of transcription profiles for Hsp27, Hsp70, and the six newly identified genes revealed that Hsp17, Hsp21, and Hsp22 are constitutively expressed under normal conditions, while under two different heat shock conditions these genes are either not activated or are even repressed (Hsp22). In contrast, Hsp23, Hsp24, and Hsp34 are significantly activated along with Hsp27 and Hsp70 during heat stress. These results strongly suggest functional differentiation within the small HSP subfamily and provide new data to help understand the coping mechanisms induced by stressful environmental stimuli.

Functional diversification and specialization of cytosolic 70-kDa heat shock proteins.

A fundamental question in molecular evolution is how protein functional differentiation alters the ability of cells and organisms to cope with stress and survive. To answer this question we used two paralogous Hsp70s from mouse and explored whether these highly similar cytosolic molecular chaperones, which apart their temporal expression have been considered functionally interchangeable, are differentiated with respect to their lipid-binding function. We demonstrate that the two proteins bind to diverse lipids with different affinities and therefore are functionally specialized. The observed lipid-binding patterns may be related with the ability of both Hsp70s to induce cell death by binding to a particular plasma-membrane lipid, and the potential of only one of them to promote cell survival by binding to a specific lysosomal-membrane lipid. These observations reveal that two seemingly identical proteins differentially modulate cellular adaptation and survival by having acquired specialized functions via sequence divergence. Therefore, this study provides an evolutionary paradigm, where promiscuity, specificity, sub- and neo-functionalization orchestrate one of the most conserved systems in nature, the cellular stress-response.

Heat shock factors in carrot: genome-wide identification, classification, and expression profiles response to abiotic stress.

Heat shock factors (HSFs) play key roles in the response to abiotic stress in eukaryotes. In this study, 35 DcHSFs were identified from carrot (Daucus carota L.) based on the carrot genome database. All 35 DcHSFs were divided into three classes (A, B, and C) according to the structure and phylogenetic relationships of four different plants, namely, Arabidopsis thaliana, Vitis vinifera, Brassica rapa, and Oryza sativa. Comparative analysis of algae, gymnosperms, and angiosperms indicated that the numbers of HSF transcription factors were related to the plant's evolution. The expression profiles of five DcHsf genes (DcHsf 01, DcHsf 02, DcHsf 09, DcHsf 10, and DcHsf 16), which selected from each subfamily (A, B, and C), were detected by quantitative real-time PCR under abiotic stresses (cold, heat, high salinity, and drought) in two carrot cultivars, D. carota L. cvs. Kurodagosun and Junchuanhong. The expression levels of DcHsfs were markedly increased by heat stress, except that of DcHsf 10, which was down regulated. The expression profiles of different DcHsfs in the same class also differed under various stress treatments. The expression profiles of these DcHsfs were also different in tissues of two carrot cultivars. This study is the first to identify and characterize the DcHSF family transcription factors in plants of Apiaceae using whole-genome analysis. The results of this study provide an in-depth understanding of the DcHSF family transcription factors' structure, function, and evolution in carrot.

The heat shock factor family from Triticum aestivum in response to heat and other major abiotic stresses and their role in regulation of heat shock protein genes.

Heat shock factors (Hsfs) play a central regulatory role in acquired thermotolerance. To understand the role of the major molecular players in wheat adaptation to heat stress, the Hsf family was investigated in Triticum aestivum. Bioinformatic and phylogenetic analyses identified 56 TaHsf members, which are classified into A, B, and C classes. Many TaHsfs were constitutively expressed. Subclass A6 members were predominantly expressed in the endosperm under non-stress conditions. Upon heat stress, the transcript levels of A2 and A6 members became the dominant Hsfs, suggesting an important regulatory role during heat stress. Many TaHsfA members as well as B1, C1, and C2 members were also up-regulated during drought and salt stresses. The heat-induced expression profiles of many heat shock protein (Hsp) genes were paralleled by those of A2 and A6 members. Transactivation analysis revealed that in addition to TaHsfA members (A2b and A4e), overexpression of TaHsfC2a activated expression of TaHsp promoter-driven reporter genes under non-stress conditions, while TaHsfB1b and TaHsfC1b did not. Functional heat shock elements (HSEs) interacting with TaHsfA2b were identified in four TaHsp promoters. Promoter mutagenesis analysis demonstrated that an atypical HSE (GAACATTTTGGAA) in the TaHsp17 promoter is functional for heat-inducible expression and transactivation by Hsf proteins. The transactivation of Hsp promoter-driven reporter genes by TaHsfC2a also relied on the presence of HSE. An activation motif in the C-terminal domain of TaHsfC2a was identified by amino residue substitution analysis. These data demonstrate the role of HsfA and HsfC2 in regulation of Hsp genes in wheat.

Ecological adaption analysis of the cotton aphid (Aphis gossypii) in different phenotypes by transcriptome comparison.

BACKGROUND: The cotton aphid, Aphis gossypii Glover, is a destructive insect pest worldwide; it directly or indirectly damages (virus transmission) 300 species of host plants. Knowledge of their ecologically adaptive mechanisms at the molecular level may provide an essential and urgent method to effectively control this pest. However, no transcriptome information is available for the cotton aphid and sequence data are scarce. Therefore, we obtained transcriptome data. RESULTS: To facilitate such a study, two cotton aphid transcriptomes at different growth stages of cotton, seedling and summer, were sequenced. A total of 161,396 and 66,668 contigs were obtained and assembled into 83,671 and 42,438 transcripts, respectively. After combining the raw date for both transcriptomes, the sequences were reassembled into 66,695 transcripts, and 52,160 were annotated based on BLASTX analyses. Comparison of the transcriptomes revealed that summer presented less challenges for the cotton aphids than the seedling stage of cotton. In total, 58 putative heat shock protein genes and 66 candidate cytochrome p450 genes were identified with BLASTX. CONCLUSIONS: Our results form a basis for exploring the molecular mechanisms of ecological adaption in the cotton aphid. Our study also provides a baseline for the exploration of abiotic stress responses. In addition, it provides large-scale sequence information for further studies on this species.

De novo assembly, gene annotation, and marker discovery in stored-product pest Liposcelis entomophila (Enderlein) using transcriptome sequences.

BACKGROUND: As a major stored-product pest insect, Liposcelis entomophila has developed high levels of resistance to various insecticides in grain storage systems. However, the molecular mechanisms underlying resistance and environmental stress have not been characterized. To date, there is a lack of genomic information for this species. Therefore, studies aimed at profiling the L. entomophila transcriptome would provide a better understanding of the biological functions at the molecular levels. METHODOLOGY/PRINCIPAL FINDINGS: We applied Illumina sequencing technology to sequence the transcriptome of L. entomophila. A total of 54,406,328 clean reads were obtained and that de novo assembled into 54,220 unigenes, with an average length of 571 bp. Through a similarity search, 33,404 (61.61%) unigenes were matched to known proteins in the NCBI non-redundant (Nr) protein database. These unigenes were further functionally annotated with gene ontology (GO), cluster of orthologous groups of proteins (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of genes potentially involved in insecticide resistance were manually curated, including 68 putative cytochrome P450 genes, 37 putative glutathione S-transferase (GST) genes, 19 putative carboxyl/cholinesterase (CCE) genes, and other 126 transcripts to contain target site sequences or encoding detoxification genes representing eight types of resistance enzymes. Furthermore, to gain insight into the molecular basis of the L. entomophila toward thermal stresses, 25 heat shock protein (Hsp) genes were identified. In addition, 1,100 SSRs and 57,757 SNPs were detected and 231 pairs of SSR primes were designed for investigating the genetic diversity in future. CONCLUSIONS/SIGNIFICANCE: We developed a comprehensive transcriptomic database for L. entomophila. These sequences and putative molecular markers would further promote our understanding of the molecular mechanisms underlying insecticide resistance or environmental stress, and will facilitate studies on population genetics for psocids, as well as providing useful information for functional genomic research in the future.

Hsp21 potentiates antifungal drug tolerance in Candida albicans.

Systemic infections of humans with the fungal pathogen Candida albicans are associated with a high mortality rate. Currently, efficient treatment of these infections is hampered by the relatively low number of available antifungal drugs. We recently identified the small heat shock protein Hsp21 in C. albicans and demonstrated its fundamental role for environmental stress adaptation and fungal virulence. Hsp21 was found in several pathogenic Candida species but not in humans. This prompted us to investigate the effects of a broad range of different antifungal drugs on an Hsp21-null C. albicans mutant strain. Our results indicate that combinatorial therapy targeting Hsp21, together with specific antifungal drug targets, has strong synergistic potential. In addition, we demonstrate that Hsp21 is required for tolerance to ethanol-induced stress and induction of filamentation in response to pharmacological inhibition of Hsp90. These findings might pave the way for the development of new treatment strategies against Candida infections.

An unusual dimeric small heat shock protein provides insight into the mechanism of this class of chaperones.

Small heat shock proteins (sHSPs) are virtually ubiquitous stress proteins that are also found in many normal tissues and accumulate in diseases of protein folding. They generally act as ATP-independent chaperones to bind and stabilize denaturing proteins that can be later reactivated by ATP-dependent Hsp70/DnaK, but the mechanism of substrate capture by sHSPs remains poorly understood. A majority of sHSPs form large oligomers, a property that has been linked to their effective chaperone action. We describe AtHsp18.5 from Arabidopsis thaliana, demonstrating that it is dimeric and exhibits robust chaperone activity, which adds support to the model that suboligomeric sHSP forms are a substrate binding species. Notably, like oligomeric sHSPs, when bound to substrate, AtHsp18.5 assembles into large complexes, indicating that reformation of sHSP oligomeric contacts is not required for assembly of sHSP-substrate complexes. Monomers of AtHsp18.5 freely exchange between dimers but fail to coassemble in vitro with dodecameric plant cytosolic sHSPs, suggesting that AtHsp18.5 does not interact by coassembly with these other sHSPs in vivo. Data from controlled proteolysis and hydrogen-deuterium exchange coupled with mass spectrometry show that the N- and C-termini of AtHsp18.5 are highly accessible and lack stable secondary structure, most likely a requirement for substrate interaction. Chaperone activity of a series of AtHsp18.5 truncation mutants confirms that the N-terminal arm is required for substrate protection and that different substrates interact differently with the N-terminal arm. In total, these data imply that the core α-crystallin domain of the sHSPs is a platform for flexible arms that capture substrates to maintain their solubility.

The family of Deg/HtrA proteases in plants.

BACKGROUND: The Deg/HtrA family of ATP-independent serine endopeptidases is present in nearly all organisms from bacteria to human and vascular plants. In recent years, multiple deg/htrA protease genes were identified in various plant genomes. During genome annotations most proteases were named according to the order of discovery, hence the same names were sometimes given to different types of Deg/HtrA enzymes in different plant species. This can easily lead to false inference of individual protease functions based solely on a shared name. Therefore, the existing names and classification of these proteolytic enzymes does not meet our current needs and a phylogeny-based standardized nomenclature is required. RESULTS: Using phylogenetic and domain arrangement analysis, we improved the nomenclature of the Deg/HtrA protease family, standardized protease names based on their well-established nomenclature in Arabidopsis thaliana, and clarified the evolutionary relationship between orthologous enzymes from various photosynthetic organisms across several divergent systematic groups, including dicots, a monocot, a moss and a green alga. Furthermore, we identified a "core set" of eight proteases shared by all organisms examined here that might provide all the proteolytic potential of Deg/HtrA proteases necessary for a hypothetical plant cell. CONCLUSIONS: In our proposed nomenclature, the evolutionarily closest orthologs have the same protease name, simplifying scientific communication when comparing different plant species and allowing for more reliable inference of protease functions. Further, we proposed that the high number of Deg/HtrA proteases in plants is mainly due to gene duplications unique to the respective organism.

Evolutionary analysis of prokaryotic heat-shock transcription regulatory protein σ³².

Heat-stress to any living cell is known to trigger a universal defense response, called heat-shock response, with rapid induction of tens of different heat-shock proteins. Bacterial heat-shock genes are transcribed by the σ(32)-bound RNA polymerase instead of the normal σ(70)-bound RNA polymerase. In this study, the diversity in sequence, variation in secondary structure and function amongst the different functional regions of the proteobacterial σ(32) family of proteins, and their phylogenetic relationships have been analyzed. Bacterial σ(32) proteins can be subdivided into different functional regions which are referred to as regions 2, 3, and 4. There is a great deal of sequence conservation among the functional regions of proteobacterial σ(32) family of proteins though some mutations are also present in these regions. Region 2 is the most conserved one, while region 4 has comparatively more variable sequences. In the present work, we tried to explore the effects of mutations in these regions. Our study suggests that the sequence diversities due to natural mutations in the different regions of proteobacterial σ(32) family lead to different functions. So far, this study is the first bioinformatic approach towards the understanding of the mechanistic details of σ(32) family of proteins using the protein sequence information only. This study therefore may help in elucidating the hitherto unknown molecular mechanism of the functionalities of σ(32)family of proteins.

Gene expression analysis of heat shock proteins in the nucleus accumbens of rats with different morphine seeking behaviours.

Heat-shock proteins play functional roles on brain regulatory processes which are deeply involved in drug addiction, such as synaptic plasticity. However, few studies have been focused on gene expression of heat-shock proteins (Hsp) as potential diagnostic tools for addiction risk. This work tries to provide new knowledge on this field by using two rat models of differential vulnerability to morphine addiction in order to study differential gene expression of a selected group of Hsp genes in the nucleus accumbens (NAc). Hsp70-1A, 84, 86 and 105 genes were similarly regulated by an acute injection of morphine in two subpopulations of Sprague Dawley (SD) rats showing different rates of extinction of morphine conditioned preference. However, Lewis and Fischer rats, two strains that differ in many aspects of drug seeking behaviours, exhibited marked differences in their expression patterns of Hsp84, 86 and 105. These results suggest that differential Hsp gene expression could be related to addiction vulnerability and recommend further work to validate these proteins as potential markers for drug addiction risk.

Characterisation of the heat shock factor of the human thermodimorphic pathogen Paracoccidioides lutzii.

Thermodimorphic fungi include most causative agents of systemic mycoses, but the molecular mechanisms that underlie their defining trait, i.e. the ability to shift between mould and yeast on temperature change alone, remain poorly understood. We hypothesised that the heat shock factor (Hsf), a protein that evolved to sense thermal stimuli quickly, might play a role in this process in addition to the known regulator Drk1 and the Ryp proteins. To test this hypothesis, we characterised the Hsf from the thermodimorph Paracoccidioides lutzii (formerly Paracoccidioides brasiliensis isolate 01). We show in the present work that PlHsf possesses regulatory domains that are exclusive of the Eurotiomycetidae family, suggesting evolutionary specialisation; that it can successfully rescue the otherwise lethal loss of the native protein of Saccharomyces cerevisiae; and that its DNA-binding domain is able to recognise regulatory elements from the promoters of both Drk1 and Ryp1. An in silico screening of all 1 kb sequences upstream of P. lutzii ORFs revealed that 7% of them possess a heat shock element. This is the first description of a heat shock factor in a thermodimorphic fungus.

Resolution-associated molecular patterns (RAMP): RAMParts defending immunological homeostasis?

The resolution of inflammation is central to the maintenance of good health and immune homeostasis. Recently, several intracellular stress proteins have been described as having extracellular properties that are anti-inflammatory or favour the resolution of inflammation. We propose that these molecules should be defined as resolution-associated molecular patterns (RAMPs). RAMPs are released at times of cellular stress and help to counterbalance the inflammatory effects of pathogen-associated (PAMPs) and damage-associated (DAMPs) molecular patterns. We propose that heat shock protein 10 (HSP10), αB-crystallin (αBC), HSP27 and binding immunoglobulin protein (BiP) should be considered founding members of the RAMP family. A greater understanding of RAMP biology may herald the development of novel immunotherapies.