Extracting protein alignment models from the sequence database

RNA-hydrolyzing activity of metallo-β-lactamase IMP-1

Metallo-β-lactamases (MBLs) hydrolyze a wide range of β-lactam antibiotics. While all MBLs share a common αβ/βα-fold, there are many other proteins with the same folding pattern that exhibit different enzymatic activities. These enzymes, together with MBLs, form the MBL superfamily. Thermotoga maritima tRNase Z, a tRNA 3′ processing endoribonuclease of MBL-superfamily, and IMP-1, a clinically isolated MBL, showed a striking similarity in tertiary structure, despite low sequence homology. IMP-1 hydrolyzed both total cellular RNA and synthetic small unstructured RNAs. IMP-1 also hydrolyzed pre-tRNA, but its cleavage site was different from those of T. maritima tRNase Z and human tRNase Z long form, indicating a key difference in substrate recognition. Single-turnover kinetic assays suggested that substrate-binding affinity of T. maritima tRNase Z is much higher than that of IMP-1.

M42 aminopeptidase catalytic site: the structural and functional role of a strictly conserved aspartate residue

The M42 aminopeptidases are a family of dinuclear aminopeptidases widely distributed in Prokaryotes. They are potentially associated to the proteasome, achieving complete peptide destruction. Their most peculiar characteristic is their quaternary structure, a tetrahedron-shaped particle made of twelve subunits. The catalytic site of M42 aminopeptidases is defined by seven conserved residues. Five of them are involved in metal ion binding which is important to maintain both the activity and the oligomeric state. The sixth conserved residue, a glutamate, is the catalytic base deprotonating the water molecule during peptide bond hydrolysis. The seventh residue is an aspartate whose function remains poorly understood. This aspartate residue, however, must have a critical role as it is strictly conserved in all MH clan enzymes. It forms some kind of catalytic triad with the histidine residue and the metal ion of the M2 binding site. We assess its role in TmPep1050, an M42 aminopeptidase of Thermotoga maritima, through a mutational approach. Asp-62 was substituted with alanine, asparagine, or glutamate residue. The Asp-62 substitutions completely abolished TmPep1050 activity and impeded dodecamer formation. They also interfered with metal ion binding as only one cobalt ion is bound per subunit instead of two. The structure of Asp62Ala variant was solved at 1.5 Å showing how the substitution has an impact on the active site fold. We propose a structural role for Asp-62, helping to stabilize a crucial loop in the active site and to position correctly the catalytic base and a metal ion ligand of the M1 site.

Molecular modeling and QM/MM calculation clarify the catalytic mechanism of β-lactamase N1

The treatment of bacterial infections is currently threatened by the emergence of pathogenic bacteria producing β-lactamase, which catalyzes the hydrolysis of β-lactams. Although the hydrolysis of the substrate nitrocefin by a metallo-β-lactamase, namely β-lactamase N1 from USA300 (a typical methicillin-resistant Staphylococcus aureus), has previously been reported in the literature, its mechanism remains elusive. Here, we show that molecular modeling and quantum-mechanical/molecular mechanics (QM/MM) calculations describing the complex of β-lactamase N1 with nitrocefin (the substrate of β-lactamase N1) can predict the catalytic mechanism of nitrocefin hydrolysis by β-lactamase N1. Molecular dynamics simulation shows that the catalytic reaction begins with hydrogen bond formation between Gln171 and a water molecule, which is thereby captured for nitrocefin hydrolysis by β-lactamase N1. In addition, the carboxyl group coordinates Zn2 in a chelating fashion. The binding energy decompositions suggest that Phe169 anchors nitrocefin by π-stacking interactions between the benzene rings. Specifically, Phe169 and Zn2 position the nitrocefin in specific orientations. The active site of β-lactamase N1 contains two residues (Gln171 and Phe169) that we expected to be crucial for guiding the nitrocefin hydrolysis reaction. Compelling evidence is provided that the mutants F169A and Q171A show lower enzymatic activity than the wild-type protein. On the basis of the QM/MM calculations, we propose that nitrocefin hydrolysis is initiated by the interaction between the oxygen atom of water and the C18 atom of nitrocefin, leading to the opening of the four-membered ring of nitrocefin and the formation of a substrate intermediate. In the next step, a hydrogen atom transfers from the nitrogen atom to the C11 atom of nitrocefin, resulting in the stable product.

Inferring joint sequence-structural determinants of protein functional specificity

Residues responsible for allostery, cooperativity, and other subtle but functionally important interactions remain difficult to detect. To aid such detection, we employ statistical inference based on the assumption that residues distinguishing a protein subgroup from evolutionarily divergent subgroups often constitute an interacting functional network. We identify such networks with the aid of two measures of statistical significance. One measure aids identification of divergent subgroups based on distinguishing residue patterns. For each subgroup, a second measure identifies structural interactions involving pattern residues. Such interactions are derived either from atomic coordinates or from Direct Coupling Analysis scores, used as surrogates for structural distances. Applying this approach to N-acetyltransferases, P-loop GTPases, RNA helicases, synaptojanin-superfamily phosphatases and nucleases, and thymine/uracil DNA glycosylases yielded results congruent with biochemical understanding of these proteins, and also revealed striking sequence-structural features overlooked by other methods. These and similar analyses can aid the design of drugs targeting allosteric sites.

Mammalian Fe-S proteins: definition of a consensus motif recognized by the co-chaperone HSC20

Iron-sulfur (Fe-S) clusters are inorganic cofactors that are fundamental to several biological processes in all three kingdoms of life. In most organisms, Fe-S clusters are initially assembled on a scaffold protein, ISCU, and subsequently transferred to target proteins or to intermediate carriers by a dedicated chaperone/co-chaperone system. The delivery of assembled Fe-S clusters to recipient proteins is a crucial step in the biogenesis of Fe-S proteins, and, in mammals, it relies on the activity of a multiprotein transfer complex that contains the chaperone HSPA9, the co-chaperone HSC20 and the scaffold ISCU. How the transfer complex efficiently engages recipient Fe-S target proteins involves specific protein interactions that are not fully understood. This mini review focuses on recent insights into the molecular mechanism of amino acid motif recognition and discrimination by the co-chaperone HSC20, which guides Fe-S cluster delivery.

The Role of Active Site Flexible Loops in Catalysis and of Zinc in Conformational Stability of Bacillus cereus 569/H/9 β-Lactamase

Metallo-β-lactamases catalyze the hydrolysis of most β-lactam antibiotics and hence represent a major clinical concern. The development of inhibitors for these enzymes is complicated by the diversity and flexibility of their substrate-binding sites, motivating research into their structure and function. In this study, we examined the conformational properties of the Bacillus cereus β-lactamase II in the presence of chemical denaturants using a variety of biochemical and biophysical techniques. The apoenzyme was found to unfold cooperatively, with a Gibbs free energy of stabilization (ΔG(0)) of 32 ± 2 kJ·mol(-1) For holoBcII, a first non-cooperative transition leads to multiple interconverting native-like states, in which both zinc atoms remain bound in an apparently unaltered active site, and the protein displays a well organized compact hydrophobic core with structural changes confined to the enzyme surface, but with no catalytic activity. Two-dimensional NMR data revealed that the loss of activity occurs concomitantly with perturbations in two loops that border the enzyme active site. A second cooperative transition, corresponding to global unfolding, is observed at higher denaturant concentrations, with ΔG(0) value of 65 ± 1.4 kJ·mol(-1) These combined data highlight the importance of the two zinc ions in maintaining structure as well as a relatively well defined conformation for both active site loops to maintain enzymatic activity.

