Swaps in protein sequences

An important question in protein evolution is to what extent proteins may have undergone swaps (switches of domain or fragment order) during evolution. Such events might have occurred in several forms: Swaps of short fragments, swaps of structural and functional motifs, or recombination of domains in multidomain proteins. This question is important for the theoretical understanding of the evolution of proteins, and has practical implications for using swaps as a design tool in protein engineering. In order to analyze the question systematically, we conducted a large scale survey of possible swaps and permutations among all pairs of protein from the Swissport database. A swap is defined as a specific kind of sequence mutation between two proteins in which two fragments that appear in both sequences have different relative order in the two sequences. For example, aXbYc and dYeXf are defined as a swap, where X and Y represent sequence fragments that switched their order. Identifying such swaps is difficult using standard sequence comparison packages. One of the main problems in the analysis stems from the fact that many sequences contain repeats, which may be identified as false-positive swaps. We have used two different approaches to detect pairs of proteins with swaps. The first approach is based on the predefined list of domains in Pfam. We identified all the proteins that share at least two domains and analyzed their relative order, looking for pairs in which the order of these domains was switched. We designed an algorithm to distinguish between real swaps and duplications. In the second approach, we used Blast to detect pairs of proteins that share several fragments. Then, we used an automatic procedure to select pairs that are likely to contain swaps. Those pairs were analyzed visually, using a graphical tool, to eliminate duplications. Combining these approaches, about 140 different cases of swaps in the Swissprot database were found (after eliminating multiple pairs within the same family). Some of the cases have been described in the literature, but many are novel examples. Although each new example identified may be interesting to analyze, our main conclusion is that cases of swaps are rare in protein evolution. This observation is at odds with the common view that proteins are very modular to the point that modules (e.g., domains) can be shuffled between proteins with minimal constraints. Our study suggests that sequential constraints, i.e., the relative order between domains, are highly conserved.     

The order of PDZ3 and TrpCage in fusion chimeras determines their properties—a biophysical characterization

Most of the structural proteins known today are composed of domains that carry their own functions while keeping their structural properties. It is supposed that such domains, when taken out of the context of the whole protein, can retain their original structure and function to a certain extent. Information on the specific functional and structural characteristics of individual domains in a new context of artificial fusion proteins may help to reveal the rules of internal and external domain communication. Moreover, this could also help explain the mechanism of such communication and address how the mutual allosteric effect plays a role in a such multi‐domain protein system. The simple model system of the two‐domain fusion protein investigated in this work consisted of a well‐folded PDZ3 domain and an artificially designed small protein domain called Tryptophan Cage (TrpCage). Two fusion proteins with swapped domain order were designed to study their structural and functional features as well as their biophysical properties. The proteins composed of PDZ3 and TrpCage, both identical in amino acid sequence but different in composition (PDZ3‐TrpCage, TrpCage‐PDZ3), were studied using circualr dichroism (CD) spectrometry, analytical ultracentrifugation, and molecular dynamic simulations. The biophysical analysis uncovered different structural and denaturation properties of both studied proteins, revealing their different unfolding pathways and dynamics.     

Elastin-like polypeptide fusions enhance the accumulation of recombinant proteins in tobacco leaves

The production of recombinant proteins in plants is an active area of research and many different high-value proteins have now been produced in plants. Tobacco leaves have many advantages for recombinant protein production particularly since they allow field production without seeds, flowers or pollen and therefore provide for contained production. Despite these biosafety advantages recombinant protein accumulation in leaves still needs to be improved. Elastin-like polypeptides are repeats of the amino acids “VPGXG” that undergo a temperature dependant phase transition and have utility in the purification of recombinant proteins but can also enhance the accumulation of recombinant proteins they are fused to. We have used a 11.3 kDa elastin-like polypeptide as a fusion partner for three different target proteins, human interleukin-10, murine interleukin-4 and the native major ampullate spidroin protein 2 gene from the spider Nephila clavipes. In both transient analyses and stable transformants the concentrations of the fusion proteins were at least an order of magnitude higher for all of the fusion proteins when compared to the target protein alone. Therefore, fusions with a small ELP tag can be used to significantly enhance the accumulation of a range of different recombinant proteins in plant leaves. An erratum to this article can be found at     

Unfolding stabilities of two structurally similar proteins as probed by temperature-induced and force-induced molecular dynamics simulations

