Pattern similarity study of functional sites in protein sequences: lysozymes and cystatins

Background  

Although it is generally agreed that topography is more conserved than sequences, proteins sharing the same fold can have different functions, while there are protein families with low sequence similarity. An alternative method for profile analysis of characteristic conserved positions of the motifs within the 3D structures may be needed for functional annotation of protein sequences. Using the approach of quantitative structure-activity relationships (QSAR), we have proposed a new algorithm for postulating functional mechanisms on the basis of pattern similarity and average of property values of side-chains in segments within sequences. This approach was used to search for functional sites of proteins belonging to the lysozyme and cystatin families.     

Structural adaptation of the subunit interface of oligomeric thermophilic and hyperthermophilic enzymes

Enzymes from thermophilic and, particularly, from hyperthermophilic organisms are surprisingly stable. Understanding of the molecular origin of protein thermostability and thermoactivity attracted the interest of many scientist both for the perspective comprehension of the principles of protein structure and for the possible biotechnological applications through application of protein engineering. Comparative studies at sequence and structure levels were aimed at detecting significant differences of structural parameters related to protein stability between thermophilic and hyperhermophilic structures and their mesophilic homologs. Comparative studies were useful in the identification of a few recurrent themes which the evolution utilized in different combinations in different protein families. These studies were mostly carried out at the monomer level. However, maintenance of a proper quaternary structure is an essential prerequisite for a functional macromolecule. At the environmental temperatures experienced typically by hyper- and thermophiles, the subunit interactions mediated by the interface must be sufficiently stable. Our analysis was therefore aimed at the identification of the molecular strategies adopted by evolution to enhance interface thermostability of oligomeric enzymes. The variation of several structural properties related to protein stability were tested at the subunit interfaces of thermophilic and hyperthermophilic oligomers. The differences of the interface structural features observed between the hyperthermophilic and thermophilic enzymes were compared with the differences of the same properties calculated from pairwise comparisons of oligomeric mesophilic proteins contained in a reference dataset. The significance of the observed differences of structural properties was measured by a t-test. Ion pairs and hydrogen bonds do not vary significantly while hydrophobic contact area increases specially in hyperthermophilic interfaces. Interface compactness also appears to increase in the hyperthermophilic proteins. Variations of amino acid composition at the interfaces reflects the variation of the interface properties.     

Computing motif correlations in proteins

Protein motifs, which are specific regions and conserved regions, are found by comparing multiple protein sequences. These conserved regions in general play an important role in protein functions and protein folds, for example, for their binding properties or enzymatic activities. The aim here is to find the existence correlations of protein motifs. The knowledge of protein motif/domain sharing should be important in shedding new light on the biologic functions of proteins and offering a basis in analyzing the evolution in the human genome or other genomes. The protein sequences used here are obtained from the PIR-NREF database and the protein motifs are retrieved from the PROSITE database. We apply data mining approach to discover the occurrence correlations of motif in protein sequences. The correlation of motifs mined can be used in evolution analyses and protein structure prediction. We discuss the latter, i.e., protein structure prediction in this study. The correlations mined are stored and maintained in a database system. The database is now available at http://bioinfo.csie.ncu.edu.tw/ProMotif/.     

Analysis of the relationships between evolvability,thermodynamics, and the functions of intrinsically disordered proteins/regions