Bayesian Top-Down Protein Sequence Alignment with Inferred Position-Specific Gap Penalties

We describe a Bayesian Markov chain Monte Carlo (MCMC) sampler for protein multiple sequence alignment (MSA) that, as implemented in the program GISMO and applied to large numbers of diverse sequences, is more accurate than the popular MSA programs MUSCLE, MAFFT, Clustal-Ω and Kalign. Features of GISMO central to its performance are: (i) It employs a "top-down" strategy with a favorable asymptotic time complexity that first identifies regions generally shared by all the input sequences, and then realigns closely related subgroups in tandem. (ii) It infers position-specific gap penalties that favor insertions or deletions (indels) within each sequence at alignment positions in which indels are invoked in other sequences. This favors the placement of insertions between conserved blocks, which can be understood as making up the proteins' structural core. (iii) It uses a Bayesian statistical measure of alignment quality based on the minimum description length principle and on Dirichlet mixture priors. Consequently, GISMO aligns sequence regions only when statistically justified. This is unlike methods based on the ad hoc, but widely used, sum-of-the-pairs scoring system, which will align random sequences. (iv) It defines a system for exploring alignment space that provides natural avenues for further experimentation through the development of new sampling strategies for more efficiently escaping from suboptimal traps. GISMO's superior performance is illustrated using 408 protein sets containing, on average, 235 sequences. These sets correspond to NCBI Conserved Domain Database alignments, which have been manually curated in the light of available crystal structures, and thus provide a means to assess alignment accuracy. GISMO fills a different niche than other MSA programs, namely identifying and aligning a conserved domain present within a large, diverse set of full length sequences. The GISMO program is available at http://gismo.igs.umaryland.edu/.

Gene Network Analysis of Metallo Beta Lactamase Family Proteins Indicates the Role of Gene Partners in Antibiotic Resistance and Reveals Important Drug Targets

Metallo Beta (β) Lactamases (MBL) are metal dependent bacterial enzymes that hydrolyze the β-lactam antibiotics. In recent years, MBL have received considerable attention because it inactivates most of the β-lactam antibiotics. Increase in dissemination of MBL encoding antibiotic resistance genes in pathogenic bacteria often results in unsuccessful treatments. Gene interaction network of MBL provides a complete understanding on the molecular basis of MBL mediated antibiotic resistance. In our present study, we have constructed the MBL network of 37 proteins with 751 functional partners from pathogenic bacterial spp. We found 12 highly interconnecting clusters. Among the 37 MBL proteins considered in the present study, 22 MBL proteins are from B3 subclass, 14 are from B1 subclass and only one is from B2 subclass. Global topological parameters are used to calculate and compare the probability of interactions in MBL proteins. Our results indicate that the proteins associated within the network have a strong influence in antibiotic resistance mechanism. Interestingly, several drug targets are identified from the constructed network. We believe that our results would be helpful for researchers exploring MBL-mediated antibiotic resistant mechanisms.

Structural and functional studies of a metallo-β-lactamase unveil a new type of structurally encoded nickel-containing heterodinuclear site

The selection of correct metal ions with high fidelity against competing cellular cations is crucial for the function of many metalloenzymes; however, the understanding of the principles that govern metal selectivity is still incomplete. In this study, the crystal structure of the Tm1162 protein from Thermotoga maritima, a metallo-β-lactamase, is reported. Several crystal structures of wild-type Tm1162 and its mutants were solved. Homologues of Tm1162 are widely distributed in bacteria and archaea, including several human pathogens. The monomer possesses an αβ/βα fold, with the core β-strands having the β-sheet sandwich structure common to the metallo-β-lactamase superfamily. Tm1162 exists as a trimer in the crystal and this trimeric unit is likely to be present in solution. In the trimer, three active sites reside at the interface between subunits, suggesting that the oligomeric assembly is crucial for catalysis. A new type of structurally encoded heterodinuclear site has been identified by confirming the identity of nickel-containing heteronuclear sites in Tm1162 via X-ray absorption spectroscopy and anomalous difference Fourier maps. The second coordination sphere, including His8 and Glu73, maintains the side-chain orientations of histidines and stabilizes the metal-binding site. Nickel coordination was crucial for the oligomerization of Tm1162. The nickel-dependent and manganese-dependent β-lactamase and phosphodiesterase activities of Tm1162 have also been characterized.

Minimum information for reporting next generation sequence genotyping (MIRING): Guidelines for reporting HLA and KIR genotyping via next generation sequencing

The development of next-generation sequencing (NGS) technologies for HLA and KIR genotyping is rapidly advancing knowledge of genetic variation of these highly polymorphic loci. NGS genotyping is poised to replace older methods for clinical use, but standard methods for reporting and exchanging these new, high quality genotype data are needed. The Immunogenomic NGS Consortium, a broad collaboration of histocompatibility and immunogenetics clinicians, researchers, instrument manufacturers and software developers, has developed the Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING) reporting guidelines. MIRING is a checklist that specifies the content of NGS genotyping results as well as a set of messaging guidelines for reporting the results. A MIRING message includes five categories of structured information - message annotation, reference context, full genotype, consensus sequence and novel polymorphism - and references to three categories of accessory information - NGS platform documentation, read processing documentation and primary data. These eight categories of information ensure the long-term portability and broad application of this NGS data for all current histocompatibility and immunogenetics use cases. In addition, MIRING can be extended to allow the reporting of genotype data generated using pre-NGS technologies. Because genotyping results reported using MIRING are easily updated in accordance with reference and nomenclature databases, MIRING represents a bold departure from previous methods of reporting HLA and KIR genotyping results, which have provided static and less-portable data. More information about MIRING can be found online at miring.immunogenomics.org.

Linkage of catalysis and 5' end recognition in ribonuclease RNase J

In diverse bacterial species, the turnover and processing of many RNAs is mediated by the ribonuclease RNase J, a member of the widely occurring metallo-β-lactamase enzyme family. We present crystal structures of Streptomyces coelicolor RNase J with bound RNA in pre- and post-cleavage states, at 2.27 Å and 2.80 Å resolution, respectively. These structures reveal snapshots of the enzyme cleaving substrate directionally and sequentially from the 5′ terminus. In the pre-cleavage state, a water molecule is coordinated to a zinc ion pair in the active site but is imperfectly oriented to launch a nucleophilic attack on the phosphate backbone. A conformational switch is envisaged that enables the in-line positioning of the attacking water and may be facilitated by magnesium ions. Adjacent to the scissile bond, four bases are stacked in a tightly sandwiching pocket, and mutagenesis results indicate that this organization helps to drive processive exo-ribonucleolytic cleavage. Like its numerous homologues, S. coelicolor RNase J can also cleave some RNA internally, and the structural data suggest how the preference for exo- versus endo-cleavage mode is linked with recognition of the chemical status of the substrate's 5′ end.

Theoretical studies of the hydrolysis of antibiotics catalyzed by a metallo-β-lactamase

