HMM Logos for visualization of protein families

BACKGROUND: Profile Hidden Markov Models (pHMMs) are a widely used tool for protein family research. Up to now, however, there exists no method to visualize all of their central aspects graphically in an intuitively understandable way. RESULTS: We present a visualization method that incorporates both emission and transition probabilities of the pHMM, thus extending sequence logos introduced by Schneider and Stephens. For each emitting state of the pHMM, we display a stack of letters. The stack height is determined by the deviation of the position's letter emission frequencies from the background frequencies. The stack width visualizes both the probability of reaching the state (the hitting probability) and the expected number of letters the state emits during a pass through the model (the state's expected contribution). A web interface offering online creation of HMM Logos and the corresponding source code can be found at the Logos web server of the Max Planck Institute for Molecular Genetics . CONCLUSIONS: We demonstrate that HMM Logos can be a useful tool for the biologist: We use them to highlight differences between two homologous subfamilies of GTPases, Rab and Ras, and we show that they are able to indicate structural elements of Ras

webPRC: the Profile Comparer for alignment-based searching of public domain databases

Profile–profile methods are well suited to detect remote evolutionary relationships between protein families. Profile Comparer (PRC) is an existing stand-alone program for scoring and aligning hidden Markov models (HMMs), which are based on multiple sequence alignments. Since PRC compares profile HMMs instead of sequences, it can be used to find distant homologues. For this purpose, PRC is used by, for example, the CATH and Pfam-domain databases. As PRC is a profile comparer, it only reports profile HMM alignments and does not produce multiple sequence alignments. We have developed webPRC server, which makes it straightforward to search for distant homologues or similar alignments in a number of domain databases. In addition, it provides the results both as multiple sequence alignments and aligned HMMs. Furthermore, the user can view the domain annotation, evaluate the PRC hits with the Jalview multiple alignment editor and generate logos from the aligned HMMs or the aligned multiple alignments. Thus, this server assists in detecting distant homologues with PRC as well as in evaluating and using the results. The webPRC interface is available at http://www.ibi.vu.nl/programs/prcwww/

MTERF3, the most conserved member of the mTERF-family, is a modular factor involved in mitochondrial protein synthesis

AbstractThe MTERF-family is a wide family of proteins identified in Metazoa and plants which includes the known mitochondrial transcription termination factors. With the aim to shed light on the function of MTERF-family members in Drosophila, we performed the cloning and characterization of D-MTERF3, a component of the most conserved group of this family. D-MTERF3 is a mitochondrial protein of 323 amino acids. Sequence analysis in seven different organisms showed that the protein contains five conserved “mTERF-motifs”, three of which include a leucine zipper-like domain. D-MTERF3 knock-down, obtained by RNAi in D.Mel-2 cells, did not affect mitochondrial replication and transcription. On the contrary, it decreased to a variable extent the rate of labelling of about half of the mitochondrial polypeptides, with ND1 being the most affected by D-MTERF3 depletion. These results indicate that D-MTERF3 is involved in mitochondrial translation. This role, likely based on protein–protein interactions, may be exerted either through a direct interaction with the translation machinery or by bridging the mitochondrial transcription and translation apparatus

Zygomycetes, Microsporidia, and the Evolutionary Ancestry of Sex Determination

Zygomycetes and their alleged sister taxon, the microsporidia, exclusively share the presence of a cluster of three genes encoding a sugar transporter, a high mobility group (HMG)-type transcription factor, and an RNA helicase. In zygomycetes, the HMG-type transcription factor acts as the sole sex determinant. This intimately ties the evolutionary history of this gene cluster to the evolution of sex determination. Here, we have unraveled the relationships of the two gene clusters by vicariously analyzing the sugar transporters and the RNA helicases. We show that if the two gene clusters share a common ancestry, it dates back to the early days of eukaryotic evolution. As a consequence, the zygomycete MAT locus would be old enough to represent the archetype of fungal and animal sex determination. However, the evolutionary scenario that has to be invoked is complex. An independent assembly of the two clusters deserves therefore consideration. In either case, shared ancestry or convergent evolution, the presence of the gene cluster in microsporidia and in zygomycetes represents at best a plesiomorphy. Hence, it is not phylogenetically informative. A further genome-wide reanalysis of gene order conservation reveals that gene order is not significantly more similar between microsporidia and zygomycetes than between microsporidia and any other fungal taxon or even humans. Consequently, the phylogenetic placement of microsporidia as sister to the zygomycetes needs to be reconsidered

