Identification of two hydrophilins that contribute to the desiccation and freezing tolerance of yeast (Saccharomyces cerevisiae) cells

Hydrophilins are a group of proteins that are present in all organisms and that have been defined as being highly hydrophilic and rich in glycine. They are assumed to play important roles in cellular dehydration tolerance. There are 12 genes in the yeast Saccharomyces cerevisiae that encode hydrophilins and most of these genes are stress responsive. However, the functional role of yeast hydrophilins, especially in desiccation and freezing tolerance, is largely unknown. Here, we selected six candidate hydrophilins for further analysis. All six proteins were predicted to be intrinsically disordered, i.e. to have no stable structure in solution. The contribution of these proteins to the desiccation and freezing tolerance of yeast was investigated in the respective knock-out strains. Only the disruption of the genes YJL144W and YMR175W (SIP18) resulted in significantly reduced desiccation tolerance, while none of the strains was affected in its freezing tolerance under our experimental conditions. Complementation experiments showed that yeast cells overexpressing these two genes were both more desiccation and freezing tolerant, confirming the role of these two hydrophilins in yeast dehydration stress tolerance.     

Functional dissection of hydrophilins during in vitro freeze protection

In plants, Late Embryogenesis Abundant (LEA) proteins typically accumulate in response to low water availability conditions imposed during development or by the environment. Analogous proteins in other organisms are induced when exposed to stress conditions. Most of this diverse set of proteins can be grouped according to properties such as high hydrophilicity and high content of glycine or other small amino acids in what we have termed hydrophilins. Previously, we showed that hydrophilins protect enzyme activities in vitro from low water availability effects. Here, we demonstrate that hydrophilins can also protect enzyme activities from the adverse effects induced by freeze-thaw cycles in vitro. We monitored conformational changes induced by freeze-thaw on the enzyme lactate dehydrogenase (LDH) using the fluorophore 1-anilinonaphthalene-8-sulfonate (ANS). Hydrophilin addition prevents enzyme inactivation and this effect is reflected in changes in the ANS-fluorescence levels determined for LDH. We further show that for selected plant hydrophilins, removal of certain conserved domains affects their protecting capabilities. Thus, we propose that hydrophilins, and in particular specific protein domains, have a role in protecting cell components from the adverse effects caused by low water availability such as those present during freezing conditions by preventing deleterious changes in protein secondary and tertiary structure.     

Characterization of protein hubs by inferring interacting motifs from protein interactions

The characterization of protein interactions is essential for understanding biological systems. While genome-scale methods are available for identifying interacting proteins, they do not pinpoint the interacting motifs (e.g., a domain, sequence segments, a binding site, or a set of residues). Here, we develop and apply a method for delineating the interacting motifs of hub proteins (i.e., highly connected proteins). The method relies on the observation that proteins with common interaction partners tend to interact with these partners through a common interacting motif. The sole input for the method are binary protein interactions; neither sequence nor structure information is needed. The approach is evaluated by comparing the inferred interacting motifs with domain families defined for 368 proteins in the Structural Classification of Proteins (SCOP). The positive predictive value of the method for detecting proteins with common SCOP families is 75% at sensitivity of 10%. Most of the inferred interacting motifs were significantly associated with sequence patterns, which could be responsible for the common interactions. We find that yeast hubs with multiple interacting motifs are more likely to be essential than hubs with one or two interacting motifs, thus rationalizing the previously observed correlation between essentiality and the number of interacting partners of a protein. We also find that yeast hubs with multiple interacting motifs evolve slower than the average protein, contrary to the hubs with one or two interacting motifs. The proposed method will help us discover unknown interacting motifs and provide biological insights about protein hubs and their roles in interaction networks.     

Plant dehydrins and stress tolerance: Versatile proteins for complex mechanisms

Dehydrins (DHNs), or group 2 LEA (Late Embryogenesis Abundant) proteins, play a fundamental role in plant response and adaptation to abiotic stresses. They accumulate typically in maturing seeds or are induced in vegetative tissues following salinity, dehydration, cold and freezing stress. The generally accepted classification of dehydrins is based on their structural features, such as the presence of conserved sequences, designated as Y, S and K segments. The K segment representing a highly conserved 15 amino acid motif forming amphiphilic a-helix is especially important since it has been found in all dehydrins. Since more than 20 y, they are thought to play an important protective role during cellular dehydration but their precise function remains unclear. This review outlines the current status of the progress made toward the structural, physico-chemical and functional characterization of plant dehydrins and how these features could be exploited in improving stress tolerance in plants.Key words: abiotic stress, dehydration stress, drought, cold acclimation, freezing tolerance, LEA proteins, dehydrins     