In this paper, hybrid QM/MM molecular dynamics (MD) simulations have been performed to explore the mechanisms of hydrolysis of two antibiotics, Imipenen (IMI), an antibiotic belonging to the subgroup of carbapenems, and the Cefotaxime (CEF), a third-generation cephalosporin antibiotic, in the active site of a mono-nuclear β-lactamase, CphA from Aeromonas hydrophila. Significant different transition state structures are obtained for the hydrolysis of both antibiotics: while the TS of the CEF is an ionic species with negative charge on nitrogen, the IMI TS presents a tetrahedral-like character with negative charge on oxygen atom of the carbonyl group of the lactam ring. Thus, dramatic conformational changes can take place in the cavity of CphA to accommodate different substrates, which would be the origin of its substrate promiscuity. Since CphA shows only activity against carbapenem antibiotic, this study sheds some light into the origin of the selectivity of the different MbL and, as a consequence, into the discovery of specific and potent MβL inhibitors against a broad spectrum of bacterial pathogens. We have finally probed that a re-parametrization of semiempirical methods should be done to properly describe the behavior the metal cation in active site, Zn(2+), when used in QM/MM calculations.

Emergence of the β-CASP ribonucleases: highly conserved and ubiquitous metallo-enzymes involved in messenger RNA maturation and degradation

The β-CASP ribonucleases, which are found in the three domains of life, have in common a core of 460 residues containing seven conserved sequence motifs involved in the tight binding of two catalytic zinc ions. A hallmark of these enzymes is their ability to catalyze both endo- and exo-ribonucleolytic degradation. Exo-ribonucleolytic degradation proceeds in the 5' to 3' direction and is sensitive to the phosphorylation state of the 5' end of a transcript. Recent phylogenomic analyses have shown that the β-CASP ribonucleases can be partitioned into two major subdivisions that correspond to orthologs of eukaryal CPSF73 and bacterial RNase J. We discuss the known functions of the CPSF73 and RNase J orthologs, their association into complexes, and their structure as it relates to mechanism of action. Eukaryal CPSF73 is part of a large multiprotein complex that is involved in the maturation of the 3' end of RNA Polymerase II transcripts and the polyadenylation of messenger RNA. RNase J1 and J2 are paralogs in Bacillus subtilis that are involved in the degradation of messenger RNA and the maturation of non-coding RNA. RNase J1 and J2 co-purify as a heteromeric complex and there is recent evidence that they interact with other enzymes to form a bacterial RNA degradosome. Finally, we speculate on the evolutionary origin of β-CASP ribonucleases and on their functions in Archaea. Orthologs of CPSF73 with endo- and exo-ribonuclease activity are strictly conserved throughout the archaea suggesting a role for these enzymes in the maturation and/or degradation of messenger RNA. This article is part of a Special Issue entitled: RNA Decay mechanisms.

Automated hierarchical classification of protein domain subfamilies based on functionally-divergent residue signatures

Background

The NCBI Conserved Domain Database (CDD) consists of a collection of multiple sequence alignments of protein domains that are at various stages of being manually curated into evolutionary hierarchies based on conserved and divergent sequence and structural features. These domain models are annotated to provide insights into the relationships between sequence, structure and function via web-based BLAST searches.

Results

Here we automate the generation of conserved domain (CD) hierarchies using a combination of heuristic and Markov chain Monte Carlo (MCMC) sampling procedures and starting from a (typically very large) multiple sequence alignment. This procedure relies on statistical criteria to define each hierarchy based on the conserved and divergent sequence patterns associated with protein functional-specialization. At the same time this facilitates the sequence and structural annotation of residues that are functionally important. These statistical criteria also provide a means to objectively assess the quality of CD hierarchies, a non-trivial task considering that the protein subgroups are often very distantly related—a situation in which standard phylogenetic methods can be unreliable. Our aim here is to automatically generate (typically sub-optimal) hierarchies that, based on statistical criteria and visual comparisons, are comparable to manually curated hierarchies; this serves as the first step toward the ultimate goal of obtaining optimal hierarchical classifications. A plot of runtimes for the most time-intensive (non-parallelizable) part of the algorithm indicates a nearly linear time complexity so that, even for the extremely large Rossmann fold protein class, results were obtained in about a day.

Conclusions

This approach automates the rapid creation of protein domain hierarchies and thus will eliminate one of the most time consuming aspects of conserved domain database curation. At the same time, it also facilitates protein domain annotation by identifying those pattern residues that most distinguish each protein domain subgroup from other related subgroups.

The CphAII protein from Aquifex aeolicus exhibits a metal-dependent phosphodiesterase activity

The CphAII protein from the hyperthermophile Aquifex aeolicus shows the five conserved motifs of the metallo-β-lactamase (MBL) superfamily and presents 28% identity with the Aeromonas hydrophila subclass B2 CphA MBL. The gene encoding CphAII was amplified by PCR from the A. aeolicus genomic DNA and overexpressed in Escherichia coli using a pLex-based expression system. The recombinant CphAII protein was purified by a combination of heating (to denature E. coli proteins) and two steps of immobilized metal affinity chromatography. The purified enzyme preparation did not exhibit a β-lactamase activity but showed a metal-dependent phosphodiesterase activity versus bis-p-nitrophenyl phosphate and thymidine 5'-monophosphate p-nitrophenyl ester, with an optimum at 85°C. The circular dichroism spectrum was in agreement with the percentage of secondary structures characteristic of the MBL αββα fold.

Structure and activity of a novel archaeal β-CASP protein with N-terminal KH domains

MTH1203, a β-CASP metallo-β-lactamase family nuclease from the archaeon Methanothermobacter thermautotrophicus, was identified as a putative nuclease that might contribute to RNA processing. The crystal structure of MTH1203 reveals that, in addition to the metallo-β-lactamase nuclease and the β-CASP domains, it contains two contiguous KH domains that are unique to MTH1203 and its orthologs. RNA-binding experiments indicate that MTH1203 preferentially binds U-rich sequences with a dissociation constant in the micromolar range. In vitro nuclease activity assays demonstrated that MTH1203 is a zinc-dependent nuclease. MTH1203 is also shown to be a dimer and, significantly, this dimerization enhances the nuclease activity. Transcription termination in archaea produces mRNA transcripts with U-rich 3' ends that could be degraded by MTH1203 considering its RNA-binding specificity. We hypothesize that this nuclease degrades mRNAs of proteins targeted for degradation and so regulates archaeal RNA turnover, possibly in concert with the exosome.

Analysis of the functional contributions of Asn233 in metallo-β-lactamase IMP-1