The evolvability of proteins is not only restricted by functional and structural importance, but also by other factors such as gene duplication, protein stability, and an organism's robustness. Recently, intrinsically disordered proteins (IDPs)/regions (IDRs) have been suggested to play a role in facilitating protein evolution. However, the mechanisms by which this occurs remain largely unknown. To address this, we have systematically analyzed the relationship between the evolvability, stability, and function of IDPs/IDRs. Evolutionary analysis shows that more recently emerged IDRs have higher evolutionary rates with more functional constraints relaxed (or experiencing more positive selection), and that this may have caused accelerated evolution in the flanking regions and in the whole protein. A systematic analysis of observed stability changes due to single amino acid mutations in IDRs and ordered regions shows that while most mutations induce a destabilizing effect in proteins, mutations in IDRs cause smaller stability changes than in ordered regions. The weaker impact of mutations in IDRs on protein stability may have advantages for protein evolvability in the gain of new functions. Interestingly, however, an analysis of functional motifs in the PROSITE and ELM databases showed that motifs in IDRs are more conserved, characterized by smaller entropy and lower evolutionary rate, than in ordered regions. This apparently opposing evolutionary effect may be partly due to the flexible nature of motifs in IDRs, which require some key amino acid residues to engage in tighter interactions with other molecules. Our study suggests that the unique conformational and thermodynamic characteristics of IDPs/IDRs play an important role in the evolvability of proteins to gain new functions.     

Molecular phylogenetic study and expression analysis of ATP-binding cassette transporter gene family in Oryza sativa in response to salt stress

ATP-binding cassette (ABC) transporter is a large gene superfamily that utilizes the energy released from ATP hydrolysis for transporting myriad of substrates across the biological membranes. Although many investigations have been done on the structural and functional analysis of the ABC transporters in Oryza sativa, much less is known about molecular phylogenetic and global expression pattern of the complete ABC family in rice. In this study, we have carried out a comprehensive phylogenetic analysis constructing neighbor-joining and maximum-likelihood trees based on various statistical methods of different ABC protein subfamily of five plant lineages including Chlamydomonas reinhardtii (green algae), Physcomitrella patens (moss), Selaginella moellendorffii (lycophyte), Arabidopsis thaliana (dicot) and O. sativa (monocot) to explore the origin and evolutionary patterns of these ABC genes. We have identified several conserved motifs in nucleotide binding domain (NBD) of ABC proteins among all plant lineages during evolution. Amongst the different ABC protein subfamilies, ‘ABCE’ has not yet been identified in lower plant genomes (algae, moss and lycophytes). The result indicated that gene duplication and diversification process acted upon these genes as a major operative force creating new groups and subgroups and functional divergence during evolution. We have demonstrated that rice ABCI subfamily consists of only half size transporters that represented highly dynamic members showing maximum sequence variations among the other rice ABC subfamilies. The evolutionary and the expression analysis contribute to a deep insight into the evolution and diversity of rice ABC proteins and their roles in response to salt stress that facilitate our further understanding on rice ABC transporters.     

The archaeal Sac10b protein family: conserved proteins with divergent functions

Here we review the present state of structural and functional studies of the Sac10b protein family, a class of highly conserved 10 kDa nucleic acid-binding proteins in archaea. Based on biochemical and structural studies, these proteins were originally assigned a role in the structural organization of chromatin; Sac10b proteins of hyperthermophilic archaea, for example, showed tight, unspecific DNA binding. More recently, however, Sac10b proteins of mesophilic archaea were found to interact preferentially with specific DNA sequences thereby affecting the expression of distinct genes. Furthermore, Sac10b proteins of hyperthermophilic, thermophilic and mesophilic archaea were also shown to bind to RNA with distinct affinities and specificities but functional consequences of RNA binding of these proteins, besides perhaps RNA stabilization, have not yet been observed. To better understand the physiological meaning of the various interactions of Sac10b proteins with nucleic acids, future work should concentrate on elucidating the molecular structures of complexes of Sac10b proteins of hyperthermophilic and mesophilic archaea with DNA and RNA. In addition, existing and new X-ray and NMR structures of individual hyperthermophilic Sac10b proteins may represent very good models for introducing thermostability especially in enzymes for industrial use.     

Subsequence-based feature map for protein function classification