Codep: maximizing co-evolutionary interdependencies to discover interacting proteins

Approaches for the determination of interacting partners from different protein families (such as ligands and their receptors) have made use of the property that interacting proteins follow similar patterns and relative rates of evolution. Interacting protein partners can then be predicted from the similarity of their phylogenetic trees or evolutionary distances matrices. We present a novel method called Codep, for the determination of interacting protein partners by maximizing co-evolutionary signals. The order of sequences in the multiple sequence alignments from two protein families is determined in such a manner as to maximize the similarity of substitution patterns at amino acid sites in the two alignments and, thus, phylogenetic congruency. This is achieved by maximizing the total number of interdependencies of amino acids sites between the alignments. Once ordered, the corresponding sequences in the two alignments indicate the predicted interacting partners. We demonstrate the efficacy of this approach with computer simulations and in analyses of several protein families. A program implementing our method, Codep, is freely available to academic users from our website: http://www.uhnresearch.ca/labs/tillier/.     

On the Regularities of the Polar Profiles of Proteins Related to Ebola Virus Infection and their Functional Domains

The number of fatalities and economic losses caused by the Ebola virus infection across the planet culminated in the havoc that occurred between August and November 2014. However, little is known about the molecular protein profile of this devastating virus. This work represents a thorough bioinformatics analysis of the regularities of charge distribution (polar profiles) in two groups of proteins and their functional domains associated with Ebola virus disease: Ebola virus proteins and Human proteins interacting with Ebola virus. Our analysis reveals that a fragment exists in each of these proteins—one named the “functional domain”—with the polar profile similar to the polar profile of the protein that contains it. Each protein is formed by a group of short sub-sequences, where each fragment has a different and distinctive polar profile and where the polar profile between adjacent short sub-sequences changes orderly and gradually to coincide with the polar profile of the whole protein. When using the charge distribution as a metric, it was observed that it effectively discriminates the proteins from their functional domains. As a counterexample, the same test was applied to a set of synthetic proteins built for that purpose, revealing that any of the regularities reported here for the Ebola virus proteins and human proteins interacting with Ebola virus were not present in the synthetic proteins. Our results indicate that the polar profile of each protein studied and its corresponding functional domain are similar. Thus, when building each protein from its functional domai—adding one amino acid at a time and plotting each time its polar profile—it was observed that the resulting graphs can be divided into groups with similar polar profiles.     

Discovery of Unique Lanthionine Synthetases Reveals New Mechanistic and Evolutionary Insights

Lantibiotic synthetases are remarkable biocatalysts generating conformationally constrained peptides with a variety of biological activities by repeatedly utilizing two simple posttranslational modification reactions: dehydration of Ser/Thr residues and intramolecular addition of Cys thiols to the resulting dehydro amino acids. Since previously reported lantibiotic synthetases show no apparent homology with any other known protein families, the molecular mechanisms and evolutionary origin of these enzymes are unknown. In this study, we present a novel class of lanthionine synthetases, termed LanL, that consist of three distinct catalytic domains and demonstrate in vitro enzyme activity of a family member from Streptomyces venezuelae. Analysis of individually expressed and purified domains shows that LanL enzymes install dehydroamino acids via phosphorylation of Ser/Thr residues by a protein kinase domain and subsequent elimination of the phosphate by a phosphoSer/Thr lyase domain. The latter has sequence homology with the phosphothreonine lyases found in various pathogenic bacteria that inactivate host mitogen activated protein kinases. A LanC-like cyclase domain then catalyzes the addition of Cys residues to the dehydro amino acids to form the characteristic thioether rings. We propose that LanL enzymes have evolved from stand-alone protein Ser/Thr kinases, phosphoSer/Thr lyases, and enzymes catalyzing thiol alkylation. We also demonstrate that the genes for all three pathways to lanthionine-containing peptides are widespread in Nature. Given the remarkable efficiency of formation of lanthionine-containing polycyclic peptides and the latter''s high degree of specificity for their cognate cellular targets, it is perhaps not surprising that (at least) three distinct families of polypeptide sequences have evolved to access this structurally and functionally diverse class of compounds.     

On the accuracy of inferring energetic coupling between distant sites in protein families from evolutionary imprints: Illustrations using lattice model

