Search

Iron homeostasis is vital for the host's defense against pathogenic invasion and the ferritinophagy is a crucial mechanism in maintaining intracellular iron homeostasis by facilitating the degradation and recycling of stored iron. The nuclear receptor coactivator 4 (NCOA4) serves as a ferritinophagy receptor, facilitating the binding and delivery of ferritin to the autophagosome and lysosome. However, NCOA4 of the sea cucumber Apostichopus japonicus (AjNCOA4) has not been reported until now. In this study, we identified and characterized AjNCOA4 in A. japonicus. This gene encodes a polypeptide containing 597 amino acids with an open reading frame of 1794 bp. The inferred amino acid sequence of AjNCOA4 comprises an ARA70 domain. Furthermore, a multiple sequence alignment demonstrated varying degrees of sequence homology between AjNCOA4 from A. japonicus and other NCOA4 orthologs. The phylogenetic tree of NCOA4 correlates with the established timeline of metazoan evolution. Expression analysis revealed that AjNCOA4 is expressed in all tested tissues, including the body wall, muscle, intestine, respiratory tree, and coelomocytes. Following challenge with Vibrio splendidus, the coelomocytes exhibited a significant increase in AjNCOA4 mRNA levels, peaking at 24 h. We successfully obtained recombinant AjNCOA4 protein through prokaryotic expression and prepared a specific polyclonal antibody. Immunofluorescence and co-immunoprecipitation experiments demonstrated an interaction between AjNCOA4 and AjFerritin in coelomocytes. RNA interference-mediated knockdown of AjNCOA4 expression resulted in elevated iron ion levels in coelomocytes. Bacterial stimulation enhanced ferritinophagy in coelomocytes, while knockdown of AjNCOA4 reduced the occurrence of ferritinophagy. These findings suggest that AjNCOA4 modulates ferritinophagy induced by V. splendidus in coelomocytes of A. japonicus., Competing Interests: Declaration of competing interest The authors declare no competing financial interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

C-type lectins (CTLs) are an important class of pattern recognition receptors (PRRs) that exhibit structural and functional diversity in invertebrates. Repetitive DNA sequences are ubiquitous in eukaryotic genomes, representing distinct modes of genome evolution and promoting new gene generation. Our study revealed a new CTL that is composed of two long tandem repeats, abundant threonine, and one carbohydrate recognition domain (CRD) in Exopalaemon carinicauda and has been designated EcTR-CTL. The full-length cDNA of EcTR-CTL was 1242 bp long and had an open reading frame (ORF) of 999 bp that encoded a protein of 332 amino acids. The genome structure of EcTR-CTL contains 4 exons and 3 introns. The length of each repeat unit in EcTR-CTL was 198 bp, which is different from the short tandem repeats reported previously in prawns and crayfish. EcTR-CTL was abundantly expressed in the intestine and hemocytes. After Vibrio parahaemolyticus and white spot syndrome virus (WSSV) challenge, the expression level of EcTR-CTL in the intestine was upregulated. Knockdown of EcTR-CTL downregulated the expression of anti-lipopolysaccharide factor, crustin, and lysozyme during Vibrio infection. The recombinant CRD of EcTR-CTL (rCRD) could bind to bacteria, lipopolysaccharides, and peptidoglycans. Additionally, rCRD can directly bind to WSSV. These findings indicate that 1) CTLs with tandem repeats may be ubiquitous in crustaceans, 2) EcTR-CTL may act as a PRR to participate in the innate immune defense against bacteria via nonself-recognition and antimicrobial peptide regulation, and 3) EcTR-CTL may play a positive or negative role in the process of WSSV infection by capturing virions., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

H2A.Z, the most evolutionarily conserved variant of histone H2A, plays a pivotal role in chromatin remodeling and contributes significantly to gene transcription and genome stability. However, the role of H2A.Z in the silkworm (Bombyx mori) remains unclear. In this study, we cloned the BmH2A.Z from B. mori. The open reading frame of BmH2A.Z is 390 bp, encoding 129 amino acids, with a confirmed molecular weight of 13.4 kDa through prokaryotic expression analysis. Sequence analysis revealed that BmH2A.Z has a conserved H2A.Z domain and is closely related to the systemic evolution of other known H2A.Zs. The expression profile of BmH2A.Z at various developmental stages of the B. mori exhibited the highest expression level in the 1st instar, followed by the grain stage and the 2nd instar, and the lowest expression level in the moth. The highest transcript level of BmH2A.Z was observed in the head, with relatively lower levels detected in the blood than in the other tissues under consideration. In addition, the upregulation of BmH2A.Z resulted in the amplified expression of B. mori nucleopolyhedrovirus (BmNPV) genes, thus facilitating the proliferation of BmNPV. This study establishes a foundation for investigating the role of BmH2A.Z in B. mori and its participation in virus-host interactions., (© 2024 Wiley Periodicals LLC.)

Echinococcus granulosus sensu lato is a platyhelminth parasite and the etiological cause of cystic echinococcosis (CE), a zoonotic and neglected disease that infects animals and humans worldwide. As a part of the biological arsenal of the parasite, cathepsin L proteases are a group of proteins that are believed to be essential for parasite penetration, immune evasion, and establishment in the tissues of the host. In this work, we have cloned and sequenced a new putative cathepsin L protease from Echinococcus canadensis (EcCLP1). The bioinformatic analysis suggests that EcCLP1 could be synthesized as a zymogen and activated after proteolytic cleavage. The multiple sequence alignment with other cathepsin proteases reveals important functional conserved features like a conserved active site, an N-linked glycosylation residue, a catalytic triad, an oxyanion hole, and three putative disulfide bonds. The phylogenetic analysis suggests that EcCLP1 could indeed be a cathepsin L cysteine protease from clade 1 as it grouped with cathepsins from other species in this clade. Modeling studies suggest that EcCLP1 has two domains forming a cleft where the active site is located and an occluding role for the propeptide. The transcriptomic analysis reveals different levels of cathepsin transcript expression along the different stages of the parasite life cycle. The whole-mount immunohistochemistry shows an interesting superficial punctate pattern of staining which suggests a secretory pattern of expression. The putative cathepsin L protease characterized here may represent an interesting tool for diagnostic purposes, vaccine design, or a new pharmacological target for antiparasitic intervention., (© A. Naidich et al., published by EDP Sciences, 2024.)

Protein language models, trained on millions of biologically observed sequences, generate feature-rich numerical representations of protein sequences. These representations, called sequence embeddings, can infer structure-functional properties, despite protein language models being trained on primary sequence alone. While sequence embeddings have been applied toward tasks such as structure and function prediction, applications toward alignment-free sequence classification have been hindered by the lack of studies to derive, quantify and evaluate relationships between protein sequence embeddings. Here, we develop workflows and visualization methods for the classification of protein families using sequence embedding derived from protein language models. A benchmark of manifold visualization methods reveals that Neighbor Joining (NJ) embedding trees are highly effective in capturing global structure while achieving similar performance in capturing local structure compared with popular dimensionality reduction techniques such as t-SNE and UMAP. The statistical significance of hierarchical clusters on a tree is evaluated by resampling embeddings using a variational autoencoder (VAE). We demonstrate the application of our methods in the classification of two well-studied enzyme superfamilies, phosphatases and protein kinases. Our embedding-based classifications remain consistent with and extend upon previously published sequence alignment-based classifications. We also propose a new hierarchical classification for the S-Adenosyl-L-Methionine (SAM) enzyme superfamily which has been difficult to classify using traditional alignment-based approaches. Beyond applications in sequence classification, our results further suggest NJ trees are a promising general method for visualizing high-dimensional data sets., (© The Author(s) 2023. Published by Oxford University Press.)

