Search for ancient patterns in protein sequences

Proteins of related functions are often similar in sequence, reflecting a common phylogenetic origin. Proteins with no known homology are probably diversified proteins, too distantly related to known sequences in databases to retain significant similarity. All proteins, however, probably share common ancestries if one moves far enough back in evolution; therefore, given the huge accumulation of protein sequences in current databases, it could be expected that some proteins with no obvious sequence resemblance to any other share some residues that could represent footprints of ancient common ancestries. To identify such putative footprints, we have searched for short stretches of amino acids present in a given protein sequence that are also found in a significant number of nonrelated proteins in the database. The significantly high frequency of occurrence of these patterns in the database would support a common evolutionary source, and a diversity of non-related proteins that contain the pattern would express their ancient origin. Using this strategy, significant patterns were found in actual exons, but not in randomized amino acid sequences, nor in translated sequences of noncoding DNA, suggesting that this strategy actually leads to the identification of patterns with a biological significance. These significant patterns are not randomly positioned along the sequences analyzed, but they tend to accumulate within specific regions, producing a profile of discrete domains. In some well-known proteins analyzed in this study, some of these domains are coincident with known motifs. Thus, the procedure described in this paper could be useful for identifying ancient patterns and domains in protein sequences, some of which could also have a functional or structural significance.     

Migration patterns of proteins in ancient human bones detected by SDS electrophoresis

Purified and ultraconcentrated extracts of 7 burial remains from different areas and various periods (from IV B.C. to 1300–1400 A.D.) were repeatedly submitted to electrophoresis in SDS. Five bands: A, B, C, D, E (from the slowest to the fastest) were identified altogether. Band B appears to be the most resistant, because it is present in all the samples with reactivity. Proteins corresponding to the five bands have molecular weights comprised between 65000 and 45000 in agreement with other authors.     

Methods for mapping of interaction networks involving membrane proteins

Nearly one-third of all genes in various organisms encode membrane-associated proteins that participate in numerous protein-protein interactions important to the processes of life. However, membrane protein interactions pose significant challenges due to the need to solubilize membranes without disrupting protein-protein interactions. Traditionally, analysis of isolated protein complexes by high-resolution 2D gel electrophoresis has been the main method used to obtain an overall picture of proteome constituents and interactions. However, this method is time consuming, labor intensive, detects only abundant proteins and is limited with respect to the coverage required to elucidate large interaction networks. In this review, we discuss the application of various methods to elucidate interactions involving membrane proteins. These techniques include methods for the direct isolation of single complexes or interactors as well as methods for characterization of entire subcellular and cellular interactomes.     

The effect of heterogeneous structural diffusion on ligand binding to heme proteins

The kinetics of ligand binding to heme proteins studied by flash photolysis display an algebraic time dependence at low temperatures in contrast exponential recombination observed under physiological conditions. This result shows that protein structures should be viewed as a time average of interconverting microstates which are frozen in at low temperatures. We propose a quasi-one-dimensional model of heterogeneous structural diffusion coupled to ligand binding which describes freezing transition as an inherent property of protein fluctuations. The structural hopping rates are derived from a temperature invariant spectrum of activation energies. The model predicts power law kinetics of the form t ^- at long times. The exponent is constant (0.5) at high temperatures but decreases below a critical temperature in the frozen regime. These results are compared to experiments performed with myoglobin and -chains of hemoglobin.     

Molecular modelling of miraculin: Structural analyses and functional hypotheses

Miraculin is a plant protein that displays the peculiar property of modifying taste by swiching sour into a sweet taste. Its monomer is flavourless at all pH as well as at high concentration; the dimer form elicits its taste-modifying activity at acidic pH; a tetrameric form is also reported as active. Two histidine residues, located in exposed regions, are the main responsible of miraculin activity, as demonstrated by mutagenesis studies. Since structural data of miraculin are not available, we have predicted its three-dimensional structure and simulated both its dimer and tetramer forms by comparative modelling and molecular docking techniques. Finally, molecular dynamics simulations at different pH conditions have indicated that at acidic pH the dimer assumes a widely open conformation, in agreement with the hypotheses coming from other studies.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Prediction of Catalytic Residues Using the Variation of Stereochemical Properties