Unfolding stabilities of two homologous proteins, cardiotoxin III and short-neurotoxin (SNTX) belonging to three-finger toxin (TFT) superfamily, have been probed by means of molecular dynamics (MD) simulations. Combined analysis of data obtained from steered MD and all-atom MD simulations at various temperatures in near physiological conditions on the proteins suggested that overall structural stabilities of the two proteins were different from each other and the MD results are consistent with experimental data of the proteins reported in the literature. Rationalization for the differential structural stabilities of the structurally similar proteins has been chiefly attributed to the differences in the structural contacts between C- and N-termini regions in their three-dimensional structures, and the findings endorse the ‘CN network’ hypothesis proposed to qualitatively analyse the thermodynamic stabilities of proteins belonging to TFT superfamily of snake venoms. Moreover, the ‘CN network’ hypothesis has been revisited and the present study suggested that ‘CN network’ should be accounted in terms of ‘structural contacts’ and ‘structural strengths’ in order to precisely describe order of structural stabilities of TFTs.     

双孢蘑菇子实体发育不同阶段蛋白质组分析

为探讨双孢蘑菇子实体发育不同阶段的蛋白质表达变化,对双孢蘑菇As2796菌株子实体原基期、幼菇期、采收期、开伞期的蛋白质组进行了二维液相色谱串联质谱（iTRAQ-MS/MS）分析,共获得不同肽段5 869个,鉴定到1 059个蛋白质,其中1 007个具有相对定量信息。与双孢蘑菇原基期相比,幼菇期、采收期和开伞期分别有差异蛋白质242、200、240个,分别占鉴定蛋白质总数的24.0%,19.9%和23.8%。对这些蛋白质及不同阶段之间的差异蛋白质进行了系列生物信息学分析,对8个上、下调表达具有连续性的差异蛋白质相关基因进行了荧光定量PCR验证,其中3个蛋白质（错配碱基识别蛋白、细胞色素C1及某推定的未知蛋白质）基因的转录与蛋白质的表达具有较为一致的趋势。这些结果为后续双孢蘑菇子实体发育相关基因的功能研究奠定了良好的基础。     

蛇毒神经毒素类似物的免疫学分析

选取眼镜蛇蛇毒和银环蛇蛇毒中 1 0个差异较大的神经毒素类似物 ,通过麦芽糖融合表达系统使它们在大肠杆菌中得到高效可溶性表达 .通过 Amylose- Sepharose6B亲和树脂纯化这些融合蛋白 .分别制备这 1 0种融合蛋白多克隆抗体 .同时将这 1 0种神经毒素类似物进行硫氧还蛋白融合表达 ,表达产物基本上都是包涵体 .以硫氧还蛋白融合表达产物为抗原和上述 1 0种神经毒素类似物抗血清进行 Western blotting.结果显示 ,大多数神经毒素类似物之间没有免疫交叉反应 ,表明这些神经毒素类似物的抗原性差别较大 ,可能存在不同的分类和进化     

酵母双杂交系统在膜蛋白相互作用研究中的应用

膜相关蛋白约占细胞总蛋白质中的1/3,它们大都参与了细胞的诸多生理、病理过程和药物反应机理。研究膜蛋白的相互作用对于揭示细胞的生命活动规律及寻找药物作用靶标都有重要的意义。由于膜蛋白本身的特性及其难以进入核内等原因,经典的酵母双杂交技术并不适用于检测膜蛋白间的相互作用。针对在活细胞中研究膜蛋白相互作用的需要,近年来国际上先后发展了一系列用于膜蛋白相互作用研究的酵母双杂交新系统,并取得了许多重要发现。     

A detailed thermodynamic analysis of ras/effector complex interfaces

Many cellular functions are based on the interaction and crosstalk of various signaling proteins. Among these, members of the Ras family of small GTP-binding proteins are important for communicating signals into different pathways. In order to answer the question of how binding affinity and specificity is achieved, we analyzed binding energetics on the molecular level, with reference to the available structural data. The interaction of two members of the Ras subfamily with two different effector proteins, namely Raf and RalGDS, were investigated using isothermal titration calorimetry and a fluorescence-based method. Experiments with alanine mutants, located in the complex interfaces, yielded an energy map for the contact areas of the Ras/effector complexes, which could be differentiated into enthalpy and entropy contributions. In addition, by using double mutant cycle analysis, we probed the energetic contribution of selected pairs of amino acid residues. The resulting energy landscapes of the Ras/effector interface areas show a highly different topology when comparing the two effectors, Raf and RalGDS, demonstrating the specificity of the respective interactions. Particularly, we observe a high degree of compensating effects between enthalpy and entropy; differences between these components are much greater than the overall free energy differences. This is observed also when using the software FOLD-X to predict the effect of point mutations on the crystal structures of the different complexes. Prediction of the free energy changes shows a very good correlation with the experimentally observed energies. Furthermore, in line with experimental data, energy decomposition indicates that many different components of large magnitude counteract each other to produce a smaller change in overall free energy, illustrating the importance of long-range electrostatic forces in complex formation.     