Malaria is a life-threatening disease caused by one of the five species of Plasmodium, among which Plasmodium falciparum cause the deadliest form of the disease. Plasmodium species are dependent on a vertebrate host and a blood-sucking insect vector to complete their life cycle. Plasmodium chitinases belonging to the GH18 family are secreted inside the mosquito midgut, during the ookinete stage of the parasite. Chitinases mediate the penetration of parasite through the peritrophic membrane, facilitating access to the gut epithelial layer. In this review, we describe Plasmodium chitinases with special emphasis on chitinases from P. falciparum and P. vivax, the representative examples of the short and long forms of this protein. In addition to the chitinase domain, chitinases belonging to the long form contain a pro-domain and chitin-binding domain. Amino acid sequence alignment of long and short form chitinase domains reveals multiple positions containing variant residues. A subset of these positions was found to be conserved or invariant within long or short forms, indicating the role of these positions in attributing form-specific activity. The reported differences in affinities to allosamidin for P. vivax and P. falciparum were predicted to be due to different residues at two amino acid positions, resulting in altered interactions with the inhibitor. Understanding the role of these amino acids in Plasmodium chitinases will help us elucidate the mechanism of catalysis and the mode of inhibition, which will be the key for identification of potent inhibitors or antibodies demonstrating transmission-blocking activity., (© 2021 The Protein Society.)

Sequence logos give a fast and concise display in visualizing consensus sequence. Protein exhibits greater complexity and diversity than DNA, which usually affects the graphical representation of the logo. Reduced amino acids perform powerful ability for simplifying complexity of sequence alignment, which motivated us to establish RaacLogo. As a new sequence logo generator by using reduced amino acid alphabets, RaacLogo can easily generate many different simplified logos tailored to users by selecting various reduced amino acid alphabets that consisted of more than 40 clustering algorithms. This current web server provides 74 types of reduced amino acid alphabet, which were manually extracted to generate 673 reduced amino acid clusters (RAACs) for dealing with protein alignment. A two-dimensional selector was proposed for easily selecting desired RAACs with underlying biology knowledge. It is anticipated that the RaacLogo web server will play more high-potential roles for protein sequence alignment, topological estimation and protein design experiments. RaacLogo is freely available at http://bioinfor.imu.edu.cn/raaclogo., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder caused by deletion or subtle variant of survival motor neuron 1 (SMN1) gene. By multiplex ligation-dependent probe amplification, genomic sequencing, and T-A cloning on cDNA level, we identified one novel SMN1 subtle variant c.835G>C (p.Gly279Arg) in a non-homozygous patient with type 1 SMA. Full-length SMN1 (fl-SMN1) transcripts in the peripheral bloods of the patient were significantly decreased compared with those in healthy individuals and the carries (p < 0.05). And two fragments of SMN1 transcripts including fl-SMN1 and △7-SMN1 were observed by RT-PCR, which indicated Exon 7 skipping of SMN1 gene. To further evaluate its splicing effects on Exon 7, we performed ex vivo splicing analysis, which showed that the mutant mini gene with c.835G>C reduced Exon 7 inclusion to 54%. In addition, self-oligomerization between mutant SMN protein with the c.835G>C (p.Gly279Arg) and wild SMN was decreased in self-interaction assays. Our study clearly demonstrates that the c.835G>C (p.Gly279Arg) variant can lead to a decrease in fl-SMN1 transcripts by interrupting correct splicing of SMN1. What is more, the variant also affects SMN self-oligomerization via amino acid substitution from Gly to Arg at amino acid position of 279. This work presents the first evidence that it does exit double-hit events for the novel variant, which is crucial to understanding a severe SMA phenotype (type 1).

Enzymes have high catalytic efficiency and low environmental impact, and are therefore potentially useful tools for various industrial processes. Crucially, however, natural enzymes do not always have the properties required for specific processes. It may be necessary, therefore, to design, engineer, and evolve enzymes with properties that are not found in natural enzymes. In particular, the creation of enzymes that are thermally stable and catalytically active at low temperature is desirable for processes involving both high and low temperatures. In the current study, we designed two ancestral sequences of 3-isopropylmalate dehydrogenase by an ancestral sequence reconstruction technique based on a phylogenetic analysis of extant homologous amino acid sequences. Genes encoding the designed sequences were artificially synthesized and expressed in Escherichia coli. The reconstructed enzymes were found to be slightly more thermally stable than the extant thermophilic homologue from Thermus thermophilus. Moreover, they had considerably higher low-temperature catalytic activity as compared with the T. thermophilus enzyme. Detailed analyses of their temperature-dependent specific activities and kinetic properties showed that the reconstructed enzymes have catalytic properties similar to those of mesophilic homologues. Collectively, our study demonstrates that ancestral sequence reconstruction can produce a thermally stable enzyme with catalytic properties adapted to low-temperature reactions.

We introduce a model of amino acid sequence evolution that accounts for the statistical behavior of real sequences induced by epistatic interactions. We base the model dynamics on parameters derived from multiple sequence alignments analyzed by using direct coupling analysis methodology. Known statistical properties such as overdispersion, heterotachy, and gamma-distributed rate-across-sites are shown to be emergent properties of this model while being consistent with neutral evolution theory, thereby unifying observations from previously disjointed evolutionary models of sequences. The relationship between site restriction and heterotachy is characterized by tracking the effective alphabet dynamics of sites. We also observe an evolutionary Stokes shift in the fitness of sequences that have undergone evolution under our simulation. By analyzing the structural information of some proteins, we corroborate that the strongest Stokes shifts derive from sites that physically interact in networks near biochemically important regions. Perspectives on the implementation of our model in the context of the molecular clock are discussed., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

We introduce a simple model that describes the average occurrence of point variations in a generic protein sequence. This model is based on the idea that mutations are more likely to be fixed at sites in contact with others that have mutated in the recent past. Therefore, we extend the usual assumptions made in protein coevolution by introducing a time dumping on the effect of a substitution on its surrounding and makes correlated substitutions happen in avalanches localized in space and time. The model correctly predicts the average correlation of substitutions as a function of their distance along the sequence. At the same time, it predicts an among-site distribution of the number of substitutions per site highly compatible with a negative binomial, consistently with experimental data. The promising outcomes achieved with this model encourage the application of the same ideas in the field of pairwise and multiple sequence alignment.

Protein target structures for the Critical Assessment of Structure Prediction round 13 (CASP13) were split into evaluation units (EUs) based on their structural domains, the domain organization of available templates, and the performance of servers on whole targets compared to split target domains. Eighty targets were split into 112 EUs. The EUs were classified into categories suitable for assessment of high accuracy modeling (or template-based modeling [TBM]) and topology (or free modeling [FM]) based on target difficulty. Assignment into assessment categories considered the following criteria: (a) the evolutionary relationship of target domains to existing fold space as defined by the Evolutionary Classification of Protein Domains (ECOD) database; (b) the clustering of target domains using eight objective sequence, structure, and performance measures; and (c) the placement of target domains in a scatter plot of target difficulty against server performance used in the previous CASP. Generally, target domains with good server predictions had close template homologs and were classified as TBM. Alternately, targets with poor server predictions represent a mixture of fast evolving homologs, structure analogs, and new folds, and were classified as FM or FM/TBM overlap., (© 2019 Wiley Periodicals, Inc.)