It is suspected that correlated motions among a subset of spatially separated residues drive conformational dynamics not only in multidomain but also in single domain proteins. Sequence and structure‐based methods have been proposed to determine covariation between two sites on a protein. The statistical coupling analysis (SCA) that compares the changes in probability at two sites in a multiple sequence alignment (MSA) and a subset of the MSA has been used to infer the network of residues that encodes allosteric signals in protein families. The structural perturbation method (SPM), that probes the response of a local perturbation at all other sites, has been used to probe the allostery wiring diagram in biological machines and enzymes. To assess the efficacy of the SCA, we used an exactly soluble two dimensional lattice model and performed double‐mutant cycle (DMC) calculations to predict the extent of physical coupling between two sites. The predictions of the SCA and the DMC results show that only residues that are in contact in the native state are accurately identified. In addition, covariations among strongly interacting residues are most easily identified by the SCA. These conclusions are consistent with the DMC experiments on the PDZ family. Good correlation between the SCA and the DMC is only obtained by performing multiple experiments that vary the nature of amino acids at a given site. In contrast, the energetic coupling found in experiments for the PDZ domain are recovered using the SPM. We also predict, using the SPM, several residues that are coupled energetically. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

The origins of specificity in polyketide synthase protein interactions

Polyketides, a diverse group of heteropolymers with antibiotic and antitumor properties, are assembled in bacteria by multiprotein chains of modular polyketide synthase (PKS) proteins. Specific protein-protein interactions determine the order of proteins within a multiprotein chain, and thereby the order in which chemically distinct monomers are added to the growing polyketide product. Here we investigate the evolutionary and molecular origins of protein interaction specificity. We focus on the short, conserved N- and C-terminal docking domains that mediate interactions between modular PKS proteins. Our computational analysis, which combines protein sequence data with experimental protein interaction data, reveals a hierarchical interaction specificity code. PKS docking domains are descended from a single ancestral interacting pair, but have split into three phylogenetic classes that are mutually noninteracting. Specificity within one such compatibility class is determined by a few key residues, which can be used to define compatibility subclasses. We identify these residues using a novel, highly sensitive co-evolution detection algorithm called CRoSS (correlated residues of statistical significance). The residue pairs selected by CRoSS are involved in direct physical interactions in a docked-domain NMR structure. A single PKS system can use docking domain pairs from multiple classes, as well as domain pairs from multiple subclasses of any given class. The termini of individual proteins are frequently shuffled, but docking domain pairs straddling two interacting proteins are linked as an evolutionary module. The hierarchical and modular organization of the specificity code is intimately related to the processes by which bacteria generate new PKS pathways.     

Identification of unique SUN-interacting nuclear envelope proteins with diverse functions in plants

Although a plethora of nuclear envelope (NE) transmembrane proteins (NETs) have been identified in opisthokonts, plant NETs are largely unknown. The only known NET homologues in plants are Sad1/UNC-84 (SUN) proteins, which bind Klarsicht/ANC-1/Syne-1 homology (KASH) proteins. Therefore, de novo identification of plant NETs is necessary. Based on similarities between opisthokont KASH proteins and the only known plant KASH proteins, WPP domain–interacting proteins, we used a computational method to identify the KASH subset of plant NETs. Ten potential plant KASH protein families were identified, and five candidates from four of these families were verified for their NE localization, depending on SUN domain interaction. Of those, Arabidopsis thaliana SINE1 is involved in actin-dependent nuclear positioning in guard cells, whereas its paralogue SINE2 contributes to innate immunity against an oomycete pathogen. This study dramatically expands our knowledge of plant KASH proteins and suggests that plants and opisthokonts have recruited different KASH proteins to perform NE regulatory functions.     

Highly hydrophilic proteins in prokaryotes and eukaryotes are common during conditions of water deficit

The late embryogenesis abundant (LEA) proteins are plant proteins that are synthesized at the onset of desiccation in maturing seeds and in vegetative organs exposed to water deficit. Here, we show that most LEA proteins are comprised in a more widespread group, which we call "hydrophilins." The defining characteristics of hydrophilins are high glycine content (>6%) and a high hydrophilicity index (>1.0). By data base searching, we show that this criterion selectively differentiates most known LEA proteins as well as additional proteins from different taxons. We found that within the genomes of Escherichia coli and Saccharomyces cerevisiae, only 5 and 12 proteins, respectively, meet our criterion. Despite their deceivingly loose definition, hydrophilins usually represent <0.2% of the proteins of a genome. Additionally, we demonstrate that the criterion that defines hydrophilins seems to be an excellent predictor of responsiveness to hyperosmosis since most of the genes encoding these proteins in E. coli and S. cerevisiae are induced by osmotic stress. Evidence for the participation of one of the E. coli hydrophilins in the adaptive response to hyperosmotic conditions is presented. Apparently, hydrophilins represent analogous adaptations to a common problem in such diverse taxons as prokaryotes and eukaryotes.     