Structural motifs in protein cores and at protein–protein interfaces are different

Structures of proteins and protein–protein complexes are determined by the same physical principles and thus share a number of similarities. At the same time, there could be differences because in order to function, proteins interact with other molecules, undergo conformations changes, and so forth, which might impose different restraints on the tertiary versus quaternary structures. This study focuses on structural properties of protein–protein interfaces in comparison with the protein core, based on the wealth of currently available structural data and new structure‐based approaches. The results showed that physicochemical characteristics, such as amino acid composition, residue–residue contact preferences, and hydrophilicity/hydrophobicity distributions, are similar in protein core and protein–protein interfaces. On the other hand, characteristics that reflect the evolutionary pressure, such as structural composition and packing, are largely different. The results provide important insight into fundamental properties of protein structure and function. At the same time, the results contribute to better understanding of the ways to dock proteins. Recent progress in predicting structures of individual proteins follows the advancement of deep learning techniques and new approaches to residue coevolution data. Protein core could potentially provide large amounts of data for application of the deep learning to docking. However, our results showed that the core motifs are significantly different from those at protein–protein interfaces, and thus may not be directly useful for docking. At the same time, such difference may help to overcome a major obstacle in application of the coevolutionary data to docking—discrimination of the intramolecular information not directly relevant to docking.     

Evolutionary analysis reveals collective properties and specificity in the C-type lectin and lectin-like domain superfamily

Members of the C-type lectin/C-type lectin-like domain (CTL/CTLD) superfamily share a common fold and are involved in a variety of functions, such as generalized defense mechanisms against foreign agents, discrimination between healthy and pathogen-infected cells, and endocytosis and blood coagulation. In this work we used ConSurf, a computer program recently developed in our lab, to perform an evolutionary analysis of this superfamily in order to further identify characteristics of all or part of its members. Given a set of homologous proteins in the form of multiple sequence alignment (MSA) and an inferred phylogenetic tree, ConSurf calculates the conservation score in every alignment position, taking into account the relationships between the sequences and the physicochemical similarity between the amino acids. The scores are then color-coded onto the three-dimensional structure of one of the homologous proteins. We provide here and at http://ashtoret.tau.ac.il/ approximately sharon a detailed analysis of the conservation pattern obtained for the entire superfamily and for two subgroups of proteins: (a) 21 CTLs and (b) 11 heterodimeric CTLD toxins. We show that, in general, proteins of the superfamily have one face that is constructed mostly of conserved residues and another that is not, and we suggest that the former face is involved in binding to other proteins or domains. In the CTLs examined we detected a region of highly conserved residues, corresponding to the known calcium- and carbohydrate-binding site of the family, which is not conserved throughout the entire superfamily, and in the CTLD toxins we found a patch of highly conserved residues, corresponding to the known dimerization region of these proteins. Our analysis also detected patches of conserved residues with yet unknown function(s).     

Matching fusion protein systems for affinity analysis of two interacting families of proteins: the cohesin-dockerin interaction

Cellulosomes are multi-enzyme complexes that orchestrate the efficient degradation of cellulose and related plant cell wall polysaccharides. The complex is maintained by the high-affinity protein-protein interaction between two complementary modules: the cohesin and the dockerin. In order to characterize the interaction between different cohesins and dockerins, we have developed matching fusion-protein systems, which harbor either the cohesin or the dockerin component. For this purpose, corresponding plasmid cassettes were designed, which encoded for the following carrier proteins: (i) a thermostable xylanase with an appended His-tag; and (ii) a highly stable cellulose-binding module (CBM). The resultant xylanase-dockerin and CBM-cohesin fusion products exhibited high expression levels of soluble protein. The expressed, affinity-purified proteins were extremely stable, and the functionality of the cohesin or dockerin component was retained. The fusion protein system was used to establish a sensitive and reliable, semi-quantitative enzyme-linked affinity assay for determining multiple samples of cohesin-dockerin interactions in microtiter plates. A variety of cohesin-dockerin systems, which had been examined previously using other methodologies, were revisited applying the affinity-based enzyme assay, the results of which served to verify the validity of the approach.     

Intrinsic disorder in protein sense‐antisense recognition