The rate at which new protein and gene sequences are being discovered has grown explosively in the omics era, which has increasingly complicated the efficient characterization and analysis of their biological properties. In this study, we propose a web-based graphical database tool, SeQuery, for intuitively visualizing proteome/genome networks by integrating the sequential, structural and functional information of sequences. As a demonstration of our tool's effectiveness, we constructed a graph database of G protein-coupled receptor (GPCR) sequences by integrating data from the UniProt, GPCRdb and RCSB PDB databases. Our tool attempts to achieve two goals: (i) given the sequence of a query protein, correctly and efficiently identify whether the protein is a GPCR, and, if so, define its sequential and functional roles in the GPCR superfamily; and (ii) present a panoramic view of the GPCR superfamily and its network centralities that allows users to explore the superfamily at various resolutions. Such a bottom-up-to-top-down view can provide the users with a comprehensive understanding of the GPCR superfamily through interactive navigation of the graph database. A test of SeQuery with the GPCR2841 dataset shows that it correctly identifies 99 out of 100 queried protein sequences. The developed tool is readily applicable to other biological networks, and we aim to expand SeQuery by including additional biological databases in the near future., (© The Author(s) 2019. Published by Oxford University Press.)

Residues in proteins that are in close spatial proximity are more prone to covariate as their interactions are likely to be preserved due to structural and evolutionary constraints. If we can detect and quantify such covariation, physical contacts may then be predicted in the structure of a protein solely from the sequences that decorate it. To carry out such predictions, and following the work of others, we have implemented a multivariate Gaussian model to analyze correlation in multiple sequence alignments. We have explored and tested several numerical encodings of amino acids within this model. We have shown that 1D encodings based on amino acid biochemical and biophysical properties, as well as higher dimensional encodings computed from the principal components of experimentally derived mutation/substitution matrices, do not perform as well as a simple twenty dimensional encoding with each amino acid represented with a vector of one along its own dimension and zero elsewhere. The optimum obtained from representations based on substitution matrices is reached by using 10 to 12 principal components; the corresponding performance is less than the performance obtained with the 20-dimensional binary encoding. We highlight also the importance of the prior when constructing the multivariate Gaussian model of a multiple sequence alignment.

African Horse Sickness Virus (AHSV) (Orbivirus genus, Reoviridae family) causes high mortality in naïve domestic horses with enormous economic and socio-emotional impact. There are nine AHSV serotypes showing limited cross neutralization. AHSV is transmitted by competent species of Culicoides biting midges. AHS is a serious threat beyond the African continent as endemic Culicoides species in moderate climates transmit the closely related prototype bluetongue virus. There is a desperate need for safe and efficacious vaccines, while DIVA (Differentiating Infected from Vaccinated) vaccines would accelerate control of AHS. Previously, we have shown that highly virulent AHSV with an in-frame deletion of 77 amino acids (aa) in NS3/NS3a is completely safe, does not cause viremia and shows protective capacity. This deletion mutant is a promising DISA (Disabled Infectious Single Animal) vaccine platform, since exchange of serotype specific virus proteins has been shown for all nine serotypes. Here, we show that a prototype NS3 competitive ELISA is DIVA compliant to AHS DISA vaccine platforms. Epitope mapping of NS3/NS3a shows that more research is needed to evaluate this prototype serological DIVA assay regarding sensitivity and specificity, in particular for AHSVs expressing antigenically different NS3/NS3a proteins. Further, an experimental panAHSV PCR test targeting genome segment 10 is developed that detects reference AHSV strains, whereas AHS DISA vaccine platforms were not detected. This DIVA PCR test completely guarantees genetic DIVA based on in silico and in vitro validation, although test validation regarding diagnostic sensitivity and specificity has not been performed yet. In conclusion, the prototype NS3 cELISA and the PCR test described here enable serological and genetic DIVA accompanying AHS DISA vaccine platforms., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