Metallo-β-lactamases, such as IMP-1, are a major global health threat, as they provide for bacterial resistance to a wide range of β-lactam antibiotics, including carbapenems. Understanding the molecular details of the enzymatic process and the sequence requirements for function are essential aids in overcoming β-lactamase-mediated resistance. An asparagine residue is conserved at position 233 in approximately 67% of all metallo-β-lactamases. Despite its conservation, the molecular basis of Asn233 function is poorly understood and remains controversial. It has previously been shown that mutations at this site exhibit context-dependent sequence requirements in that the importance of a given amino acid depends on the antibiotic being tested. To provide a more thorough examination as to the function and sequence requirements at this position, a collection of IMP-1 mutants encoding each of the 19 possible amino acid substitutions was generated. The resistance levels toward four β-lactam antibiotics were measured for Escherichia coli containing each of these mutants. The sequence requirements at position 233 for wild-type levels of resistance toward two cephalosporins were the most relaxed, while there were more stringent sequence requirements for resistance to ampicillin or imipenem. Enzyme kinetic analysis and determinations of steady-state protein levels indicated that the effects of the substitutions on resistance are due to changes in the kinetic parameters of the enzyme. Taken together, the results indicate that substitutions at position 233 significantly alter the kinetic parameters of the enzyme, but most substituted enzymes are able to provide for a high level of resistance to a broad range of β-lactams.

Surveying the manifold divergence of an entire protein class for statistical clues to underlying biochemical mechanisms

Certain residues have no known function yet are co-conserved across distantly related protein families and diverse organisms, suggesting that they perform critical roles associated with as-yet-unidentified molecular properties and mechanisms. This raises the question of how to obtain additional clues regarding these mysterious biochemical phenomena with a view to formulating experimentally testable hypotheses. One approach is to access the implicit biochemical information encoded within the vast amount of genomic sequence data now becoming available. Here, a new Gibbs sampling strategy is formulated and implemented that can partition hundreds of thousands of sequences within a major protein class into multiple, functionally-divergent categories based on those pattern residues that best discriminate between categories. The sampler precisely defines the partition and pattern for each category by explicitly modeling unrelated, non-functional and related-yet-divergent proteins that would otherwise obscure the analysis. To aid biological interpretation, auxiliary routines can characterize pattern residues within available crystal structures and identify those structures most likely to shed light on the roles of pattern residues. This approach can be used to define and annotate automatically subgroup-specific conserved domain profiles based on statistically-rigorous empirical criteria rather than on the subjective and labor-intensive process of manual curation. Incorporating such profiles into domain database search sites (such as the NCBI BLAST site) will provide biologists with previously inaccessible molecular information useful for hypothesis generation and experimental design. Analyses of P-loop GTPases and of AAA+ ATPases illustrate the sampler's ability to obtain such information.

Structural biology of poly(A) site definition

3' processing is an essential step in the maturation of all messenger RNAs (mRNAs) and is a tightly coupled two-step reaction: endonucleolytic cleavage at the poly(A) site is followed by the addition of a poly(A) tail, except for metazoan histone mRNAs, which are cleaved but not polyadenylated. The recognition of a poly(A) site is coordinated by the sequence elements in the mRNA 3' UTR and associated protein factors. In mammalian cells, three well-studied sequence elements, UGUA, AAUAAA, and GU-rich, are recognized by three multisubunit factors: cleavage factor I(m) (CFI(m) ), cleavage and polyadenylation specificity factor (CPSF), and cleavage stimulation factor (CstF), respectively. In the yeast Saccharomyces cerevisiae, UA repeats and A-rich sequence elements are recognized by Hrp1p and cleavage factor IA. Structural studies of protein-RNA complexes have helped decipher the mechanisms underlying sequence recognition and shed light on the role of protein factors in poly(A) site selection and 3' processing machinery assembly. In this review we focus on the interactions between the mRNA cis-elements and the protein factors (CFI(m) , CPSF, CstF, and homologous factors from yeast and other eukaryotes) that define the poly(A) site. WIREs RNA 2011 2 732-747 DOI: 10.1002/wrna.88 For further resources related to this article, please visit the WIREs website.

Identification of an amphipathic helix important for the formation of ectopic septin spirals and axial budding in yeast axial landmark protein Bud3p

Correct positioning of polarity axis in response to internal or external cues is central to cellular morphogenesis and cell fate determination. In the budding yeast Saccharomyces cerevisiae, Bud3p plays a key role in the axial bud-site selection (axial budding) process in which cells assemble the new bud next to the preceding cell division site. Bud3p is thought to act as a component of a spatial landmark. However, it is not clear how Bud3p interacts with other components of the landmark, such as the septins, to control axial budding. Here, we report that overexpression of Bud3p causes the formation of small septin rings (∼1 µm in diameter) and arcs aside from previously reported spiral-like septin structures. Bud3p closely associates with the septins in vivo as Bud3p colocalizes with these aberrant septin structures and forms a complex with two septins, Cdc10p and Cdc11p. The interaction of Bud3p with the septins may involve multiple regions of Bud3p including 1–858, 850–1220, and 1221–1636 a.a. since they all target to the bud neck but exhibit different effects on septin organization when overexpressed. In addition, our study reveals that the axial budding function of Bud3p is mediated by the N-terminal region 1–858. This region shares an amphipathic helix (850–858) crucial for bud neck targeting with the middle portion 850–1103 involved in the formation of ectopic septin spirals and rings. Interestingly, the Dbl-homology domain located in 1–858 is dispensable for axial bud-site selection. Our findings suggest that multiple regions of Bud3p ensure efficient targeting of Bud3p to the bud neck in the assembly of the axial landmark and distinct domains of Bud3p are involved in axial bud-site selection and other cellular processes.

Giardia duodenalis 14-3-3 protein is polyglycylated by a tubulin tyrosine ligase-like member and deglycylated by two metallocarboxypeptidases

The flagellated protozoan Giardia duodenalis is a parasite of the upper part of the small intestine of mammals, including humans, and an interesting biological model. Giardia harbors a single 14-3-3 isoform, a multifunctional protein family, that is modified at the C terminus by polyglycylation, an unusual post-translational modification consisting of the covalent addition of one or multiple glycines on the γ-carboxyl groups of specific glutamic acids. Polyglycylation affects the intracellular localization of g14-3-3, as the shortening of the polyglycine chain is correlated with a partial relocalization of 14-3-3 inside the nuclei during encystation. In this work we demonstrate that the gTTLL3, a member of the tubulin tyrosine ligase-like family, is the enzyme responsible for the 14-3-3 polyglycylation. We also identify two metallopeptidases of the M20 family, here termed gDIP1 (giardial dipeptidase 1) and gDIP2, as enzymes able to shorten the g14-3-3 polyglycine tail both in vivo and in vitro. Finally, we show that the ectopic expression of gDIP2 alters the g14-3-3 localization and strongly hampers the cyst formation. In conclusion, we have identified a polyglycylase and two deglycylases that act in concert to modulate the stage-dependent glycylation status of the multifunctional regulatory g14-3-3 protein in G. duodenalis.

Molecular architecture of the Mn2+-dependent lactonase UlaG reveals an RNase-like metallo-beta-lactamase fold and a novel quaternary structure