Automated classification of proteins is indispensable for further in vivo investigation of excessive number of unknown sequences generated by large scale molecular biology techniques. This study describes a discriminative system based on feature space mapping, called subsequence profile map (SPMap) for functional classification of protein sequences. SPMap takes into account the information coming from the subsequences of a protein. A group of protein sequences that belong to the same level of classification is decomposed into fixed-length subsequences and they are clustered to obtain a representative feature space mapping. Mapping is defined as the distribution of the subsequences of a protein sequence over these clusters. The resulting feature space representation is used to train discriminative classifiers for functional families. The aim of this approach is to incorporate information coming from important subregions that are conserved over a family of proteins while avoiding the difficult task of explicit motif identification. The performance of the method was assessed through tests on various protein classification tasks. Our results showed that SPMap is capable of high accuracy classification in most of these tasks. Furthermore SPMap is fast and scalable enough to handle large datasets.     

Water in Ras Superfamily Evolution

The Ras GTPase superfamily of proteins coordinates a diverse set of cellular outcomes, including cell morphology, vesicle transport, and cell proliferation. Primary amino acid sequence analysis has identified Specificity determinant positions (SDPs) that drive diversified functions specific to the Ras, Rho, Rab, and Arf subfamilies (Rojas et al. 2012, J Cell Biol 196 :189–201). The inclusion of water molecules in structural and functional adaptation is likely to be a major response to the selection pressures that drive evolution, yet hydration patterns are not included in phylogenetic analysis. This article shows that conserved crystallographic water molecules coevolved with SDP residues in the differentiation of proteins within the Ras superfamily of small GTPases. The patterns of water conservation between protein subfamilies parallel those of sequence-based evolutionary trees. Thus, hydration patterns have the potential to help elucidate functional significance in the evolution of amino acid residues observed in phylogenetic analysis of homologous proteins. © 2019 Wiley Periodicals, Inc.     

The chaperonin of the archaeonSulfolobus solfataricus

The chaperonin of the hyperthermophilic archaeonSulfolobus solfataricus, briefly Ssocpn, was purified by a fast and high-yield procedure. Ssocpn, a 920 kDa-complex of two different subunits, displays a potassium-dependent ATPase activity with a temperature optimum at 80°C. The ability of Ssocpn to function in vitro was investigated using different protein substrates. Ssocpn promotes correct refolding of thermophilic and mesophilic enzymes from their chemically unfolded state; moreover, Ssocpn prevents the irreversible inactivation of native proteins by suppressing their precipitation upon heating. Both the activity in assisting refolding of unfolded proteins and that in preventing heat denaturation of native proteins require the hydrolysis of ATP. The chaperone-based strategies in different technological fields are discussed, and the advantages in using archaeal chaperonins are underlined.     

Analysis of sequence repeats of proteins in the PDB

Internal repeats in protein sequences play a significant role in the evolution of protein structure and function. Applications of different bioinformatics tools help in the identification and characterization of these repeats. In the present study, we analyzed sequence repeats in a non-redundant set of proteins available in the Protein Data Bank (PDB). We used RADAR for detecting internal repeats in a protein, PDBeFOLD for assessing structural similarity, PDBsum for finding functional involvement and Pfam for domain assignment of the repeats in a protein. Through the analysis of sequence repeats, we found that identity of the sequence repeats falls in the range of 20–40% and, the superimposed structures of the most of the sequence repeats maintain similar overall folding. Analysis sequence repeats at the functional level reveals that most of the sequence repeats are involved in the function of the protein through functionally involved residues in the repeat regions. We also found that sequence repeats in single and two domain proteins often contained conserved sequence motifs for the function of the domain.     

The sequence homologies of cytochromes P-450 and active-site geometries

Summary The amino acid sequence alignment of 16 cytochrome P-450 proteins representative of the major families is reported. The sequence matching process has been carried out on the basis of maximum homology by residue type, retention of secondary structure and minimization of deletions/insertions except where additional loop regions exist. From the starting point of known reported sequence homology matching from the literature, a realignment on the basis of conserved residues involved in both structure and function gives rise to a self-consistent set of sequences which correlates with known mechanistic and structural data. Once fitted, these archetypal sequences form a straightforward template for the alignment of all P-450 subfamilies. Computer modelling of the active-site regions constructed from homology with the bacterial form of the enzyme (P-450_CAM) evinces the correct substrate specificity. Furthermore, the construction of the macromolecular assembly of components of the cytochrome P-450 system on the microsomal endoplasmic reticular membrane is presented from the evidence of site-directed mutagenesis, analysis by molecular probes, X-ray crystallography and molecular modelling.     