In addition to the well‐established sense‐antisense complementarity abundantly present in the nucleic acid world and serving as a basic principle of the specific double‐helical structure of DNA, production of mRNA, and genetic code‐based biosynthesis of proteins, sense‐antisense complementarity is also present in proteins, where sense and antisense peptides were shown to interact with each other with increased probability. In nucleic acids, sense‐antisense complementarity is achieved via the Watson‐Crick complementarity of the base pairs or nucleotide pairing. In proteins, the complementarity between sense and antisense peptides depends on a specific hydropathic pattern, where codons for hydrophilic and hydrophobic amino acids in a sense peptide are complemented by the codons for hydrophobic and hydrophilic amino acids in its antisense counterpart. We are showing here that in addition to this pattern of the complementary hydrophobicity, sense and antisense peptides are characterized by the complementary order‐disorder patterns and show complementarity in sequence distribution of their disorder‐based interaction sites. We also discuss how this order‐disorder complementarity can be related to protein evolution.     

Effect of microtubule-associated proteins on the protofilament number of microtubules assembled in vitro

In order to demonstrate the effect of microtubule-associated proteins on the protofilament number of microtubules, we used different systems of microtubule formation in vitro in which these proteins are either functionally eliminated (by DNA or glycerol) or absent (purified tubulin). The results obtained by electron microscopy of ultrathin-sectioned material indicate that under standard conditions in the presence of microtubule-associated proteins microtubules are formed consisting predominantly of 14 protofilaments. In cases of deficiency of microtubule-associated proteins, the mean value of the protofilament number is lower, and the protofilament number within the microtubule population varies remarkably. On the other hand, the action of microtubule-associated proteins is enhanced by histones resulting in increased protofilament numbers. A model is proposed illustrating that the quality and the quantity of microtubule-associated proteins bound to microtubules determine the curvature between the protofilaments and restrict the variety of their binding angles. In this way the microtubule-associated proteins may be regarded as an important factor in determining the structural fidelity of microtubules.     

Insights into protein structure and function from disorder-complexity space

Intrinsically disordered proteins have a wide variety of important functional roles. However, the relationship between sequence and function in these proteins is significantly different than that for well-folded proteins. In a previous work, we showed that the propensity to be disordered can be recognized based on sequence composition alone. Here that analysis is furthered by examining the relationship of disorder propensity to sequence complexity, where the metrics for these two properties depend only on composition. The distributions of 40 amino acid peptides from both ordered and disordered proteins are graphed in this disorder-complexity space. An analysis of Swiss-Prot shows that most peptides have high complexity and relatively low disorder. However, there are also an appreciable number of low complexity-high disorder peptides in the database. In contrast, there are no low complexity-low disorder peptides. A similar analysis for peptides in the PDB reveals a much narrower distribution, with few peptides of low complexity and high disorder. In this case, the bounds of the disorder-complexity distribution are well defined and might be used to evaluate the likelihood that a peptide can be crystallized with current methods. The disorder-complexity distributions of individual proteins and sets of proteins grouped by function are also examined. Among individual proteins, there is an enormous variety of distributions that in some cases can be rationalized with regard to function. Groups of functionally related proteins are found to have distributions that are similar within each group but show notable differences between groups. Finally, a pattern matching algorithm is used to search for proteins with particular disorder-complexity distributions. The results suggest that this approach might be used to identify relationships between otherwise dissimilar proteins.     

嗜热与嗜常温微生物的蛋白质氨基酸组成比较

嗜热微生物的嗜热特性与其蛋白质的高度热稳定性紧密相关。为了探索嗜热蛋白质的热稳定机制,比较嗜热和嗜常温微生物的蛋白质在氨基酸组成上的差别,收集110对分别来自嗜热和嗜常温微生物的同源蛋白质序列,比较两组蛋白质各种氨基酸含量以及疏水性氨基酸组成、疏水性指数和荷电氨基酸组成的差别,结果两者在多种氨基酸含量上存在微小但统计学上显著的差别,嗜热蛋白质比嗜常温蛋白质具有较高的平均疏水性和荷电氨基酸组成。对两组蛋白质的“脂肪族氨基酸指数”进行分析,证明嗜热蛋白质之所以具有较高的脂肪族氨基酸指数是由于其亮氨酸含量较高,与影响该指数的其它几种氨基酸无关;从而认为该指数的意义值得怀疑。通过对大量同源嗜热蛋白质和嗜常温蛋白质氨基酸组成的比较,能够揭示一些有关蛋白质热稳定性的普遍规律。     

Homology modeling in a dynamical world