The ulaG gene, located in the ula regulon, is crucial for the catabolism of l-ascorbate under anaerobic conditions and it has been proposed to encode for the putative l-ascorbate-6-P lactonase. The ulaG gene is widespread among eubacteria, including human commensal and pathogenic genera such as Escherichia, Shigella, Klebsiella and Salmonella. Here, we report the three-dimensional structures of the apoenzyme and Mn(2+) holoenzyme of UlaG from E. coli to 2.6 A resolution, determined using single-wavelength anomalous diffraction phasing and molecular replacement, respectively. The structures reveal a highly specialized metallo-beta-lactamase-like fold derived from an ancient structural template that was involved in RNA maturation and DNA repair. This fold has a novel quaternary architecture consisting of a hexameric ring formed by a trimer of UlaG dimers. A mononuclear Mn(2)(+)-binding site resides at the core of the active site, which displays micromolar affinity for Mn(2+) and a distorted trigonal bipyramidal coordination. The active site Mn(2+) ion can be replaced by Co(2+) or Zn(2+), but not by Fe(3+). We further show that the Mn(2+) or Co(2)(+)-loaded enzyme exhibits lactonase activity towards l-ascorbate 6-P, thereby providing the first direct evidence of its catalytic role in the L-ascorbate catabolic pathway. Guided by the structural homology, we show that UlaG is able to cleave phosphodiester linkages in cyclic nucleotides, suggesting that the conservation of the fold and of the key catalytic residues allows for the evolutionary acquisition of substrate specificity for novel but related substrates.

Rapid detection, classification and accurate alignment of up to a million or more related protein sequences

MOTIVATION

The patterns of sequence similarity and divergence present within functionally diverse, evolutionarily related proteins contain implicit information about corresponding biochemical similarities and differences. A first step toward accessing such information is to statistically analyze these patterns, which, in turn, requires that one first identify and accurately align a very large set of protein sequences. Ideally, the set should include many distantly related, functionally divergent subgroups. Because it is extremely difficult, if not impossible for fully automated methods to align such sequences correctly, researchers often resort to manual curation based on detailed structural and biochemical information. However, multiply-aligning vast numbers of sequences in this way is clearly impractical.

RESULTS

This problem is addressed using Multiply-Aligned Profiles for Global Alignment of Protein Sequences (MAPGAPS). The MAPGAPS program uses a set of multiply-aligned profiles both as a query to detect and classify related sequences and as a template to multiply-align the sequences. It relies on Karlin-Altschul statistics for sensitivity and on PSI-BLAST (and other) heuristics for speed. Using as input a carefully curated multiple-profile alignment for P-loop GTPases, MAPGAPS correctly aligned weakly conserved sequence motifs within 33 distantly related GTPases of known structure. By comparison, the sequence- and structurally based alignment methods hmmalign and PROMALS3D misaligned at least 11 and 23 of these regions, respectively. When applied to a dataset of 65 million protein sequences, MAPGAPS identified, classified and aligned (with comparable accuracy) nearly half a million putative P-loop GTPase sequences.

AVAILABILITY

A C++ implementation of MAPGAPS is available at http://mapgaps.igs.umaryland.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The charge-dipole pocket: a defining feature of signaling pathway GTPase on/off switches

Ras-like GTPases function as on/off switches in intracellular signaling pathways. Their on or off state is communicated through conformational changes in the so-called switch I and II regions. It is commonly believed that the distinguishing molecular features of these GTPases are well known. Here, however, I identify-through a Bayesian iterative analysis of GTPase evolutionary divergence-a previously undescribed switch II structural component that (along with previously described, functionally critical residues) most distinguish these signaling pathway on/off switches from other GTPases. In certain Ras-like GTPases this newly-identified component forms an aromatic pocket around the negative-dipole moment at the end of a switch II helix with a positively charged residue inserted into the pocket. This helix is oriented in a specific direction away from the GTPase core, but is reoriented dramatically upon disruption of the charge-dipole pocket. The charge-dipole pocket occurs in both the on and off states and both the charge-dipole pocket and an alternative configuration occur within the unit cell of a single crystal structure of Rab5a GTPase in the off state. Thus, the charge-dipole pocket configuration is closely associated, not with the on or off state, but rather with formation of the outward-oriented helix and, as a result, with restructuring of the switch II N-terminal region, which has a critical role both in sensing the on/off state and in mediating GTP hydrolysis and nucleotide exchange.

The glycine brace: a component of Rab, Rho, and Ran GTPases associated with hinge regions of guanine- and phosphate-binding loops

Background

Ras-like GTPases function as on-off switches in intracellular signalling pathways and include the Rab, Rho/Rac, Ran, Ras, Arf, Sar and Gα families. How these families have evolutionarily diverged from each other at the sequence level provides clues to underlying mechanisms associated with their functional specialization.

Results

Bayesian analysis of divergent patterns within a multiple alignment of Ras-like GTPase sequences identifies a structural component, termed here the glycine brace, as the feature that most distinguishes Rab, Rho/Rac, Ran and (to some degree) Ras family GTPases from other Ras-like GTPases. The glycine brace consists of four residues: An aromatic residue that forms a stabilizing CH-π interaction with a conserved glycine at the start of the guanine-binding loop; a second aromatic residue, which is nearly always a tryptophan, that likewise forms stabilizing CH-π and NH-π interactions with a glycine at the start of the phosphate-binding P-loop; and two other residues (typically an aspartate and a serine or threonine) that, together with a conserved buried water molecule, form a network of interactions connecting the two aromatic residues.

Conclusion

It is proposed that the two glycine residues function as hinges and that the glycine brace influences guanine nucleotide binding and release by interacting with these hinges.

Protein stabilization by the rational design of surface charge-charge interactions

The design of proteins with increased stability has many important applications in biotechnology. In recent years, strategies involving directed evolution, sequence-based design, or computational design have proven successful for generating stabilized proteins. A brief overview of the various methods that have been used to increase protein stability is presented, followed by a detailed example of how the rational design of surface charge-charge interactions has provided a robust method for protein stabilization.

ARCS-Motif: discovering correlated motifs from unaligned biological sequences

MOTIVATION

The goal of motif discovery is to detect novel, unknown, and important signals from biology sequences. In most models, the importance of a motif is equal to the sum of the similarity of every single position. In 2006, Song et al. introduced Aggregated Related Column Score (ARCS) measure which includes correlation information to the evaluation of motif importance. The paper showed that the ARCS measure is superior to other measures. Due to the complicated nature of the ARCS motif model, we cannot directly apply existing sequential motif discovery methods to find motifs with high ARCS values.

RESULTS

This article presents a novel mining algorithm, ARCS-Motif, to discover related sequential motifs in biological sequences. ARCS-Motif is applied to 400 PROSITE datasets and compared with five alternative methods (CONSENSUS, Gibbs sampler, MEME, SPLASH and DIALIGN-TX). ARCS-Motif outperforms all the methods in accuracy, and most of the methods in efficiency. Although SPLASH has better efficiency than ARCS-Motif, ARCS-Motif has much better accuracy than SPLASH. On average, ARCS-Motif is able to produce the motifs which are at least 10% better than the best of the alternative methods. Among the 400 PROSITE datasets, ARCS-Motif produces the best motifs for more than 200 families. Other than SPLASH, the execution time of ARCS-Motif is less than a third of that of the fastest alternative method and its execution time grows at the slowest rate with respect to the number of sequences and the average sequence among all methods.

Structural and functional diversity of bacterial membrane fusion proteins

