Distribution and complementarity of hydropathy in multisubunit proteins

A survey of 40 multisubunit proteins and 2 protein-protein complexes was performed to assay quantitatively the distribution of hydropathy among the exterior surface, interior, contact surface, and noncontact exterior surface of the isolated subunits. We suggest a useful way to present this distribution by using a "hydropathy level diagram." Additionally, we have devised a function called "hydropathy complementarity" to quantitate the degree to which interacting surfaces have matching hydropathy distributions. Our survey revealed the following patterns: (1) The difference in hydropathy between the interior and exterior of subunits is a fairly invariant quantity. (2) On average, the hydropathy of the contact surface is higher than that of the exterior surface, but is not greater than that of the protein as a whole. There was variation, however, among the proteins. In some instances, the contact surface was more hydrophilic than the noncontact exterior, and in a few cases the contact surface was as hydrophobic as the protein interior. (3) The average interface manifests significant hydropathy complementarity, signifying that proteins interact by placing hydrophobic centers of one surface against hydrophobic centers of the other surface, and by similarly matching hydrophilic centers. As a measure of recognition and specificity, hydropathy complementarity could be a useful tool for predicting correct docking of interacting proteins. We suggest that high hydropathy complementarity is associated with static inflexible interactions. (4) We have found that some subunits that bind predominantly through hydrophilic forces, such as hydrogen bonds, ionic pairs, and water and metal bridges, are involved in dynamic quaternary organization and allostery.     

CARd: Carbon distribution analysis program for protein sequences

Carbon distribution is responsible for stability and structure of proteins. Arrangement of carbon along the protein sequence is depends on how the amino acids are organized and is guided by mRNAs. An atomic level revision is important for understanding these codes. This will ultimately help in identification of disorders and suggest mutations. For this purpose a carbon distribution analysis program has been developed. This program captures the hydrophobic / hydrophilic / disordered regions in a protein. The program gives accurate results. The calculations are precise and sensitive to single amino acid resolution. This program is to help in mutational studies leading to protein stabilisation.     

Wavelet transforms for the characterization and detection of repeating motifs

The role of repeating motifs in protein structures is thought to be as modular building blocks which allow an economic way of constructing complex proteins. In this work novel wavelet transform analysis techniques are used to detect and characterize repeating motifs in protein sequence and structure data, where the Kyte-Doolittle hydrophobicity scale (Eta Phi) and relative accessible surface area (rASA) data provide residue information about the protein sequence and structure, respectively. We analyze a variety of repeating protein motifs, TIM barrels, propellor blades, coiled coils and leucine-rich repeat structures. Detection and characterization of these motifs is performed using techniques based on the continuous wavelet transform (CWT). Results indicate that the wavelet transform techniques developed herein are a promising approach for the detection and characterization of repeating motifs for both structural and in some instances sequence data.     

SVA elements: a hominid-specific retroposon family

SVA is a composite repetitive element named after its main components, SINE, VNTR and Alu. We have identified 2762 SVA elements from the human genome draft sequence. Genomic distribution analysis indicates that the SVA elements are enriched in G+C-rich regions but have no preferences for inter- or intragenic regions. A phylogenetic analysis of the elements resulted in the recovery of six subfamilies that were named SVA_A to SVA_F. The composition, age and genomic distribution of the subfamilies have been examined. Subfamily age estimates based upon nucleotide divergence indicate that the expansion of four SVA subfamilies (SVA_A, SVA_B, SVA_C and SVA_D) began before the divergence of human, chimpanzee and gorilla, while subfamilies SVA_E and SVA_F are restricted to the human lineage. A survey of human genomic diversity associated with SVA_E and SVA_F subfamily members showed insertion polymorphism frequencies of 37.5% and 27.6%, respectively. In addition, we examined the amplification dynamics of SVA elements throughout the primate order and traced their origin back to the beginnings of hominid primate evolution, approximately 18 to 25 million years ago. This makes SVA elements the youngest family of retroposons in the primate order.     

Ribokinase family evolution and the role of conserved residues at the active site of the PfkB subfamily representative, Pfk-2 from Escherichia coli

Phosphofructokinase-2 (Pfk-2) belongs to the ribokinase family and catalyzes the ATP-dependent phosphorylation of fructose-6-phosphate, showing allosteric inhibition by a second ATP molecule. Several structures have been deposited on the PDB for this family of enzymes. A structure-based multiple sequence alignment of a non-redundant set of these proteins was used to infer phylogenetic relationships between family members with different specificities and to dissect between globally conserved positions and those common to phosphosugar kinases. We propose that phosphosugar kinases appeared early in the evolution of the ribokinase family. Also, we identified two conserved sequence motifs: the TR motif, not described previously, present in phosphosugar kinases but not in other members of the ribokinase family, and the globally conserved GXGD motif. Site-directed mutagenesis of R90 and D256 present in these motifs, indicate that R90 participates in the binding of the phosphorylated substrate and that D256 is involved in the phosphoryl transfer mechanism.     