Comparative analysis of carboxysome shell proteins

Carboxysomes are metabolic modules for CO2 fixation that are found in all cyanobacteria and some chemoautotrophic bacteria. They comprise a semi-permeable proteinaceous shell that encapsulates ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) and carbonic anhydrase. Structural studies are revealing the integral role of the shell protein paralogs to carboxysome form and function. The shell proteins are composed of two domain classes: those with the bacterial microcompartment (BMC; Pfam00936) domain, which oligomerize to form (pseudo)hexamers, and those with the CcmL/EutN (Pfam03319) domain which form pentamers in carboxysomes. These two shell protein types are proposed to be the basis for the carboxysome’s icosahedral geometry. The shell proteins are also thought to allow the flux of metabolites across the shell through the presence of the small pore formed by their hexameric/pentameric symmetry axes. In this review, we describe bioinformatic and structural analyses that highlight the important primary, tertiary, and quaternary structural features of these conserved shell subunits. In the future, further understanding of these molecular building blocks may provide the basis for enhancing CO2 fixation in other organisms or creating novel biological nanostructures

GPS-ARM: Computational Analysis of the APC/C Recognition Motif by Predicting D-Boxes and KEN-Boxes

Anaphase-promoting complex/cyclosome (APC/C), an E3 ubiquitin ligase incorporated with Cdh1 and/or Cdc20 recognizes and interacts with specific substrates, and faithfully orchestrates the proper cell cycle events by targeting proteins for proteasomal degradation. Experimental identification of APC/C substrates is largely dependent on the discovery of APC/C recognition motifs, e.g., the D-box and KEN-box. Although a number of either stringent or loosely defined motifs proposed, these motif patterns are only of limited use due to their insufficient powers of prediction. We report the development of a novel GPS-ARM software package which is useful for the prediction of D-boxes and KEN-boxes in proteins. Using experimentally identified D-boxes and KEN-boxes as the training data sets, a previously developed GPS (Group-based Prediction System) algorithm was adopted. By extensive evaluation and comparison, the GPS-ARM performance was found to be much better than the one using simple motifs. With this powerful tool, we predicted 4,841 potential D-boxes in 3,832 proteins and 1,632 potential KEN-boxes in 1,403 proteins from H. sapiens, while further statistical analysis suggested that both the D-box and KEN-box proteins are involved in a broad spectrum of biological processes beyond the cell cycle. In addition, with the co-localization information, we predicted hundreds of mitosis-specific APC/C substrates with high confidence. As the first computational tool for the prediction of APC/C-mediated degradation, GPS-ARM is a useful tool for information to be used in further experimental investigations. The GPS-ARM is freely accessible for academic researchers at: http://arm.biocuckoo.org

Insights from the Structure of Mycobacterium tuberculosis Topoisomerase I with a novel protein fold