In phylogenomic profiling, the genomic context based methods are based on the observation that two or more proteins having the same pattern of presence or absence in many diverse genomes most likely have a functional link. In this research work, a tool (PATSIM) has been developed to predict the protein patterns based on the SOPM tool. In this tool, the secondary structure for CATH database protein sequences, predicted by the SOPM (Self Optimized Prediction Method) server is passed as input to fulfill objectives such as, (i) Predict the Amino Acid Pattern using the proposed Hybrid KMP and BM algorithm, (ii) Predict the physiochemical properties such as Hydrophobic Non-Polar ALKYL Amino Acid groups, Hydrophobic Non-Polar AROMATIC Amino Acid groups, Hydrophilic Polar Neutral Amino Acid groups, Hydrophilic Polar Acidic Amino Acid groups and Hydrophilic Polar Basic Amino Acid groups of protein sequence, (iii) Predict the secondary structure of protein where the structure of protein sequence is unknown, and (iv) Similarity analysis of protein sequence (structure unknown) with the CATH database. From the results, it is inferred that this tool effectively predicts the similarity between the sequences and also identifies the protein patterns for four secondary structural classes, namely Alpha Helix (h), Beta Sheet (e), Turn (t) and Coil (c). Based on the experimental results, it is inferred that this tool identifies the physiochemical properties of the protein sequence in an effective manner. The source code and its documentation for the PATSIM tool is freely available in the GitHub public repository (https://github.com/manimkn89/Protein-Sequence-Analysis)., (Copyright © 2018 Elsevier B.V. All rights reserved.)

In this study, we delineate an unsupervised clustering algorithm, minimum span clustering (MSC), and apply it to detect G-protein coupled receptor (GPCR) sequences and to study the GPCR network using a base dataset of 2770 GPCR and 652 non-GPCR sequences. High detection accuracy can be achieved with a proper dataset. The clustering results of GPCRs derived from MSC show a strong correlation between their sequences and functions. By comparing our level 1 MSC results with the GPCRdb classification, the consistency is 87.9% for the fourth level of GPCRdb, 89.2% for the third level, 98.4% for the second level, and 100% for the top level (the lowest resolution level of GPCRdb). The MSC results of GPCRs can be well explained by estimating the selective pressure of GPCRs, as exemplified by investigating the largest two subfamilies, peptide receptors (PRs) and olfactory receptors (ORs), in class A GPCRs. PRs are decomposed into three groups due to a positive selective pressure, whilst ORs remain as a single group due to a negative selective pressure. Finally, we construct and compare phylogenetic trees using distance-based and character-based methods, a combination of which could convey more comprehensive information about the evolution of GPCRs.

Bax∆2 is a functional pro-apoptotic Bax isoform having alterations in its N-terminus, but sharing the rest of its sequence with Baxα. Bax∆2 is unable to target mitochondria due to the loss of helix α1. Instead, it forms cytosolic aggregates and activates caspase 8. However, the functional domain(s) responsible for BaxΔ2 behavior have remained elusive. Here we show that disruption of helix α1 makes Baxα mimic the behavior of Bax∆2. However, the other alterations in the Bax∆2 N-terminus have no significant impact on aggregation or cell death. We found that the hallmark BH3 domain is necessary but not sufficient for aggregation-mediated cell death. We also noted that the core region shared by Baxα and Bax∆2 is required for the formation of large aggregates, which is essential for BaxΔ2 cytotoxicity. However, aggregation by itself is unable to trigger cell death without the C-terminus. Interestingly, the C-terminal helical conformation, not its primary sequence, appears to be critical for caspase 8 recruitment and activation. As Bax∆2 shares core and C-terminal sequences with most Bax isoforms, our results not only reveal a structural basis for Bax∆2-induced cell death, but also imply an intrinsic potential for aggregate-mediated caspase 8-dependent cell death in other Bax family members., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Mitochondrial endonuclease G from Leishmania infantum (LiEndoG) participates in the degradation of double-stranded DNA (dsDNA) during parasite cell death and is catalytically inactive at a pH of 8.0 or above. The presence, in the primary sequence, of an acidic amino acid-rich insertion exclusive to trypanosomatids and its spatial position in a homology-built model of LiEndoG led us to postulate that this peptide stretch might act as a pH sensor for self-inhibition. We found that a LiEndoG variant lacking residues 145-180 is indeed far more active than its wild-type counterpart at pH values >7.0. In addition, we discovered that (i) LiEndoG exists as a homodimer, (ii) replacement of Ser211 in the active-site SRGH motif with the canonical aspartate from the DRGH motif of other nucleases leads to a catalytically deficient enzyme, (iii) the activity of the S211D variant can be restored upon the concomitant replacement of Ala247 with Arg and (iv) a C-terminal extension is responsible for the observed preferential cleavage of single-stranded DNA (ssDNA) and ssDNA-dsDNA junctions. Taken together, our results support the view that LiEndoG is a multidomain molecular machine whose nuclease activity can be subtly modulated or even abrogated through architectural changes brought about by environmental conditions and interaction with other binding partners., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Receptor interactions with short linear peptide fragments (ligands) are at the base of many biological signaling processes. Conserved and information-rich amino acid patterns, commonly called sequence motifs, shape and regulate these interactions. Because of the properties of a receptor-ligand system or of the assay used to interrogate it, experimental data often contain multiple sequence motifs. GibbsCluster is a powerful tool for unsupervised motif discovery because it can simultaneously cluster and align peptide data. The GibbsCluster 2.0 presented here is an improved version incorporating insertion and deletions accounting for variations in motif length in the peptide input. In basic terms, the program takes as input a set of peptide sequences and clusters them into meaningful groups. It returns the optimal number of clusters it identified, together with the sequence alignment and sequence motif characterizing each cluster. Several parameters are available to customize cluster analysis, including adjustable penalties for small clusters and overlapping groups and a trash cluster to remove outliers. As an example application, we used the server to deconvolute multiple specificities in large-scale peptidome data generated by mass spectrometry. The server is available at http://www.cbs.dtu.dk/services/GibbsCluster-2.0., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Gonadotropin-releasing hormone (GnRH) is a key neuropeptide regulating reproduction in humans and other vertebrates. Recently, GnRH-like cDNAs and peptides were reported in marine mollusks, implying that GnRH-mediated reproduction is an ancient neuroendocrine system that arose prior to the divergence of protostomes and deuterostomes. Here, we evaluated the reproductive control system mediated by GnRH in the Pacific abalone Haliotis discus hannai. We cloned a prepro-GnRH cDNA (Hdh-GnRH) from the pleural-pedal ganglion (PPG) in H. discus hannai, and analyzed its spatiotemporal gene expression pattern. The open reading frame of Hdh-GnRH encodes a protein of 101 amino acids, consisting of a signal peptide, a GnRH dodecapeptide, a cleavage site, and a GnRH-associated peptide. This structure and sequence are highly similar to GnRH-like peptides reported for mollusks and other invertebrates. Quantitative polymerase chain reaction demonstrated that Hdh-GnRH mRNA was more strongly expressed in the ganglions (PPG and cerebral ganglion [CG]) than in other tissues (gonads, gills, intestine, hemocytes, muscle, and mantle) in both sexes. In females, the expression levels of Hdh-GnRH mRNA in the PPG and branchial ganglion (BG) were significantly higher at the ripe and partial spent stages than at the early and late active stages. In males, Hdh-GnRH mRNA levels in the BG showed a significant increase in the partial spent stage. Unexpectedly, Hdh-GnRH levels in the CG were not significantly different among the examined stages in both sexes. These results suggest that Hdh-GnRH mRNA expression profiles in the BG and possibly the PPG are tightly correlated with abalone reproductive activities., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Hypotheses about the functions of ancient proteins and the effects of historical mutations on them are often tested using ancestral protein reconstruction (APR)-phylogenetic inference of ancestral sequences followed by synthesis and experimental characterization. Usually, some sequence sites are ambiguously reconstructed, with two or more statistically plausible states. The extent to which the inferred functions and mutational effects are robust to uncertainty about the ancestral sequence has not been studied systematically. To address this issue, we reconstructed ancestral proteins in three domain families that have different functions, architectures, and degrees of uncertainty; we then experimentally characterized the functional robustness of these proteins when uncertainty was incorporated using several approaches, including sampling amino acid states from the posterior distribution at each site and incorporating the alternative amino acid state at every ambiguous site in the sequence into a single "worst plausible case" protein. In every case, qualitative conclusions about the ancestral proteins' functions and the effects of key historical mutations were robust to sequence uncertainty, with similar functions observed even when scores of alternate amino acids were incorporated. There was some variation in quantitative descriptors of function among plausible sequences, suggesting that experimentally characterizing robustness is particularly important when quantitative estimates of ancient biochemical parameters are desired. The worst plausible case method appears to provide an efficient strategy for characterizing the functional robustness of ancestral proteins to large amounts of sequence uncertainty. Sampling from the posterior distribution sometimes produced artifactually nonfunctional proteins for sequences reconstructed with substantial ambiguity., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Misfolding of the prion protein (PrP) is the key step in the transmission of spongiform pathologies in humans and several animals. Although PrP is highly conserved in mammals, a few changes in the sequence of endogenous PrP are proposed to confer protection to dogs, which were highly exposed to prion during the mad-cow epidemics. D159 is a unique amino acid found in PrP from dogs and other canines that was shown to alter surface charge, but its functional relevance has never been tested in vivo. Here, we show in transgenic Drosophila that introducing the N159D substitution on mouse PrP decreases its turnover. Additionally, mouse PrP-N159D demonstrates no toxicity and accumulates no pathogenic conformations, suggesting that a single D159 substitution is sufficient to prevent PrP conformational change and pathogenesis. Understanding the mechanisms mediating the protective activity of D159 is likely to lessen the burden of prion diseases in humans and domestic animals., (Copyright © 2016 Elsevier Inc. All rights reserved.)

GlycoSiteAlign is a tool designed to align amino acid sequences of variable length surrounding glycosylation sites depending on the knowledge of glycan structure. It is an exploratory resource intended for the identification of characteristic amino acid patterns of unique glycan-protein interactions. GlycoSiteAlign uses data from the UniCarbKB and UniProtKB databases, and it is hosted on ExPASy, the Swiss Institute of Bioinformatics resource portal. The user can select either specific or general glycan features, set the length of the protein fragments, and trigger an alignment with the option of including 90% homologous proteins. The tool previews and downloads alignments that may reveal amino acid patterns corresponding to selected features (e.g., "fucosylated" versus "non-fucosylated"). GlycoSiteAlign will integrate new data as they become available to confirm and expand results. It is presented as a promising tool for assessing and refining the knowledge about the constraints that link a particular glycan structure to a particular glycosite, and in the long term, this application could help improve prediction tools.

Structural properties such as solvent accessibility and contact number predict site-specific sequence variability in many proteins. However, the strength and significance of these structure-sequence relationships vary widely among different proteins, with absolute correlation strengths ranging from 0 to 0.8. In particular, two recent works have made contradictory observations. Yeh et al. (Mol. Biol. Evol. 31:135-139, 2014) found that both relative solvent accessibility (RSA) and weighted contact number (WCN) are good predictors of sitewise evolutionary rate in enzymes, with WCN clearly out-performing RSA. Shahmoradi et al. (J. Mol. Evol. 79:130-142, 2014) considered these same predictors (as well as others) in viral proteins and found much weaker correlations and no clear advantage of WCN over RSA. Because these two studies had substantial methodological differences, however, a direct comparison of their results is not possible. Here, we reanalyze the datasets of the two studies with one uniform analysis pipeline, and we find that many apparent discrepancies between the two analyses can be attributed to the extent of sequence divergence in individual alignments. Specifically, the alignments of the enzyme dataset are much more diverged than those of the virus dataset, and proteins with higher divergence exhibit, on average, stronger structure-sequence correlations. However, the highest structure-sequence correlations are observed at intermediate divergence levels, where both highly conserved and highly variable sites are present in the same alignment., (© 2016 The Protein Society.)

Continual advancements in computing power and sophistication, coupled with rapid increases in protein sequence and structural information, have made bioinformatic tools an invaluable resource for the molecular and structural biologist. With the degree of sequence information continuing to expand at an almost exponential rate, it is essential that scientists today have a basic understanding of how to utilise, manipulate and analyse this information for the benefit of their own experiments. In the context of Toll-Interleukin I Receptor domain containing proteins, we describe here a series of the more common and user-friendly bioinformatic tools available as Internet-based resources. These will enable the identification and alignment of protein sequences; the identification of functional motifs; the characterisation of protein secondary structure; the identification of protein structural folds and distantly homologous proteins; and the validation of the structural geometry of modelled protein structures.

Annotation on the reference genome of the C57BL6/J mouse has been an ongoing project ever since the draft genome was first published. Initially, the principle focus was on the identification of all protein-coding genes, although today the importance of describing long non-coding RNAs, small RNAs, and pseudogenes is recognized. Here, we describe the progress of the GENCODE mouse annotation project, which combines manual annotation from the HAVANA group with Ensembl computational annotation, alongside experimental and in silico validation pipelines from other members of the consortium. We discuss the more recent incorporation of next-generation sequencing datasets into this workflow, including the usage of mass-spectrometry data to potentially identify novel protein-coding genes. Finally, we will outline how the C57BL6/J genebuild can be used to gain insights into the variant sites that distinguish different mouse strains and species.

Two native peptides with disulfide-directed hairpin (DDH) fold, LaIT1 and LITX, were recently isolated from scorpion venom, a development that offered insights into exploring the evolutionary linkage between DDH and inhibitor cystine knot (ICK) peptides. In this work, we isolated and identified the full-length cDNAs of LaTI1, a representative member with DDH fold, and further determined its complete gene structure. The precursor organization of LaIT1 is similar to that of ICK peptides. The LaIT1 gene contains four exons interrupted by three unique introns and differed from ICK peptides, suggesting divergent genomic organizations of DDH peptides and ICK peptides. Phylogenetic analysis further showed that the "simple" DDH peptide originates from the "complex" ICK peptide, rather than the reverse. To the best of our knowledge, this is the first report on the genomic organization of DDH-fold peptides, and it presents new evidence of an evolutionary linkage between ICK and DDH peptides., (© 2015 Wiley Periodicals, Inc.)

Many protein interactions are conserved among organisms despite changes in the amino acid sequences that comprise their contact sites, a property that has been used to infer the location of these sites from protein homology. In an inter-species complementation experiment, a sequence present in a homologue is substituted into a protein and tested for its ability to support function. Therefore, substitutions that inhibit function can identify interaction sites that changed over evolution. However, most of the sequence differences within a protein family remain unexplored because of the small-scale nature of these complementation approaches. Here we use existing high throughput mutational data on the in vivo function of the RRM2 domain of the Saccharomyces cerevisiae poly(A)-binding protein, Pab1, to analyze its sites of interaction. Of 197 single amino acid differences in 52 Pab1 homologues, 17 reduce the function of Pab1 when substituted into the yeast protein. The majority of these deleterious mutations interfere with the binding of the RRM2 domain to eIF4G1 and eIF4G2, isoforms of a translation initiation factor. A large-scale mutational analysis of the RRM2 domain in a two-hybrid assay for eIF4G1 binding supports these findings and identifies peripheral residues that make a smaller contribution to eIF4G1 binding. Three single amino acid substitutions in yeast Pab1 corresponding to residues from the human orthologue are deleterious and eliminate binding to the yeast eIF4G isoforms. We create a triple mutant that carries these substitutions and other humanizing substitutions that collectively support a switch in binding specificity of RRM2 from the yeast eIF4G1 to its human orthologue. Finally, we map other deleterious substitutions in Pab1 to inter-domain (RRM2-RRM1) or protein-RNA (RRM2-poly(A)) interaction sites. Thus, the combined approach of large-scale mutational data and evolutionary conservation can be used to characterize interaction sites at single amino acid resolution.

Tumor necrosis factor (TNF) is a potent inflammatory cytokine produced during inflammation. In this study, two Crucian carp TNF-alpha genes, TNFalpha-1 and TNFalpha-2 were cloned and sequenced. The TNFalpha-1 was 720 bp and-consisted of a 699 bp opening reading frame encoding 232 amino acids and TNFalpha-2 was 793 bp and contained an open reading frame of 687 bp which encoded 228 amino acids. The genomic structure of both genes consists of4 exons and 3 introns similar to other known TNF-alpha genes. The amino acid sequence ofcrucian carp TNFalpha-1 and TNFalpha-2 shared the highest identity with common carp, goldfish and the lowest with flounder and human. The phylogenetic analysis grouped crucian carp TNFalpha-1 and TNFalpha-2 with other cypriniformes, which are more closely related to the goldfish and common carp TNFEa-1 and TNFalpha-2. Real time PCR analysis showed a con- stitutive expression of crucian carp TNFalpha-1 and TNFalpha-2 in all seven tissues examined. The TNFalpha-1 mRNA was ex- pressed significantly higher in the liver and kidney than those of TNFalpha-2 while TNFalpha-2 was expressed significantly higher in the muscle than that of TNFalpha-1. Aeromonas hydrophila BSK-10 strain upregulated the expression level of TNFalpha-1 and TNFalpha-2 in all tissues tested. At 6 h, the expression level of TNFalpha-1 increased significantly higher in muscle, skin and liver, while the expression level of TNFalpha-2 increased significantly higher in muscle and gill. TNFalpha-1 expressed much stronger than TNFalpha-2. At 12 h the expression level started declining and was more reduced at 24 h. These results imply that both TNFalpha-1 and TNFalpha-2 mRNA are distributed differently in tissues and are implicated in the immune response to bacterial infection.

Serine protease inhibitors have been widely discovered from different animal venoms, but most of them belong to Kunitz-type toxin subfamily. Here, by screening scorpion venom gland cDNA libraries, we identified four new non-Kunitz serine protease inhibitors with a conserved Ascaris-type structural fold: Ascaris-type toxins Lychas mucronatus Ascaris-type protease inhibitor (LmAPI), Pandinus cavimanus Ascaris-type protease inhibitor (PcAPI), Pandinus cavimanus Ascaris-type protease inhibitor 2 (PcAPI-2), and Hottentotta judaicus Ascaris-type protease inhibitor (HjAPI). The detailed characterization of one Ascaris-type toxin LmAPI was further carried out, which contains 60 residues and possesses a classical Ascaris-type cysteine framework reticulated by five disulfide bridges. Enzyme and inhibitor reaction kinetics experiments showed that recombinant LmAPI inhibits the activity of chymotrypsin potently with a Ki value of 15.5 nM, but has little effect on trypsin and elastase. Bioinformatics analyses suggested that LmAPI contains unique functional residues "TQD" and might be a useful template to produce specific protease inhibitors. Our results indicated that animal venoms are a natural source of new type of protease inhibitors, which will accelerate the development of diagnostic and therapeutic agents for human diseases that target diverse proteases., (Copyright © 2014 Elsevier B.V. All rights reserved.)

An inducible and aromatic nitrilase from Gordonia terrae was purified with a yield of 19%. The enzyme had turnover number of 63 s⁻¹ x 10⁻¹, Km 1.4 mM and Vmax 95 Umg⁻¹ protein for benzonitrile. The nitrilase of G. terrae was active at basic pH (7-10), moderate temperature (20-45 °C) and has a half-life of 4 h at 35 °C. MALDI analysis and amino acid sequence deduced from cloned nucleotide fragment showed 97% homology with putative amidohydrolase of Gordonia sputi NBRC 100414 and G. namibiensis. The enzyme showed regioselectivity towards hydroxybenzonitriles, as different position of hydroxyl group i.e. meta-, para- and orthosubstitutions on benzonitrile effect enzyme activity. The in-silico interactions of these substrates with the predicted 3D model of this enzyme also showed differential interaction between hydroxyl group of substrates and the polar amino acids surrounding enzyme's active site. This leads to different proximity and orientation of substrates vis-a-vis their interaction with catalytic residues.

Alignment-free methods, in which shared properties of sub-sequences (e.g. identity or match length) are extracted and used to compute a distance matrix, have recently been explored for phylogenetic inference. However, the scalability and robustness of these methods to key evolutionary processes remain to be investigated. Here, using simulated sequence sets of various sizes in both nucleotides and amino acids, we systematically assess the accuracy of phylogenetic inference using an alignment-free approach, based on D2 statistics, under different evolutionary scenarios. We find that compared to a multiple sequence alignment approach, D2 methods are more robust against among-site rate heterogeneity, compositional biases, genetic rearrangements and insertions/deletions, but are more sensitive to recent sequence divergence and sequence truncation. Across diverse empirical datasets, the alignment-free methods perform well for sequences sharing low divergence, at greater computation speed. Our findings provide strong evidence for the scalability and the potential use of alignment-free methods in large-scale phylogenomics.

Background: Proteins that share a high sequence homology while exhibiting drastically different 3D structures are investigated in this study. Recently, artificial proteins related to the sequences of the GA and IgG binding GB domains of human serum albumin have been designed. These artificial proteins, referred to as GA and GB, share 98% amino acid sequence identity but exhibit different 3D structures, namely, a 3α bundle versus a 4β + α structure. Discriminating between their 3D structures based on their amino acid sequences is a very difficult problem. In the present work, in addition to using bioinformatics techniques, an analysis based on inter-residue average distance statistics is used to address this problem., Results: It was hard to distinguish which structure a given sequence would take only with the results of ordinary analyses like BLAST and conservation analyses. However, in addition to these analyses, with the analysis based on the inter-residue average distance statistics and our sequence tendency analysis, we could infer which part would play an important role in its structural formation., Conclusions: The results suggest possible determinants of the different 3D structures for sequences with high sequence identity. The possibility of discriminating between the 3D structures based on the given sequences is also discussed.

Higher order taxonomic assignments (family level and above) in the phylum Actinobacteria are currently based only on 16S-rRNA gene sequence analyses. Additional molecular markers need to be identified to increase the number of reference points for defining actinobacterial families and other higher taxa. Furthermore, since most novel actinobacterial taxa are defined at the level of species and genera, it is necessary to define molecular signatures at the genus level to enhance the robustness of genus descriptions. The current use of chemotaxonomic markers to define genera could be improved by the identification of genus-specific molecular signatures. In this study, GyrB amino acid sequences for members of the family Streptosporangiaceae were analysed for molecular signatures. Phylogenetic analyses showed that the gyrB gene tree supported the composition of the currently recognised genera in this family. The catalytically important amino acids were identified in the GyrB sequences, as were the GHKL superfamily motifs. Examination of GyrB protein sequence alignments revealed that there are genus-specific sequences for most of the multi-species genera and genus-defining amino acid insertions for the genera Herbidospora and Microbispora. Furthermore, there are GyrB signature amino acids which distinguish the family Streptosporangiaceae from the family Nocardiopsaceae., (Copyright © 2013 Elsevier GmbH. All rights reserved.)

The estimate of the root-mean-square deviation (r.m.s.d.) in coordinates between the model and the target is an essential parameter for calibrating likelihood functions for molecular replacement (MR). Good estimates of the r.m.s.d. lead to good estimates of the variance term in the likelihood functions, which increases signal to noise and hence success rates in the MR search. Phaser has hitherto used an estimate of the r.m.s.d. that only depends on the sequence identity between the model and target and which was not optimized for the MR likelihood functions. Variance-refinement functionality was added to Phaser to enable determination of the effective r.m.s.d. that optimized the log-likelihood gain (LLG) for a correct MR solution. Variance refinement was subsequently performed on a database of over 21,000 MR problems that sampled a range of sequence identities, protein sizes and protein fold classes. Success was monitored using the translation-function Z-score (TFZ), where a TFZ of 8 or over for the top peak was found to be a reliable indicator that MR had succeeded for these cases with one molecule in the asymmetric unit. Good estimates of the r.m.s.d. are correlated with the sequence identity and the protein size. A new estimate of the r.m.s.d. that uses these two parameters in a function optimized to fit the mean of the refined variance is implemented in Phaser and improves MR outcomes. Perturbing the initial estimate of the r.m.s.d. from the mean of the distribution in steps of standard deviations of the distribution further increases MR success rates.

The p53 transactivation domain (p53TAD) is an intrinsically disordered protein (IDP) domain that undergoes coupled folding and binding when interacting with partner proteins like the E3 ligase, MDM2, and the 70 kDa subunit of replication protein A, RPA70. The secondary structure and dynamics of six closely related mammalian homologues of p53TAD were investigated using nuclear magnetic resonance (NMR) spectroscopy. Differences in both transient secondary structure and backbone dynamics were observed for the homologues. Many of these differences were localized to the binding sites for MDM2 and RPA70. The amount of transient helical secondary structure observed for the MDM2 binding site was lower for the dog and mouse homologues, compared with human, and the amount of transient helical secondary structure observed for the RPA70 binding site was higher for guinea pig and rabbit, compared with human. Differences in the amount of transient helical secondary structure observed for the MDM2 binding site were directly related to amino acid substitutions occurring on the solvent exposed side of the amphipathic helix that forms during the p53TAD/MDM2 interaction. Differences in the amount of transient helical secondary structure were not as easily explained for the RPA70 binding site because of its extensive sequence divergence. Clustering analysis shows that the divergence in the transient secondary structure of the p53TAD homologues exceeds the amino acid sequence divergence. In contrast, strong correlations were observed between the backbone dynamics of the homologues and the sequence identity matrix, suggesting that the dynamic behavior of IDPs is a conserved evolutionary feature., (Copyright © 2013 Wiley Periodicals, Inc.)

We have recently showed that the weighted contact number profiles (or the packing density profiles) of proteins are well correlated with those of the corresponding sequence conservation profiles. The results suggest that a protein structure may contain sufficient information about sequence conservation comparable to that derived from multiple homologous sequences. However, there are ambiguities concerning how to compute the packing density of the subunit of a protein complex. For the subunits of a complex, there are different ways to compute its packing density--one including the packing contributions of the other subunits and the other one excluding their contributions. Here we selected two sets of enzyme complexes. Set A contains complexes with the active sites comprising residues from multiple subunits, while set B contains those with the active sites residing on single subunits. In Set A, if the packing density profile of a subunit is computed considering the contributions of the other subunits of the complex, it will agree better with the sequence conservation profile. But in Set B the situations are reversed. The results may be due to the stronger functional and structural constraints on the evolution processes on the complexes of Set A than those of Set B to maintain the enzymatic functions of the complexes. The comparison of the packing density and the sequence conservation profiles may provide a simple yet potentially useful way to understanding the structural and evolutionary couplings between the subunits of protein complexes., (Copyright © 2013 Wiley Periodicals, Inc.)

Since its first release in 2001 as mainly a software package for phylogenetic analysis, data analysis for molecular biology and evolution (DAMBE) has gained many new functions that may be classified into six categories: 1) sequence retrieval, editing, manipulation, and conversion among more than 20 standard sequence formats including MEGA, NEXUS, PHYLIP, GenBank, and the new NeXML format for interoperability, 2) motif characterization and discovery functions such as position weight matrix and Gibbs sampler, 3) descriptive genomic analysis tools with improved versions of codon adaptation index, effective number of codons, protein isoelectric point profiling, RNA and protein secondary structure prediction and calculation of minimum folding energy, and genomic skew plots with optimized window size, 4) molecular phylogenetics including sequence alignment, testing substitution saturation, distance-based, maximum parsimony, and maximum-likelihood methods for tree reconstructions, testing the molecular clock hypothesis with either a phylogeny or with relative-rate tests, dating gene duplication and speciation events, choosing the best-fit substitution models, and estimating rate heterogeneity over sites, 5) phylogeny-based comparative methods for continuous and discrete variables, and 6) graphic functions including secondary structure display, optimized skew plot, hydrophobicity plot, and many other plots of amino acid properties along a protein sequence, tree display and drawing by dragging nodes to each other, and visual searching of the maximum parsimony tree. DAMBE features a graphic, user-friendly, and intuitive interface and is freely available from http://dambe.bio.uottawa.ca (last accessed April 16, 2013).