Polarity and hydropathic properties of natural amino acids

A new classification of amino acids according to their polarity and symmetric location in the spatial structure of the genetic code is suggested. The polar amino acids are: R, S (codons AGC and AGU), K, N, Q, H, W, C, Y, G, E, D; apolar ones are: T, M, I, P, L, S (codons UCN). Polar and apolar amino acids are grouped into three families whose members possess complementarity with respect to the symmetric structure of the genetic code. Interaction of these complementary polar and apolar amino acids encodes formation of the space structures and ligand-receptor complexes of proteins. Correlation between the polar and hydropathic properties of amino acids is investigated. Normalization of 38 hydrophobicity scales of natural amino acids is carried out. A discrepancy between structures of polar/hydrophilic and apolar/hydrophobic groups of amino acids is demonstrated. According to the signature principle this discrepancy is due to different properties of amino acid side radicals which, in turn, depend on the second component of the reaction and on environmental conditions.     

Protein stability to the effect of denaturing agents and elevated temperature is directly proportional to the fraction of hydrophobic residues in its molecule

The dependence of protein stability against the action of urea, guanidine hydrochloride and elevated temperature on the fraction of hydrophobic residues in their molecules was studied. It was shown that proteins can be divided into several stability classes. The stability of proteins within each class is directly proportional to the fraction of hydrophobic residues in their molecules and this dependence is a linear one. As shown on several examples proteins can pass from one class to another and such discrete transitions are determined by the formation or disruption of bonds with cofactors. The effect of point mutations on the stability of proteins is discussed from the viewpoint of polar group dehydration.     

Two novel heat-soluble protein families abundantly expressed in an anhydrobiotic tardigrade

Tardigrades are able to tolerate almost complete dehydration by reversibly switching to an ametabolic state. This ability is called anhydrobiosis. In the anhydrobiotic state, tardigrades can withstand various extreme environments including space, but their molecular basis remains largely unknown. Late embryogenesis abundant (LEA) proteins are heat-soluble proteins and can prevent protein-aggregation in dehydrated conditions in other anhydrobiotic organisms, but their relevance to tardigrade anhydrobiosis is not clarified. In this study, we focused on the heat-soluble property characteristic of LEA proteins and conducted heat-soluble proteomics using an anhydrobiotic tardigrade. Our heat-soluble proteomics identified five abundant heat-soluble proteins. All of them showed no sequence similarity with LEA proteins and formed two novel protein families with distinct subcellular localizations. We named them Cytoplasmic Abundant Heat Soluble (CAHS) and Secretory Abundant Heat Soluble (SAHS) protein families, according to their localization. Both protein families were conserved among tardigrades, but not found in other phyla. Although CAHS protein was intrinsically unstructured and SAHS protein was rich in β-structure in the hydrated condition, proteins in both families changed their conformation to an α-helical structure in water-deficient conditions as LEA proteins do. Two conserved repeats of 19-mer motifs in CAHS proteins were capable to form amphiphilic stripes in α-helices, suggesting their roles as molecular shield in water-deficient condition, though charge distribution pattern in α-helices were different between CAHS and LEA proteins. Tardigrades might have evolved novel protein families with a heat-soluble property and this study revealed a novel repertoire of major heat-soluble proteins in these anhydrobiotic animals.     

Structure model of core proteins in photosystem I inferred from the comparison with those in photosystem II and bacteria; an application of principal component analysis to detect the similar regions between distantly related families of proteins.

A principal component analysis based on the physico-chemical properties of amino acid residues is developed to assign similar regions between distantly related families of proteins, taking account of the species diversities in respective families. The most important advantage of this analysis should be that it reflects different physico-chemical properties and thus can predict more detailed structural properties, including the transmembrane helices, than the hydropathy analysis. Its first application reconfirms the similarity between the core proteins of photosynthetic reaction center in purple bacteria and those of photosystem II, indicating that the low percentage of identical amino acid residues estimated previously between them is due to much allowance for amino acid substitutions in purple bacteria. The application of this analysis to the core proteins of photosystem I reveals that any of these proteins includes two domains, each showing high similarity to the amino acid sequences of core proteins in photosystem II and purple bacteria. A core structure model of A1 and A2 proteins folded into four layers of sheets of transmembrane helices is proposed to provide a molecular basis for the electron pathway suggested by spectroscopic experiments as well as for the interaction sites with plastocyanin, 9 kDa protein and LHC proteins.     

Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.