A key concept in template‐based modeling (TBM) is the high correlation between sequence and structural divergence, with the practical consequence that homologous proteins that are similar at the sequence level will also be similar at the structural level. However, conformational diversity of the native state will reduce the correlation between structural and sequence divergence, because structural variation can appear without sequence diversity. In this work, we explore the impact that conformational diversity has on the relationship between structural and sequence divergence. We find that the extent of conformational diversity can be as high as the maximum structural divergence among families. Also, as expected, conformational diversity impairs the well‐established correlation between sequence and structural divergence, which is nosier than previously suggested. However, we found that this noise can be resolved using a priori information coming from the structure‐function relationship. We show that protein families with low conformational diversity show a well‐correlated relationship between sequence and structural divergence, which is severely reduced in proteins with larger conformational diversity. This lack of correlation could impair TBM results in highly dynamical proteins. Finally, we also find that the presence of order/disorder can provide useful beforehand information for better TBM performance.     

Insights into the local residual entropy of proteins provided by NMR relaxation.

A simple model is used to illustrate the relationship between the dynamics measured by NMR relaxation methods and the local residual entropy of proteins. The expected local dynamic behavior of well-packed extended amino acid side chains are described by employing a one-dimensional vibrator that encapsulates both the spatial and temporal character of the motion. This model is then related to entropy and to the generalized order parameter of the popular "model-free" treatment often used in the analysis of NMR relaxation data. Simulations indicate that order parameters observed for the methyl symmetry axes in, for example, human ubiquitin correspond to significant local entropies. These observations have obvious significance for the issue of the physical basis of protein structure, dynamics, and stability.     

Insertion of carrier proteins into hydrophilic loops of the Escherichia coli lactose permease

We describe the design and characterization of a set of fusion proteins of the Escherichia coli lactose (lac) permease in which a set of five different soluble “carrier” proteins (cytochrome_b562, flavodoxin, T4 lysozyme, β-lactamase and 70 kDa heat shock ATPase domain) were systematically inserted into selected loop positions of the transporter. The design goal was to increase the exposed hydrophilic surface area of the permease, while minimizing the internal flexibility of the resulting fusion proteins in order to improve the crystallization properties of the membrane protein. Fusion proteins with insertions into the central hydrophilic loop of the lac permease were active in transport lactose, although only the fusion proteins with E. coli cytochrome_b562, E. coli flavodoxin or T4 lysozyme were expressed at near wild-type lac permease levels. Eight other loop positions were tested with these three carriers, leading to the identification of additional fusion proteins that were active and well-expressed. By combining the results from the single carrier insertions, we have expressed functional “double fusion” proteins containing cytochrome_b562 domains inserted in two different loop positions.     

Aggregation of polyQ‐extended proteins is promoted by interaction with their natural coiled‐coil partners

Polyglutamine (polyQ) diseases are genetically inherited neurodegenerative disorders. They are caused by mutations that result in polyQ expansions of particular proteins. Mutant proteins form intranuclear aggregates, induce cytotoxicity and cause neuronal cell death. Protein interaction data suggest that polyQ regions modulate interactions between coiled‐coil (CC) domains. In the case of the polyQ disease spinocerebellar ataxia type‐1 (SCA1), interacting proteins with CC domains further enhance aggregation and toxicity of mutant ataxin‐1 (ATXN1). Here, we suggest that CC partners interacting with the polyQ region of a mutant protein, increase its aggregation while partners that interact with a different region reduce the formation of aggregates. Computational analysis of genetic screens revealed that CC‐rich proteins are highly enriched among genes that enhance pathogenicity of polyQ proteins, supporting our hypothesis. We therefore suggest that blocking interactions between mutant polyQ proteins and their CC partners might constitute a promising preventive strategy against neurodegeneration.     

Evolution of the isoelectric point of mammalian proteins as a consequence of indels and adaptive evolution

Although important shifts in the isoelectric point of prokaryotic proteins, mainly due to adaptation to environmental pH, have been widely reported, such studies have not covered mammalian proteins, where pH changes may relate to changes in subcellular or tissue compartmentalization. We explored the isoelectric point of the proteome of 13 mammalian species. We detected proteins that have shifted their pI the most among 13 mammalian species, and investigated if these differences reflect adaptations of the orthologous proteins to different conditions. We find that proteins exhibiting a high isoelectric point change are enriched in certain GO terms, including immune defense, and mitochondrial proteins. We show that the shift in pI between orthologous proteins is not strongly associated with the overall rate of protein evolution, nor with protein length. Our results reveal that insertions/deletions are the main reason behind the shift of pI. However, for some proteins we find evidence of selection shifting the pI of the protein through amino acid replacement. Finally, we argue that shifts in pI might relate to the gain of additional activities, such as new interacting partners, in one ortholog as opposed to the other, and may potentially relate to functional differences between mammals.