Microsporidia from the Encephalitozoonidae are obligate intracellular parasites with highly conserved and compacted nuclear genomes: they have few introns, short intergenic regions, and almost identical gene complements and chromosome arrangements. Comparative genomics of Encephalitozoon and microsporidia in general have focused largely on the genomic diversity between different species, and we know very little about the levels of genetic diversity within species. Polymorphism studies with Encephalitozoon are so far restricted to a small number of genes, and a few genetically distinct strains have been identified; most notably, three genotypes (ECI, ECII, and ECIII) of the model species E. cuniculi have been identified based on variable repeats in the rRNA internal transcribed spacer (ITS). To determine if E. cuniculi genotypes are genetically distinct lineages across the entire genome and at the same time to examine the question of intraspecies genetic diversity in microsporidia in general, we sequenced de novo genomes from each of the three genotypes and analyzed patterns of single nucleotide polymorphisms (SNPs) and insertions/deletions across the genomes. Although the strains have almost identical gene contents, they harbor large numbers of SNPs, including numerous nonsynonymous changes, indicating massive intraspecies variation within the Encephalitozoonidae. Based on this diversity, we conclude that the recognized genotypes are genetically distinct and propose new molecular markers for microsporidian genotyping.

The aim of the work is to develop the common method for estimating the pairwise alignment quality versus the evolution distance (a degree of homology) between the sequences being comparing, and the type of the alignment procedure. 3D alignments as well as any data on 3D protein structure are not used in the study. Based on the accepted protein sequences evolution model it is possible to estimate the capability of the concrete alignment algorithm to recover the genuine alignment. In this study a classical Needleman and Wunsch global alignment algorithm has been tested on a set of sequences from the Prefab database. Dependences of accuracy and confidence of a global alignment procedure were calculated as consequences of insertions/deletions and mutations shares.