In this paper, we investigate a simple protein sequence conservation measure which takes amino acid similarity into account. Instead of grouping 20 amino acids into disjoint sets in previous methods, we consider ten overlapping classes. The method is based on the assumption that a column in a multiple sequence alignment is evolved from an identical column in the evolutionary history. Two ten-dimensional vectors are constructed for each position to denote frequencies of ten classes in a column and the corresponding hypothetical identical column. Then the cosine function of the angle between these two vectors is considered as a measure of divergence of stereochemical properties at this position. This divergence, combining with other conservation scores, is used as conservation measure of the column. Finally, we evaluate our methods by identifying catalytic sites, using rank analysis criterion and receiver operator characteristic analysis criterion.     

Mutational and sequence analysis of transmembrane segment 6 orientation in TetA proteins

The packing orientations of the 8 transmembrane (TM) segments that line the central, aqueous transport channel within tetracycline resistance proteins (TetA) have been established. However, the orientations of the remaining 4 segments, TMs 3, 6, 9, and 12, located at the periphery, and away from the transport channel, have not yet been determined. In this study, the packing orientation of TM6 within the class C TetA protein encoded by plasmid pBR322 was evaluated by substitution mutagenesis and analysis of sequence conservation and amphipathicity. The combined data support a model in which the conserved and polar face of the TM6 alpha-helix containing Asn170 and Asn173 orients towards channel-lining TM segments, and the relatively non-conserved and hydrophobic face of TM6 points towards membrane lipids.     

Origins of globular structure in proteins

Since natural proteins are the products of a long evolutionary process, the structural properties of present-day proteins should depend not only on physico-chemical constraints, but also on evolutionary constraints. Here we propose a model for protein evolution, in which membranes play a key role as a scaffold for supporting the gradual evolution from flexible polypeptides to well-folded proteins. We suggest that the folding process of present-day globular proteins is a relic of this putative evolutionary process. To test the hypothesis that membranes once acted as a cradle for the folding of globular proteins, extensive research on membrane proteins and the interactions of globular proteins with membranes will be required.     

CLAP: A web-server for automatic classification of proteins with special reference to multi-domain proteins

Background

The function of a protein can be deciphered with higher accuracy from its structure than from its amino acid sequence. Due to the huge gap in the available protein sequence and structural space, tools that can generate functionally homogeneous clusters using only the sequence information, hold great importance. For this, traditional alignment-based tools work well in most cases and clustering is performed on the basis of sequence similarity. But, in the case of multi-domain proteins, the alignment quality might be poor due to varied lengths of the proteins, domain shuffling or circular permutations. Multi-domain proteins are ubiquitous in nature, hence alignment-free tools, which overcome the shortcomings of alignment-based protein comparison methods, are required. Further, existing tools classify proteins using only domain-level information and hence miss out on the information encoded in the tethered regions or accessory domains. Our method, on the other hand, takes into account the full-length sequence of a protein, consolidating the complete sequence information to understand a given protein better.

Results

Our web-server, CLAP (Classification of Proteins), is one such alignment-free software for automatic classification of protein sequences. It utilizes a pattern-matching algorithm that assigns local matching scores (LMS) to residues that are a part of the matched patterns between two sequences being compared. CLAP works on full-length sequences and does not require prior domain definitions.Pilot studies undertaken previously on protein kinases and immunoglobulins have shown that CLAP yields clusters, which have high functional and domain architectural similarity. Moreover, parsing at a statistically determined cut-off resulted in clusters that corroborated with the sub-family level classification of that particular domain family.

Conclusions

CLAP is a useful protein-clustering tool, independent of domain assignment, domain order, sequence length and domain diversity. Our method can be used for any set of protein sequences, yielding functionally relevant clusters with high domain architectural homogeneity. The CLAP web server is freely available for academic use at http://nslab.mbu.iisc.ernet.in/clap/.     