Membrane Fusion Proteins (MFPs) are functional subunits of multi-component transporters that perform diverse physiological functions in both Gram-positive and Gram-negative bacteria. MFPs associate with transporters belonging to Resistance-Nodulation-cell Division (RND), ATP-Binding Cassette (ABC) and Major Facilitator (MF) superfamilies of proteins. Recent studies suggested that MFPs interact with substrates and play an active role in transport reactions. In addition, the MFP-dependent transporters from Gram-negative bacteria recruit the outer membrane channels to expel various substrates across the outer membrane into external medium. This review is focused on the diversity, structure and molecular mechanism of MFPs that function in multidrug efflux. Using phylogenetic approaches we analyzed diversity and representation of multidrug MFPs in sequenced bacterial genomes. In addition to previously characterized MFPs from Gram-negative bacteria, we identified MFPs that associate with RND-, MF- and ABC-type transporters in Gram-positive bacteria. Sequence analyses showed that MFPs vary significantly in size (200-650 amino acid residues) with some of them lacking the signature alpha-helical domain of multidrug MFPs. Furthermore, many transport operons contain two- or three genes encoding distinct MFPs. We further discuss the diversity of MFPs in the context of current views on the mechanism and structure of MFP-dependent transporters.

Discovering sequence motifs with arbitrary insertions and deletions

BIOLOGY IS ENCODED IN MOLECULAR SEQUENCES: deciphering this encoding remains a grand scientific challenge. Functional regions of DNA, RNA, and protein sequences often exhibit characteristic but subtle motifs; thus, computational discovery of motifs in sequences is a fundamental and much-studied problem. However, most current algorithms do not allow for insertions or deletions (indels) within motifs, and the few that do have other limitations. We present a method, GLAM2 (Gapped Local Alignment of Motifs), for discovering motifs allowing indels in a fully general manner, and a companion method GLAM2SCAN for searching sequence databases using such motifs. glam2 is a generalization of the gapless Gibbs sampling algorithm. It re-discovers variable-width protein motifs from the PROSITE database significantly more accurately than the alternative methods PRATT and SAM-T2K. Furthermore, it usefully refines protein motifs from the ELM database: in some cases, the refined motifs make orders of magnitude fewer overpredictions than the original ELM regular expressions. GLAM2 performs respectably on the BAliBASE multiple alignment benchmark, and may be superior to leading multiple alignment methods for "motif-like" alignments with N- and C-terminal extensions. Finally, we demonstrate the use of GLAM2 to discover protein kinase substrate motifs and a gapped DNA motif for the LIM-only transcriptional regulatory complex: using GLAM2SCAN, we identify promising targets for the latter. GLAM2 is especially promising for short protein motifs, and it should improve our ability to identify the protein cleavage sites, interaction sites, post-translational modification attachment sites, etc., that underlie much of biology. It may be equally useful for arbitrarily gapped motifs in DNA and RNA, although fewer examples of such motifs are known at present. GLAM2 is public domain software, available for download at http://bioinformatics.org.au/glam2.

The CHAIN program: forging evolutionary links to underlying mechanisms

Proteins evolve new functions by modifying and extending the molecular machinery of an ancestral protein. Such changes show up as divergent sequence patterns, which are conserved in descendent proteins that maintain the divergent function. After multiply-aligning a set of input sequences, the CHAIN program partitions the sequences into two functionally divergent groups and then outputs an alignment that is annotated to reveal the selective pressures imposed on divergent residue positions. If atomic coordinates are also provided, hydrogen bonds and other atomic interactions associated with various categories of divergent residues are graphically displayed. Such analyses establish links between protein evolutionary divergence and functionally crucial atomic features and, as a result, can suggest plausible molecular mechanisms for experimental testing. This is illustrated here by its application to bacterial clamp-loader ATPases.

The global trace graph, a novel paradigm for searching protein sequence databases

MOTIVATION

Propagating functional annotations to sequence-similar, presumably homologous proteins lies at the heart of the bioinformatics industry. Correct propagation is crucially dependent on the accurate identification of subtle sequence motifs that are conserved in evolution. The evolutionary signal can be difficult to detect because functional sites may consist of non-contiguous residues while segments in-between may be mutated without affecting fold or function.

RESULTS

Here, we report a novel graph clustering algorithm in which all known protein sequences simultaneously self-organize into hypothetical multiple sequence alignments. This eliminates noise so that non-contiguous sequence motifs can be tracked down between extremely distant homologues. The novel data structure enables fast sequence database searching methods which are superior to profile-profile comparison at recognizing distant homologues. This study will boost the leverage of structural and functional genomics and opens up new avenues for data mining a complete set of functional signature motifs.

AVAILABILITY

http://www.bioinfo.biocenter.helsinki.fi/gtg.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Metallo-beta-lactamases (classification, activity, genetic organization, structure, zinc coordination) and their superfamily

One strategy employed by bacterial strains to resist beta-lactam antibiotics is the expression of metallo-beta-lactamases requiring Zn(2+) for activity. In the last few years, many new zinc beta-lactamases have been described and several pathogens are now known to synthesize members of this class. Metallo-beta-lactamases are especially worrisome due to: (1) their broad activity profiles that encompass most beta-lactam antibiotics, including the carbapenems; (2) potential for horizontal transference; and (3) the absence of clinically useful inhibitors. On the basis of the known sequences, three different lineages, identified as subclasses B1, B2, and B3 have been characterized. The three-dimensional structure of at least one metallo-beta-lactamase of each subclass has been solved. These very similar 3D structures are characterized by the presence of an alphabetabetaalpha-fold. In addition to metallo-beta-lactamases which cleave the amide bond of the beta-lactam ring, the metallo-beta-lactamase superfamily includes enzymes which hydrolyze thiol-ester, phosphodiester and sulfuric ester bonds as well as oxydoreductases. Most of the 6000 members of this superfamily share five conserved motifs, the most characteristic being the His116-X-His118-X-Asp120-His121 signature. They all exhibit an alphabetabetaalpha-fold, similar to that found in the structure of zinc beta-lactamases. Many members of this superfamily are involved in mRNA maturation and DNA reparation. This fact suggests the hypothesis that metallo-beta-lactamases may be the result of divergent evolution starting from an ancestral protein which did not have a beta-lactamase activity.

COBALT: constraint-based alignment tool for multiple protein sequences

MOTIVATION

A tool that simultaneously aligns multiple protein sequences, automatically utilizes information about protein domains, and has a good compromise between speed and accuracy will have practical advantages over current tools.

RESULTS

We describe COBALT, a constraint based alignment tool that implements a general framework for multiple alignment of protein sequences. COBALT finds a collection of pairwise constraints derived from database searches, sequence similarity and user input, combines these pairwise constraints, and then incorporates them into a progressive multiple alignment. We show that using constraints derived from the conserved domain database (CDD) and PROSITE protein-motif database improves COBALT's alignment quality. We also show that COBALT has reasonable runtime performance and alignment accuracy comparable to or exceeding that of other tools for a broad range of problems.

AVAILABILITY

COBALT is included in the NCBI C++ toolkit. A Linux executable for COBALT, and CDD and PROSITE data used is available at: ftp://ftp.ncbi.nlm.nih.gov/pub/agarwala/cobalt

Identification and characterization of archaeal and fungal tRNA methyltransferases