Background: Development of sensitive sequence search procedures for the detection of distant relationships between proteins at superfamily/fold level is still a big challenge. The intermediate sequence search approach is the most frequently employed manner of identifying remote homologues effectively. In this study, examination of serine proteases of prolyl oligopeptidase, rhomboid and subtilisin protein families were carried out using plant serine proteases as queries from two genomes including A. thaliana and O. sativa and 13 other families of unrelated folds to identify the distant homologues which could not be obtained using PSI-BLAST., Methodology/principal Findings: We have proposed to start with multiple queries of classical serine protease members to identify remote homologues in families, using a rigorous approach like Cascade PSI-BLAST. We found that classical sequence based approaches, like PSI-BLAST, showed very low sequence coverage in identifying plant serine proteases. The algorithm was applied on enriched sequence database of homologous domains and we obtained overall average coverage of 88% at family, 77% at superfamily or fold level along with specificity of ~100% and Mathew's correlation coefficient of 0.91. Similar approach was also implemented on 13 other protein families representing every structural class in SCOP database. Further investigation with statistical tests, like jackknifing, helped us to better understand the influence of neighbouring protein families., Conclusions/significance: Our study suggests that employment of multiple queries of a family for the Cascade PSI-BLAST searches is useful for predicting distant relationships effectively even at superfamily level. We have proposed a generalized strategy to cover all the distant members of a particular family using multiple query sequences. Our findings reveal that prior selection of sequences as query and the presence of neighbouring families can be important for covering the search space effectively in minimal computational time. This study also provides an understanding of the 'bridging' role of related families.