GnRH-Bik/Bax/Bak chimeric proteins target and kill adenocarcinoma cells; the general use of pro-apoptotic proteins of the Bcl-2 family as novel killing components of targeting chimeric proteins

In recent years chimeric proteins carrying bacterial toxins as their killing moiety, have been developed to selectively recognize and kill cell populations expressing speciific receptors. The involvement of Gonadotropin releasing hormone (GnRH) has been demonstrated in several adenocarcinomas and a GnRH-bacterial toxin chimeric protein (GnRH-PE66) was thus developed and found to specifically target and kill adenocarcinoma cells both in vitro and in vivo. Because of the immunogenicity and the non-specific toxicity of the bacterial toxins, we have developed new chimeric proteins, introducing apoptosis inducing proteins of the Bcl-2 family as novel killing components. Sequences encoding the human Bik, Bak or Bax proteins were fused to the GnRH coding sequence at the DNA level and were expressed in E. coli. GnRH-Bik, GnRH-Bak and GnRH-Bax new chimeric proteins efficiently and specifically inhibited the cell growth of adenocarcinoma cell lines and eventually led to cell death. All three Bcl2-proteins-based chimeric proteins seem to induce apoptosis within the target cells, without any additional cell death stimulus. Apoptosis-inducing-proteins of the Bcl-2 family targeted by the GnRH are novel potential therapeutic reagents for adenocarcinoma treatment in humans. This novel approach could be widely applied, using any molecule that binds a specific cell type, fused to an apoptosis-inducing protein.     

Periplasmic proteins of Escherichia coli are highly resistant to aggregation: reappraisal for roles of molecular chaperones in periplasm

Periplasmic proteins of Gram-negative bacteria like Escherichia coli are subjected to immediate affect of environmental fluctuation that may unfold proteins, due to the permeability of the outer membrane to small molecules. They are thus supposedly protected by certain molecular chaperones. Nevertheless, no homologues of typical molecular chaperones have so far been found in periplasm, and the recently reported chaperone activities of periplasmic protein disulfide isomerase (PDI) and peptidyl prolyl isomerase (PPI) seem to be too weak to satisfy such assumed needs. In an attempt to reveal whether periplasmic proteins exhibit certain unusual properties, we discovered that such proteins as a whole are highly resistant to aggregation under a wide variety of denaturing conditions. Furthermore, in an effort to unveil the nature behind this phenomenon we purified and examined four prominent periplasmic proteins. Our results demonstrate that these proteins unfold at rather mild denaturing conditions and expose hydrophobic surfaces during such unfolding process, but hardly form complexes with a typical molecular chaperone. Based on these observations, we propose that the periplasmic proteins have been evolved to resist the formation of aggregates when subjected to various denaturing conditions and molecular chaperones may thus not be needed in periplasm.     

Differential degradation for small heat shock proteins IbpA and IbpB is synchronized in Escherichia coli: Implications for their functional cooperation in substrate refolding

Small heat shock proteins (sHSPs), as a conserved family of ATP-independent molecular chaperones, are known to bind non-native substrate proteins and facilitate the substrate refolding in cooperation with ATP-dependent chaperones (e.g., DnaK and ClpB). However, how different sHSPs function in coordination is poorly understood. Here we report that IbpA and IbpB, the two sHSPs of Escherichia coli, are coordinated by synchronizing their differential in vivo degradation. Whereas the individually expressed IbpA and IbpB are respectively degraded slowly and rapidly in cells cultured under both heat shock and normal conditions, their simultaneous expression leads to a synchronized degradation at a moderate rate. Apparently, such synchronization is linked to their hetero-oligomerization and cooperation in binding substrate proteins. In addition, truncation of the flexible N- and C-terminal tails dramatically suppresses the IbpB degradation, and somehow accelerates the IbpA degradation. In view of these in vivo data, we propose that the synchronized degradation for IbpA and IbpB are crucial for their synergistic promoting effect on DnaK/ClpB-mediated substrate refolding, conceivably via the formation of IbpA–IbpB-substrate complexes. This scenario may be common for different sHSPs that interact with each other in cells.     