All organisms modify their tRNAs by use of evolutionarily conserved enzymes. Members of the Archaea contain an extensive set of modified nucleotides that were early evidence of the fundamental evolutionary divergence of the Archaea from Bacteria and Eucarya. However, the enzymes responsible for these posttranscriptional modifications were largely unknown before the advent of genome sequencing. This chapter explains methods to identify tRNA methyltransferases in genome sequences, emphasizing the identification and characterization of six enzymes from the hyperthermophilic archaeon Methanocaldococcus jannaschii. We describe methods to express these proteins, purify or synthesize tRNA substrates, measure methyltransferase activity, and map tRNA modifications. Comparison of the archaeal methyltransferases with their yeast homologs suggests that the common ancestor of the archaeal and eucaryal organismal lineages already had extensive tRNA modifications.

Convergent evolution of receptors for protein import into mitochondria

BACKGROUND

Mitochondria evolved from intracellular bacterial symbionts. Establishing mitochondria as organelles required a molecular machine to import proteins across the mitochondrial outer membrane. This machinery, the TOM complex, is composed of at least seven component parts, and its creation and evolution represented a sizeable challenge. Although there is good evidence that a core TOM complex, composed of three subunits, was established in the protomitochondria, we suggest that the receptor component of the TOM complex arose later in the evolution of this machine.

RESULTS

We have solved by nuclear magnetic resonance the structure of the presequence binding receptor from the TOM complex of the plant Arabidopsis thaliana. The protein fold suggests that this protein, AtTom20, belongs to the tetratricopeptide repeat (TPR) superfamily, but it is unusual in that it contains insertions lengthening the helices of each TPR motif. Peptide titrations map the presequence binding site to a groove of the concave surface of the receptor. In vitro functional assays and peptide titrations suggest that the plant Tom20 is functionally equivalent to fungal and animal Tom20s.

CONCLUSIONS

Comparison of the sequence and structure of Tom20 from plants and animals suggests that these two presequence binding receptors evolved from two distinct ancestral genes following the split of the animal and plant lineages. The need to bind equivalent mitochondrial targeting sequences and to make similar interactions within an equivalent protein translocation machine has driven the convergent evolution of two distinct proteins to a common structure and function.

MUSTANG: a multiple structural alignment algorithm

Multiple structural alignment is a fundamental problem in structural genomics. In this article, we define a reliable and robust algorithm, MUSTANG (MUltiple STructural AligNment AlGorithm), for the alignment of multiple protein structures. Given a set of protein structures, the program constructs a multiple alignment using the spatial information of the C(alpha) atoms in the set. Broadly based on the progressive pairwise heuristic, this algorithm gains accuracy through novel and effective refinement phases. MUSTANG reports the multiple sequence alignment and the corresponding superposition of structures. Alignments generated by MUSTANG are compared with several handcurated alignments in the literature as well as with the benchmark alignments of 1033 alignment families from the HOMSTRAD database. The performance of MUSTANG was compared with DALI at a pairwise level, and with other multiple structural alignment tools such as POSA, CE-MC, MALECON, and MultiProt. MUSTANG performs comparably to popular pairwise and multiple structural alignment tools for closely related proteins, and performs more reliably than other multiple structural alignment methods on hard data sets containing distantly related proteins or proteins that show conformational changes.

Structural and functional evolution of transthyretin and transthyretin-like proteins

Transthyretin (TTR) is a tetrameric protein involved in the distribution of thyroid hormones in vertebrates. The amino acid sequence of TTR is highly conserved across vertebrates. Hypothetical TTR-like proteins (TLPs) were inferred from the identification of genes in nonvertebrate species. Here, we identified five motifs defining TLPs and three motifs defining both TTRs and TLPs. These motifs were mapped onto structurally conserved and functionally important regions of TTRs. These motifs were used to build hidden Markov models for accurate identification of TLPs in other organisms. TLPs were divided into three main groups based on their N-terminal regions. Most TLPs are cytosolic, but in plants and slime mold, we predict they are peroxisomal. We verified that the TLPs from enterobacteria were periplasmic. We demonstrated that TLP genes are expressed in a bacterium (E. coli), an invertebrate animal (C. elegans), and a plant (A. thaliana). These TLPs have similar subunit molecular weights to TTRs, are tetramers, and are predicted to have similar three-dimensional (3D) structures to TTRs, but do not bind thyroid hormones or similar ligands. We suggest that like TTRs, the N-terminal and C-terminal regions of TLPs are integral in defining the function of TLPs in nonvertebrate species and that the TLP gene duplicated in primitive vertebrates to produce the TTR gene. TLP/TTR has retained its overall structure, but changed function and localization during evolution in bacteria, invertebrates, plants, and vertebrates.

Crystal structure of TDP-fucosamine acetyltransferase (WecD) from Escherichia coli, an enzyme required for enterobacterial common antigen synthesis