The DNA topoisomerase I enzyme of Mycobacterium tuberculosis (MtTOP1) is essential for the viability of the organism and survival in a murine model. This topoisomerase is being pursued as a novel target for the discovery of new therapeutic agents for the treatment of drug-resistant tuberculosis. In this study, we succeeded in obtaining a structure of MtTOP1 by first predicting that the C-terminal region of MtTOP1 contains four repeated domains that do not involve the Zn-binding tetracysteine motifs seen in the C-terminal domains of Escherichiacoli topoisomerase I. A construct (amino acids A2–T704), MtTOP1-704t, that includes the N-terminal domains (D1–D4) and the first predicted C-terminal domain (D5) of MtTOP1 was expressed and found to retain DNA cleavage–religation activity and catalyze single-stranded DNA catenation. MtTOP1-704t was crystallized, and a structure of 2.52 Å resolution limit was obtained. The structure of the MtTOP1 N-terminal domains has features that have not been observed in other previously available bacterial topoisomerase I crystal structures. The first C-terminal domain D5 forms a novel protein fold of a four-stranded antiparallel β-sheet stabilized by a crossing-over α-helix. Since there is only one type IA topoisomerase present in Mycobacteriaceae and related Actinobacteria, this subfamily of type IA topoisomerase may be required for multiple functions in DNA replication, transcription, recombination, and repair. The unique structural features observed for MtTOP1 may allow these topoisomerase I enzymes to carry out physiological functions associated with topoisomerase III enzyme in other bacteria

Detecting species-site dependencies in large multiple sequence alignments

Multiple sequence alignments (MSAs) are one of the most important sources of information in sequence analysis. Many methods have been proposed to detect, extract and visualize their most significant properties. To the same extent that site-specific methods like sequence logos successfully visualize site conservations and sequence-based methods like clustering approaches detect relationships between sequences, both types of methods fail at revealing informational elements of MSAs at the level of sequence–site interactions, i.e. finding clusters of sequences and sites responsible for their clustering, which together account for a high fraction of the overall information of the MSA. To fill this gap, we present here a method that combines the Fisher score-based embedding of sequences from a profile hidden Markov model (pHMM) with correspondence analysis. This method is capable of detecting and visualizing group-specific or conflicting signals in an MSA and allows for a detailed explorative investigation of alignments of any size tractable by pHMMs. Applications of our methods are exemplified on an alignment of the Neisseria surface antigen LP2086, where it is used to detect sites of recombinatory horizontal gene transfer and on the vitamin K epoxide reductase family to distinguish between evolutionary and functional signals

High-throughput screening for substrate specificity-adapted mutants of the nisin dehydratase NisB

Microbial lanthipeptides are formed by a two-step enzymatic introduction of (methyl)lanthionine rings. A dehydratase catalyzes the dehydration of serine and threonine residues, yielding dehydroalanine and dehydrobutyrine, respectively. Cyclase-catalyzed coupling of the formed dehydroresidues to cysteines forms (methyl)lanthionine rings in a peptide. Lanthipeptide biosynthetic systems allow discovery of target-specific, lanthionine-stabilized therapeutic peptides. However, the substrate specificity of existing modification enzymes impose limitations on installing lanthionines in non-natural substrates. The goal of the present study was to obtain a lanthipeptide dehydratase with the capacity to dehydrate substrates that are unsuitable for the nisin dehydratase NisB. We report high-throughput screening for tailored specificity of intracellular, genetically encoded NisB dehydratases. The principle is based on the screening of bacterially displayed lanthionine-constrained streptavidin ligands, which have a much higher affinity for streptavidin than linear ligands. The designed NisC-cyclizable high-affinity ligands can be formed via mutant NisB-catalyzed dehydration but less effectively via wild-type NisB activity. In Lactococcus lactis, a cell surface display precursor was designed comprising DSHPQFC. The Asp residue preceding the serine in this sequence disfavors its dehydration by wild-type NisB. The cell surface display vector was coexpressed with a mutant NisB library and NisTC. Subsequently, mutant NisB-containing bacteria that display cyclized strep ligands on the cell surface were selected via panning rounds with streptavidin-coupled magnetic beads. In this way, a NisB variant with a tailored capacity of dehydration was obtained, which was further evaluated with respect to its capacity to dehydrate nisin mutants. These results demonstrate a powerful method for selecting lanthipeptide modification enzymes with adapted substrate specificity