HomologyPlot: Searching for homology to a family of proteins using a database of unique conserved patterns

Summary A new database of conserved amino acid residues is derived from the multiple sequence alignment of over 84 families of protein sequences that have been reported in the literature. This database contains sequences of conserved hydrophobic core patterns which are probably important for structure and function, since they are conserved for most sequences in that family. This database differs from other single-motif or signature databases reported previously, since it contains multiple patterns for each family. The new database is used to align a new sequence with the conserved regions of a family. This is analogous to reports in the literature where multiple sequence alignments are used to improve a sequence alignment. A program called Homology-Plot (suitable for IBM or compatible computers) uses this database to find homology of a new sequence to a family of protein sequences. There are several advantages to using multiple patterns. First, the program correctly identifies a new sequence as a member of a known family. Second, the search of the entire database is rapid and requires less than one minute. This is similar to performing a multiple sequence alignment of a new sequence to all of the known protein family sequences. Third, the alignment of a new sequence to family members is reliable and can reproduce the alignment of conserved regions already described in the literature. The speed and efficiency of this method is enhanced, since there is no need to score for insertions or deletions as is done in the more commonly used sequence alignment methods. In this method only the patterns are aligned. HomologyPlot also provides general information on each family, as well as a listing of patterns in a family.     

Molecular evolution of the plant ECERIFERUM1 and ECERIFERUM3 genes involved in aliphatic hydrocarbon production

The Arabidopsis ECERIFERUM1 (CER1) protein is a decarbonylase that converts fatty acid metabolites into alkanes. Alkanes are components of waxes in the plant cuticle, a waterproof barrier serving to protect land plants from both biotic and abiotic stimuli. CER1 enzymes can be used to produce alternative and sustainable hydrocarbons in eukaryotic systems. In this report we identified 193 CER1 and 128 CER3 sequences from 56 land plants respectively. CER1 and CER3 proteins have high amino acid similarity and both are involved in alkane synthesis in Arabidopsis. The common homologues of CER1 and CER3 genes were identified in three species of chlorophytes, which may be one of the earliest plant taxa that possess CER1 and CER3 genes. To facilitate potential applications, the 3-dimensional structure and conserved motifs of CER1 proteins were also characterized. CER1 and CER3 proteins are structurally similar, but CER1 proteins have more conserved histidine-containing motifs common to fatty acid hydroxylases and stearoyl-CoA desaturases. There was no significant loss or gain of protein motifs after ancient and recent duplications, suggesting that varied properties of CER1 proteins may be associated with less-conserved regions. Among 56 land plants, the codon-based assessments of selection modes revealed that neither entire proteins nor individual amino acids of CER1 proteins were significantly subjected to positive selection, indicating that CER1 proteins are highly conserved throughout evolution.     

基于氨基酸模糊聚类分析的跨膜区域预测

跨膜蛋白在进化过程中,序列保守性较差,即使是同源蛋白序列的一致性程度也较低,因而在跨膜区预测算法中,通过序列的一致性程度来选取训练集并不能有效地消除预测结果对训练集的过度适应性．本文提出了一种基于氨基酸模糊聚类分析的预测算法,通过氨基酸在各个区域分布的相似性程度进行模糊聚类,从而根据一类氨基酸的分布特性而不是各个氨基酸的分布特性进行跨膜区预测．结果表明,该方法能在一定程度上消除训练集的选取对测试结果的影响,提高跨膜蛋白拓扑结构预测的准确度,特别是提高对目前知之甚少的跨膜蛋白的预测准确度．     

Targeting the Small GTPase Superfamily through Their Regulatory Proteins