Enterobacterial common antigen (ECA) is a polysaccharide found on the outer membrane of virtually all gram-negative enteric bacteria and consists of three sugars, N-acetyl-d-glucosamine, N-acetyl-d-mannosaminuronic acid, and 4-acetamido-4,6-dideoxy-d-galactose, organized into trisaccharide repeating units having the sequence -->3)-alpha-d-Fuc4NAc-(1-->4)-beta-d-ManNAcA-(1-->4)-alpha-d-GlcNAc-(1-->. While the precise function of ECA is unknown, it has been linked to the resistance of Shiga-toxin-producing Escherichia coli (STEC) O157:H7 to organic acids and the resistance of Salmonella enterica to bile salts. The final step in the synthesis of 4-acetamido-4,6-dideoxy-d-galactose, the acetyl-coenzyme A (CoA)-dependent acetylation of the 4-amino group, is carried out by TDP-fucosamine acetyltransferase (WecD). We have determined the crystal structure of WecD in apo form at a 1.95-Angstrom resolution and bound to acetyl-CoA at a 1.66-Angstrom resolution. WecD is a dimeric enzyme, with each monomer adopting the GNAT N-acetyltransferase fold, common to a number of enzymes involved in acetylation of histones, aminoglycoside antibiotics, serotonin, and sugars. The crystal structure of WecD, however, represents the first structure of a GNAT family member that acts on nucleotide sugars. Based on this cocrystal structure, we have used flexible docking to generate a WecD-bound model of the acetyl-CoA-TDP-fucosamine tetrahedral intermediate, representing the structure during acetyl transfer. Our structural data show that WecD does not possess a residue that directly functions as a catalytic base, although Tyr208 is well positioned to function as a general acid by protonating the thiolate anion of coenzyme A.

Motif discoveries in unaligned molecular sequences using self-organizing neural networks

In this paper, we study the problem of motif discoveries in unaligned DNA and protein sequences. The problem of motif identification in DNA and protein sequences has been studied for many years in the literature. Major hurdles at this point include computational complexity and reliability of the search algorithms. We propose a self-organizing neural network structure for solving the problem of motif identification in DNA and protein sequences. Our network contains several layers, with each layer performing classifications at different levels. The top layer divides the input space into a small number of regions and the bottom layer classifies all input patterns into motifs and nonmotif patterns. Depending on the number of input patterns to be classified, several layers between the top layer and the bottom layer are needed to perform intermediate classifications. We maintain a low computational complexity through the use of the layered structure so that each pattern's classification is performed with respect to a small subspace of the whole input space. Our self-organizing neural network will grow as needed (e.g., when more motif patterns are classified). It will give the same amount of attention to each input pattern and will not omit any potential motif patterns. Finally, simulation results show that our algorithm outperforms existing algorithms in certain aspects. In particular, simulation results show that our algorithm can identify motifs with more mutations than existing algorithms. Our algorithm works well for long DNA sequences as well.

Protein classification based on text document classification techniques

The need for accurate, automated protein classification methods continues to increase as advances in biotechnology uncover new proteins. G-protein coupled receptors (GPCRs) are a particularly difficult superfamily of proteins to classify due to extreme diversity among its members. Previous comparisons of BLAST, k-nearest neighbor (k-NN), hidden markov model (HMM) and support vector machine (SVM) using alignment-based features have suggested that classifiers at the complexity of SVM are needed to attain high accuracy. Here, analogous to document classification, we applied Decision Tree and Naive Bayes classifiers with chi-square feature selection on counts of n-grams (i.e. short peptide sequences of length n) to this classification task. Using the GPCR dataset and evaluation protocol from the previous study, the Naive Bayes classifier attained an accuracy of 93.0 and 92.4% in level I and level II subfamily classification respectively, while SVM has a reported accuracy of 88.4 and 86.3%. This is a 39.7 and 44.5% reduction in residual error for level I and level II subfamily classification, respectively. The Decision Tree, while inferior to SVM, outperforms HMM in both level I and level II subfamily classification. For those GPCR families whose profiles are stored in the Protein FAMilies database of alignments and HMMs (PFAM), our method performs comparably to a search against those profiles. Finally, our method can be generalized to other protein families by applying it to the superfamily of nuclear receptors with 94.5, 97.8 and 93.6% accuracy in family, level I and level II subfamily classification respectively.

A phylogenomic profile of globins

Background

Globins occur in all three kingdoms of life: they can be classified into single-domain globins and chimeric globins. The latter comprise the flavohemoglobins with a C-terminal FAD-binding domain and the gene-regulating globin coupled sensors, with variable C-terminal domains. The single-domain globins encompass sequences related to chimeric globins and «truncated» hemoglobins with a 2-over-2 instead of the canonical 3-over-3 α-helical fold.

Results

A census of globins in 26 archaeal, 245 bacterial and 49 eukaryote genomes was carried out. Only ~25% of archaea have globins, including globin coupled sensors, related single domain globins and 2-over-2 globins. From one to seven globins per genome were found in ~65% of the bacterial genomes: the presence and number of globins are positively correlated with genome size. Globins appear to be mostly absent in Bacteroidetes/Chlorobi, Chlamydia, Lactobacillales, Mollicutes, Rickettsiales, Pastorellales and Spirochaetes. Single domain globins occur in metazoans and flavohemoglobins are found in fungi, diplomonads and mycetozoans. Although red algae have single domain globins, including 2-over-2 globins, the green algae and ciliates have only 2-over-2 globins. Plants have symbiotic and nonsymbiotic single domain hemoglobins and 2-over-2 hemoglobins. Over 90% of eukaryotes have globins: the nematode Caenorhabditis has the most putative globins, ~33. No globins occur in the parasitic, unicellular eukaryotes such as Encephalitozoon, Entamoeba, Plasmodium and Trypanosoma.

Conclusion

Although Bacteria have all three types of globins, Archaeado not have flavohemoglobins and Eukaryotes lack globin coupled sensors. Since the hemoglobins in organisms other than animals are enzymes or sensors, it is likely that the evolution of an oxygen transport function accompanied the emergence of multicellular animals.

Diverse membrane‐associated proteins contain a novel SMP domain

We have analyzed the sequence of a mitochondrial integral membrane protein, Mdm12, and found that it forms the prototype for a novel domain, designated the SMP domain, that is common to an extended family of membrane-associated proteins. Comprehensive sequence searches using protein alignment models of SMP proteins were cross-validated by statistical resampling; providing strong support for these relationships. No consensus of 3-dimensional structure was reached upon threading sequences through known folds. SMP proteins are widespread amongst eukaryotic species with a particular enrichment in plants and features suggestive of species-specific functional variations. Members of 2 SMP families, the mitochore and tricalbin proteins, are essential components of protein complexes involved in mitochondrial inheritance and receptor endocytosis while a third SMP protein family, HT008, is associated with the Rvs161-Rvs167 complex, a known regulator of sphingolipid metabolism. In addition, HT008 and PDZK8 SMP proteins possess additional protein-protein interaction domains in domain architectures that are typical of protein scaffolds and adaptors. We therefore predict that the SMP domain is an important link between these distinct membrane-associated proteins and a key regulatory hub for unidentified global regulators.-Lee, I., Hong, W. Diverse membrane-associated proteins contain a novel SMP domain.

MKLP1 requires specific domains for its dendritic targeting

The mitotic kinesin-like protein 1 (MKLP1) is specifically localized to the dendrite of the developed neuron, but its targeting mechanism is still unclear. In this study, the role of distinct domains of MKLP1 in dendritic targeting was investigated by producing a series of enhanced green fluorescent protein (eGFP)-tagged MKLP1 and its variant mutations, and studying the distribution of these molecules in cultured primary hippocampal neurons using fluorescence microscopy. We have found that: (a) full-length MKLP1(1-856)-eGFP was distributed to the dendrite of hippocampal neurons, not the axon; (b) deletion of the ;motor' domain prevented the dendritic distribution; (c) deletion of the ;tail' domain caused axonal appearance; (d) mutants devoid of the ;stalk' domain were still specifically distributed to the dendrite. The results indicate that the motor and tail regions of MKLP1 are important and significant for its localization to the dendrites. We also discuss the difference between the targeting of membrane-anchoring proteins and the kinesin-like protein.

Structure and function of lineage-specific sequence insertions in the bacterial RNA polymerase beta' subunit

The large beta and beta' subunits of the bacterial core RNA polymerase (RNAP) are highly conserved throughout evolution. Nevertheless, large sequence insertions in beta and beta' characterize specific evolutionary lineages of bacteria. The Thermus aquaticus RNAP beta' subunit contains a 283 residue insert between conserved regions A and B that is found in only four bacterial species. The Escherichia coli RNAP beta' subunit contains a 188 residue insert in the middle of conserved region G that is found in a wide range of bacterial species. Here, we present structural studies of these two beta' insertions. We show that the inserts comprise repeats of a previously characterized fold, the sandwich-barrel hybrid motif (as predicted from previous sequence analysis) and that the inserts serve significant roles in facilitating protein/protein and/or protein/nucleic acid interactions.

A robust method to detect structural and functional remote homologues

With currently available sequence data, it is feasible to conduct extensive comparisons among large sets of protein sequences. It is still a much more challenging task to partition the protein space into structurally and functionally related families solely based on sequence comparisons. The ProtoNet system automatically generates a treelike classification of the whole protein space. It stands to reason that this classification reflects evolutionary relationships, both close and remote. In this article, we examine this hypothesis. We present a semiautomatic procedure that singles out certain inner nodes in the ProtoNet tree that should ideally correspond to structurally and functionally defined protein families. We compare the performance of this method against several expert systems. Some of the competing methods incorporate additional extraneous information on protein structure or on enzymatic activities. The ProtoNet-based method performs at least as well as any of the methods with which it was compared. This article illustrates the ProtoNet-based method on several evolutionarily diverse families. Using this new method, an evolutionary divergence scheme can be proposed for a large number of structural and functional related superfamilies.