In silico inference of inclusion membrane protein family in obligate intracellular parasites chlamydiae.

Chlamydiae are obligate intracellular pathogens that proliferate only inside a vacuole, called an inclusion. Chlamydial Inc proteins are known to be a major component of the inclusion membrane, but little is known about the gene number and function. The Inc proteins share very low sequence similarity but a similar hydropathy profile among them. Using the hydropathy profile, we computationally searched the open reading frames (ORFs) having a similar profile and predicted 90 and 36 ORFs (Inc-like ORFs) as candidates for Inc proteins in Chlamydia pneumoniae J138 and Chlamydia trachomatis serovar D, respectively. On the other hand, only a few Inc-like ORFs were found in organisms other than chlamydiae, suggesting that the Inc-like ORFs are specific to chlamydiae. Comparative genome analysis also revealed that the Inc-like ORFs have multiplied and diverged as paralogues and orthologues in the chlamydial genomes, and that some Inc-like ORFs lacked the N-terminal portion or encoded the split form. The data suggest that these gene products constitute a large protein family and may play an important role in chlamydial infection, growth and survival in the host cell.     

Classification of 29 families of secondary transport proteins into a single structural class using hydropathy profile analysis

A classification scheme for membrane proteins is proposed that clusters families of proteins into structural classes based on hydropathy profile analysis. The averaged hydropathy profiles of protein families are taken as fingerprints of the 3D structure of the proteins and, therefore, are able to detect more distant evolutionary relationships than amino acid sequences. A procedure was developed in which hydropathy profile analysis is used initially as a filter in a BLAST search of the NCBI protein database. The strength of the procedure is demonstrated by the classification of 29 families of secondary transporters into a single structural class, termed ST[3]. An exhaustive search of the database revealed that the 29 families contain 568 unique sequences. The proteins are predominantly from prokaryotic origin and most of the characterized transporters in ST[3] transport organic and inorganic anions and a smaller number are Na(+)/H(+) antiporters. All modes of energy coupling (symport, antiport, uniport) are found in structural class ST[3]. The relevance of the classification for structure/function prediction of uncharacterised transporters in the class is discussed.     

CARBANA: Carbon analysis program forprotein sequences

There are lots of works gone into proteins to understand the nature of proteins. Hydrophobic interaction is the dominant force that drives the proteins to carry out the biochemical reactions in all living system. Carbon is the only element that contributes towards this hydrophobic interaction. Studies find that globular proteins prefer to have 31.45% of carbon for its stability. Taking this as standard, a carbon analysis program has been developed to study the carbon distribution profile of protein sequences. This carbon analysis program has been made available online. This can be accessed at www.rajasekaran.net.in/tools/carbana.html. This new program is hoped to help in identification and development of active sites, study of protein stability, evolutionary understating of proteins, gene identification, ligand binding site identification, and to solve the long-standing problem of protein-protein and protein-DNA interactions.     

Prediction of function divergence in protein families using the substitution rate variation parameter alpha

Protein families typically embody a range of related functions and may thus be decomposed into subfamilies with, for example, distinct substrate specificities. Detection of functionally divergent subfamilies is possible by methods for recognizing branches of adaptive evolution in a gene tree. As the number of genome sequences is growing rapidly, it is highly desirable to automatically detect subfamily function divergence. To this end, we here introduce a method for large-scale prediction of function divergence within protein families. It is called the alpha shift measure (ASM) as it is based on detecting a shift in the shape parameter (alpha [alpha]) of the substitution rate gamma distribution. Four different methods for estimating alpha were investigated. We benchmarked the accuracy of ASM using function annotation from Enzyme Commission numbers within Pfam protein families divided into subfamilies by the automatic tree-based method BETE. In a test using 563 subfamily pairs in 162 families, ASM outperformed functional site-based methods using rate or conservation shifting (rate shift measure [RSM] and conservation shift measure [CSM]). The best results were obtained using the "GZ-Gamma" method for estimating alpha. By combining ASM with RSM and CSM using linear discriminant analysis, the prediction accuracy was further improved.     

Novel approach for structural identification of protein family: glyoxalase I