Unique purine-rich mRNA sequences embedded in the coding sequences of a distinct group of gammaherpesvirus maintenance proteins underlie the ability of the latently infected cell to minimize immune recognition. The Epstein-Barr virus nuclear antigen, EBNA1, a well characterized lymphocryptovirus maintenance protein has been shown to inhibit in cis antigen presentation, due in part to a large internal repeat domain encoding glycine and alanine residues (GAr) encoded by a purine-rich mRNA sequence. Recent studies have suggested that it is the purine-rich mRNA sequence of this repeat region rather than the encoded GAr polypeptide that directly inhibits EBNA1 self-synthesis and contributes to immune evasion. To test this hypothesis, we generated a series of EBNA1 internal repeat frameshift constructs and assessed their effects on cis-translation and endogenous antigen presentation. Diverse peptide sequences resulting from alternative repeat reading frames did not alleviate the translational inhibition characteristic of EBNA1 self-synthesis or the ensuing reduced surface presentation of EBNA1-specific peptide-MHC class I complexes. Human cells expressing the EBNA1 frameshift variants were also poorly recognized by antigen-specific T-cells. Furthermore, a comparative analysis of the mRNA sequences of the corresponding repeat regions of different viral maintenance homologues highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. Based on these combined observations, we propose that the cis-translational inhibitory effect of the EBNA1 internal repeat sequence operates mechanistically at the nucleotide level, potentially through RNA secondary structural elements, and is unlikely to be mediated through the GAr polypeptide. The demonstration that the EBNA1 repeat mRNA sequence and not the encoded protein sequence underlies immune evasion in this class of virus suggests a novel approach to therapeutic development through the use of anti-sense strategies or small molecules targeting EBNA1 mRNA structure.

Industrial dyes such as azodyes are potential environmental pollutants causing deleterious health hazards complications. These dyes are potentially degraded by azoreductase enzyme which is widely distributed in bacteria and also cyanobacteria. The azoreductase enzymes from cyanobacteria have not been explored in detail. Hence this enzyme from Nostoc sp. PCC 7120 has been addressed in detail in the present study considering to explore the physico-chemical properties, evolutionary relationships, functional sites and structural properties using various bioinformatics tools. Four conserved regions were obtained from the multiple sequence analysis. The multiple sequence alignment showed conserved regions at different stretches from 1-11, 40-57, 82-120 and 161-177 amino acid residues. These regions could be used for designing degenerate primers or probes for PCR-based amplification or hybridization-based detection of azoreductase sequences from different source organisms. Domain analysis and functional site prediction showed the presence of functional sites and domain such as flavodoxin like fold responsible for enzyme activity. 3D model was constructed and the best model was selected and validated. Superimposition of the final structure and the template showed variations in certain regions which might be involved in the accommodation of various dyes. Our results may be helpful for further investigations like docking studies as well as in vivo and in vitro conditions although these predictions still need to be studied.