Amplification of the ancient murine Lx family of long interspersed repeated DNA occurred during the murine radiation

We identified and characterized the relics of an ancient rodent Ll family, referred to as Lx, which was extensively amplified at the time of the murine radiation about 12 million years ago, and which we showed was ancestral to the modern L1 families in rat and mouse. Here we have extended our analysis of the Lx amplification by examining more murine and nonmurine species for Lx sequences using both blot hybridization and the polymerase chain reaction for a total of 36 species. In addition we have determined the relative copy number and sequence divergence, or age, of Lx elements in representative murine genera. Our results show that while Lx sequences are confined to murine genera, the extent of the amplification was different in the different murine lineages, indicating that the amplification of Lx did not precede, but was coincident with, the murine radiation. The implications of our findings for the evolutionary dynamics of L1 families and the utility of ancestral amplification events for systematics are discussed. Correspondence to: A.V. Furano     

Interaction surface and topology of Get3-Get4-Get5 protein complex, involved in targeting tail-anchored proteins to endoplasmic reticulum

Recent work has uncovered the "GET system," which is responsible for endoplasmic reticulum targeting of tail-anchored proteins. Although structural information and the individual roles of most components of this system have been defined, the interactions and interplay between them remain to be elucidated. Here, we investigated the interactions between Get3 and the Get4-Get5 complex from Saccharomyces cerevisiae. We show that Get3 interacts with Get4-Get5 via an interface dominated by electrostatic forces. Using isothermal titration calorimetry and small-angle x-ray scattering, we further demonstrate that the Get3 homodimer interacts with two copies of the Get4-Get5 complex to form an extended conformation in solution.     

Purification of antifreeze proteins by adsorption to ice

Antifreeze proteins (AFPs) can protect organisms from freezing injury by adsorbing to ice and inhibiting its growth. We describe here a method where ice, grown on a cold finger, is used to selectively adsorb and purify these ice-binding proteins from a crude mixture. Type III recombinant AFP was enriched approximately 50-fold after one round of partitioning into ice and purified to homogeneity by a second round. This method can also be used to purify non-ice-binding proteins by linkage to AFP domains as demonstrated by the recovery of a 50 kDa maltose-binding protein-AFP fusion from a crude lysate of Escherichia coli.     

Self-assembly of repeat proteins: Concepts and design of new interfaces

In nature, assembled protein structures offer the most complex functional structures. The understanding of the mechanisms ruling protein–protein interactions opens the door to manipulate protein assemblies in a rational way. Proteins are versatile scaffolds with great potential as tools in nanotechnology and biomedicine because of their chemical, structural, and functional versatility. Currently, bottom-up self-assembly based on biomolecular interactions of small and well-defined components, is an attractive approach to biomolecular engineering and biomaterial design. Specifically, repeat proteins are simplified systems for this purpose.In this work, we provide an overview of fundamental concepts of the design of new protein interfaces. We describe an experimental approach to form higher order architectures by a bottom-up assembly of repeated building blocks. For this purpose, we use designed consensus tetratricopeptide repeat proteins (CTPRs). CTPR arrays contain multiple identical repeats that interact through a single inter-repeat interface to form elongated superhelices. Introducing a novel interface along the CTPR superhelix allows two CTPR molecules to assemble into protein nanotubes. We apply three approaches to form protein nanotubes: electrostatic interactions, hydrophobic interactions, and π-π interactions. We isolate and characterize the stability and shape of the formed dimers and analyze the nanotube formation considering the energy of the interaction and the structure in the three different models. These studies provide insights into the design of novel protein interfaces for the control of the assembly into more complex structures, which will open the door to the rational design of nanostructures and ordered materials for many potential applications in nanotechnology.