The Ras superfamily of small GTPases are guanine‐nucleotide‐dependent switches essential for numerous cellular processes. Mutations or dysregulation of these proteins are associated with many diseases, but unsuccessful attempts to target the small GTPases directly have resulted in them being classed as “undruggable”. The GTP‐dependent signaling of these proteins is controlled by their regulators; guanine nucleotide exchange factors (GEFs), GTPase activating proteins (GAPs), and in the Rho and Rab subfamilies, guanine nucleotide dissociation inhibitors (GDIs). This review covers the recent small molecule and biologics strategies to target the small GTPases through their regulators. It seeks to critically re‐evaluate recent chemical biology practice, such as the presence of PAINs motifs and the cell‐based readout using compounds that are weakly potent or of unknown specificity. It highlights the vast scope of potential approaches for targeting the small GTPases in the future through their regulatory proteins.     

Complete thermal-unfolding profiles of oxidized and reduced cytochromes C

The complete thermal-unfolding profiles of both oxidized and reduced forms of cytochrome c551 (PA) from mesophilic Pseudomonas aeruginosa and cytochrome c552 (HT) from thermophilic Hydrogenobacter thermophilus were obtained by the newly developed pressure-proof cell compartment installed in a circular dichroic spectrometer, which facilitates protein thermal-unfolding experiments up to 180 degrees C. The thermodynamic cycle, which relates protein stability and redox function, indicated that the redox potentials of PA and HT in the native state are regulated by the stability of the oxidized proteins rather than by that of the reduced ones.     

Protein splicing in cis and in trans

Intein-mediated protein splicing is a self-catalytic process in which the intervening intein sequence is removed from a precursor protein and the flanking extein segments are ligated with a native peptide bond. Splice junction proximal residues and internal residues within the intein direct these reactions. The identity of these residues varies in each intein, as groups of related residues populate conserved motifs. Although the basics of the four-step protein splicing pathway are known, mechanistic details are still unknown. Structural and kinetic analyses are beginning to shed some light. Several structures were reported for precursor proteins with mutations in catalytic residues, which stabilize the precursors for crystallographic study. Progress is being made despite limitations inherent in using mutated precursors. However, no uniform mechanism has emerged. Kinetic parameters were determined using conditional trans-splicing (splicing of split precursor fragments after intein reassembly). Several groups concluded that the rate of the initial acyl rearrangement step is rapid and Asn cyclization (step 3) is slow, suggesting that this latter step is rate limiting. Understanding the protein splicing pathway has allowed scientists to harness inteins for numerous applications.     

SiteBinder: an improved approach for comparing multiple protein structural motifs

There is a paramount need to develop new techniques and tools that will extract as much information as possible from the ever growing repository of protein 3D structures. We report here on the development of a software tool for the multiple superimposition of large sets of protein structural motifs. Our superimposition methodology performs a systematic search for the atom pairing that provides the best fit. During this search, the RMSD values for all chemically relevant pairings are calculated by quaternion algebra. The number of evaluated pairings is markedly decreased by using PDB annotations for atoms. This approach guarantees that the best fit will be found and can be applied even when sequence similarity is low or does not exist at all. We have implemented this methodology in the Web application SiteBinder, which is able to process up to thousands of protein structural motifs in a very short time, and which provides an intuitive and user-friendly interface. Our benchmarking analysis has shown the robustness, efficiency, and versatility of our methodology and its implementation by the successful superimposition of 1000 experimentally determined structures for each of 32 eukaryotic linear motifs. We also demonstrate the applicability of SiteBinder using three case studies. We first compared the structures of 61 PA-IIL sugar binding sites containing nine different sugars, and we found that the sugar binding sites of PA-IIL and its mutants have a conserved structure despite their binding different sugars. We then superimposed over 300 zinc finger central motifs and revealed that the molecular structure in the vicinity of the Zn atom is highly conserved. Finally, we superimposed 12 BH3 domains from pro-apoptotic proteins. Our findings come to support the hypothesis that there is a structural basis for the functional segregation of BH3-only proteins into activators and enablers.