Pro-sequences were swapped in cis between keratinases from Bacillus licheniformis (Ker BL) and Bacillus pumilus (Ker BP) to construct Ker ProBP-BL and Ker ProBL-BP, respectively. Expression of these keratinases was carried out constitutively by E. coli HB101-pEZZ18 system. They were characterized with respect to their parent enzymes, Ker BL and Ker BP, respectively. Ker ProBP-BL became more thermostable with a t(1/2) of 45 min at 80°C contrary to Ker BL which was not stable beyond 60°C. Similarly, the activity of Ker ProBP-BL on keratin and casein substrate, i.e. K:C ratio increased to 1.2 in comparison to 0.1 for Ker BL. Hydrolysis of insulin B-chain revealed that the cleavage sites increased to six from four in case of Ker ProBP-BL in comparison to Ker BL. However, cleavage sites decreased from seven to four in case of Ker ProBL-BP in comparison to the parent keratinase, Ker BP. Likewise, Ker ProBL-BP revealed altered pH and temperature kinetics with optima at pH 10 and 60°C in comparison to Ker BP which had optima at pH 9 and 70°C. It also cleaved soluble substrates with better efficiency in comparison to Ker BP with K:C ratio of 1.6. Pro-sequence mediated conformational changes were also observed in trans and were almost similar to the features acquired by the chimeras constructed in cis by swapping the pro-sequence region., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Based on the order of 6-bit binary Gray code, a cyclic order of 20 amino acids is introduced. A novel 3D graphical representation of protein sequences is proposed according to the CGR of DNA sequences. Furthermore, the mathematical descriptor is suggested to characterize the graphical representation curve. The efficiency of our approach can be illustrated by performing the comparison of similarities/dissimilarities among sequences of the ND5 proteins of nine different species. With the correlation and significance analysis, the comparisons of both our results and results of other graphical representation with the ClustalW's results can show the utility of our approach., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Score-based pairwise alignments are widely used in bioinformatics in particular with molecular database search tools, such as the BLAST family. Due to sophisticated heuristics, such algorithms are usually fast but the underlying scoring model unfortunately lacks a statistical description of the reliability of the reported alignments. In particular, close to gaps, in low-score or low-complexity regions, a huge number of alternative alignments arise which results in a decrease of the certainty of the alignment. ppALIGN is a software package that uses hidden Markov Model techniques to compute position-wise reliability of score-based pairwise alignments of DNA or protein sequences. The design of the model allows for a direct connection between the scoring function and the parameters of the probabilistic model. For this reason it is suitable to analyze the outcomes of popular score based aligners and search tools without having to choose a complicated set of parameters. By contrast, our program only requires the classical score parameters (the scoring function and gap costs). The package comes along with a library written in C++, a standalone program for user defined alignments (ppALIGN) and another program (ppBLAST) which can process a complete result set of BLAST. The main algorithms essentially exhibit a linear time complexity (in the alignment lengths), and they are hence suitable for on-line computations. We have also included alternative decoding algorithms to provide alternative alignments. ppALIGN is a fast program/library that helps detect and quantify questionable regions in pairwise alignments. Due to its structure, the input/output interface it can to be connected to other post-processing tools. Empirically, we illustrate its usefulness in terms of correctly predicted reliable regions for sequences generated using the ROSE model for sequence evolution, and identify sensor-specific regions in the denitrifying betaproteobacterium Aromatoleum aromaticum.

The TGA transcription factors belong to the subfamily of bZIP group D that play a major role in disease resistance and development. Most of the TGA identified in Arabidopsis interact with the master regulator of SAR, NPR1 that controls the expression of PR genes. As a first approach to determine the possible involvement of these transcription factors in papaya defense, we characterized Arabidopsis TGA orthologs from the genome of Carica papaya cv. SunUp. Six orthologs CpTGA1 to CpTGA6, were identified. The predicted CpTGA proteins were highly similar to AtTGA sequences and probably share the same DNA binding properties and transcriptional regulation features. The protein sequences alignment evidenced the presence of conserved domains, characteristic of this group of transcription factors. The phylogeny showed that CpTGA evolved into three different subclades associated with defense and floral development. This is the first report of basal expression patterns assessed by RT-PCR, from the whole subfamily of CpTGA members in different tissues from papaya cv. Maradol mature plants. Overall, CpTGA1, CpTGA3 CpTGA6 and CpTGA4 showed a basal expression in all tissues tested; CpTGA2 expressed strongly in all tissues except in petioles while CpTGA5 expressed only in petals and to a lower extent in petioles. Although more detailed studies in anthers and other floral structures are required, we suggest that CpTGA5 might be tissue-specific, and it might be involved in papaya floral development. On the other hand, we report here for the first time, the expression of the whole family of CpTGA in response to salicylic acid (SA). The expression of CpTGA3, CpTGA4 and CpTGA6 increased in response to SA, what would suggest its involvement in the SAR response in papaya., (Copyright © 2012 Elsevier Masson SAS. All rights reserved.)

The phosphoprotein scaffold Dishevelled is an essential component of both Wnt signalling and of the signalsome that constitutes the supermolecular 'punctae' of assembled proteins often observed in fluorescence microscopy. The C-terminal region beyond the DEP domain displays unique and interesting character, exploited herein by careful analysis of the primary structure. Human Dishevelled-1, -2, -3 and fly Dishevelled (Dsh) sequences were downloaded and interrogated in silico. The C-terminus of Dishevelled-3 is revealed by FoldIndex(®) to be rich in ordered structure. It displays primary sequence that is unique and divergent in important ways from vertebrate isoforms as well as from the fly Dsh. The region is amphipathic, high in prolyl content, and harbours polyprolines. Dishevelled-3 displays some regions, where the proline content is >40%. Polyprolyl sequences (2-4 residues) likely constitute important sites of interaction with other Dishevelled isoforms. Several histidine-single amino acid repeats are notable. The 637,638/647,648 repeats of Dvl3 are essential for Wnt non-canonical, but not canonical signalling. Mutagenesis reveals that the C-terminal sequence is essential for the formation of punctae, made visible by fluorescence microscopy. These Dvl3-based signalsomes are very large (25-35 MDa-MW), supermolecular complexes that display dynamic reorganization in response to Wnt stimulation. Dishevelled-3 C-terminus is rich in structure and unique motifs, worthy of detailed analysis with modern molecular tools., (© 2011 The Authors. Acta Physiologica © 2011 Scandinavian Physiological Society.)

Background: Protein sequences are subject to a mosaic of constraint. Changes to functional domains and buried residues, for example, are more apt to disrupt protein structure and function than are changes to residues participating in loops or exposed to solvent. Regions of constraint on the tertiary structure of a protein often result in loose segmentation of its primary structure into stretches of slowly- and rapidly-evolving amino acids. This clustering can be exploited, and existing methods have done so by relying on local sequence conservation as a signature of selection to help identify functionally important regions within proteins. We invert this paradigm by leveraging the regional nature of protein structure and function to both illuminate and make use of genome-wide patterns of local sequence conservation., Results: Our hypothesis is that the regional nature of structural and functional constraints will assert a positive autocorrelation on the evolutionary rates of neighboring sites, which, in a pairwise comparison of orthologous proteins, will manifest itself as the clustering of non-synonymous changes across the amino acid sequence. We introduce a dispersion ratio statistic to test this and related hypotheses. Using genome-wide interspecific comparisons of orthologous protein pairs, we reveal a strong log-linear relationship between the degree of clustering and the intensity of constraint. We further demonstrate how this relationship varies with the evolutionary distance between the species being compared. We provide some evidence that proteins with a history of positive selection deviate from genome-wide trends., Conclusions: We find a significant association between the evolutionary rate of a protein and the degree to which non-synonymous changes cluster along its primary sequence. We show that clustering is a non-redundant predictor of evolutionary rate, and we speculate that conflicting signals of clustering and constraint may be indicative of a historical period of relaxed selection.