Glyoxalase is one of two enzymes of the glyoxalase detoxification system against methylglyoxal and other aldehydes, the metabolites derived from glycolysis. The glyoxalase system is found almost in all living organisms: bacteria, protozoa, plants, and animals, including humans, and is related to the class of ‘life essential proteins’. The enzyme belongs to the expanded Glyoxalase/Bleomycin resistance protein/Dioxygenase superfamily. At present the GenBank contains about 700 of amino acid sequences of this enzyme type, and the Protein Data Bank includes dozens of spatial structures. We have offered a novel approach for structural identification of glyoxalase I protein family, which is based on the selecting of basic representative proteins with known structures. On this basis, six new subfamilies of these enzymes have been derived. Most populated subfamilies A1 and A2 were based on representative human Homo sapiens and bacterial Escherichia coli enzymes. We have found that the principle feature, which defines the subfamilies’ structural differences, is conditioned by arrangement of N- and C-domains inside the protein monomer. Finely, we have deduced the structural classification for the glyoxalase I and assigned about 460 protein sequences distributed among six new subfamilies. Structural similarities and specific differences of all the subfamilies have been presented. This approach can be used for structural identification of thousands of the so-called hypothetical proteins with the known PDB structures allowing to identify many of already existing atomic coordinate entrees.     

Signal-anchor sequences are an essential factor for the Golgi-plasma membrane localization of type II membrane proteins

Despite studies of the mechanism underlying the intracellular localization of membrane proteins, the specific mechanisms by which each membrane protein localizes to the endoplasmic reticulum, Golgi apparatus, and plasma membrane in the secretory pathway are unclear. In this study, a discriminant analysis of endoplasmic reticulum, Golgi apparatus and plasma membrane-localized type II membrane proteins was performed using a position-specific scoring matrix derived from the amino acid propensity of the sequences around signal-anchors. The possibility that the sequence around the signal-anchor is a factor for identifying each localization group was evaluated. The discrimination accuracy between the Golgi apparatus and plasma membrane-localized type II membrane proteins was as high as 90%, indicating that, in addition to other factors, the sequence around signal-anchor is an essential component of the selection mechanism for the Golgi and plasma membrane localization. These results may improve the use of membrane proteins for drug delivery and therapeutic applications.     

2020年中国蜘蛛分类年报

本文总结了2020年发表的中国蜘蛛新分类单元和中国分类学者发表的国外蜘蛛新分类单元。基于国内标本, 2020年共发表中国蜘蛛223个新分类单元, 其中中国学者在12种期刊的56篇文章中命名了1新亚科、12新属和201新种, 国外学者命名了9新种, 合计210新种。云南是发现新种数量最多的省份, 占全国发现新种总数的38.1%。基于缅甸等12个国家标本, 2020年中国分类学者发表了外国蜘蛛95个新分类单元, 包括2新科、7新属和86新种。本文还总结了过去5年中国分类学者对全球蜘蛛分类的贡献, 2016‒2020年有1位中国学者发表771个新分类单元。     

Evolution of protein superfamilies and bacterial genome size

We present the structural annotation of 56 different bacterial species based on the assignment of genes to 816 evolutionary superfamilies in the CATH domain structure database. These assignments have enabled us to analyse the recurrence of specific superfamilies within and across the genomes. We have selected the superfamilies that have a very broad representation and therefore appear to be universally distributed in a significant number of bacterial lineages. Occurrence profiles of these universally distributed superfamilies are compared with genome size in order to estimate the correlation between superfamily duplication and the increase in proteome size. This distinguishes between those size-dependent superfamilies where frequency of occurrence is highly correlated with increase in genome size, and size-independent superfamilies where no correlation is observed. Consideration of the size correlation and the ratio between the mean and the standard deviations for all the superfamily profiles allows more detailed subdivisions and classification of superfamilies. For example, within the size-independent superfamilies, we distinguished a group that are distributed evenly amongst all the genomes. Within the size-dependent superfamilies we differentiated two groups: linearly distributed and non-linearly distributed. Functional annotation using the COG database was performed for all superfamilies in each of these groups, and this revealed significant differences amongst the three sets of superfamilies. Evenly distributed, size-independent domains are shown to be involved primarily in protein translation and biosynthesis. For the size-dependent superfamilies, linearly distributed superfamilies are involved mainly in metabolism, and non-linearly distributed superfamily domains are involved principally in gene regulation.     

Discrimination of Mammalian GPI-Anchored Proteins by Hydropathy and Amino Acid Propensities

A 1.7-kb DNA fragment cloned from Zymomonas mobilis genomic DNA complemented the inability to grow on sucrose of a Sue ^? mutant of Z. mobilis that was deficient in the production of both extracellular levansucrase and invertase. Analysis of the nucleotide sequence of the fragment found two open reading frames (ORFs), both of which did not correspond to the structural gene for the levansucrase or the invertase. By subcloning each ORF into two different Suc ^? mutants of Z. mobilis, it has been found that the first ORF (gene zliE) activates the production of the extracellular levansucrase and invertase, and the second ORF (gene zliS) stimulates the secretion of the two enzymes. Gene zliS might contribute to the secretion of proteins having no signal peptide. The expression of zliE and zliS seemed to be under the control of the same promoter.