Crystal structure of conserved hypothetical protein Aq1575 from Aquifex aeolicus

Gene Regulation and Transcriptomics

Borrelia (Borreliella) burgdorferi, along with closely related species, is the etiologic agent of Lyme disease. The spirochete subsists in an enzootic cycle that encompasses acquisition from a vertebrate host to a tick vector and transmission from a tick vector to a vertebrate host. To adapt to its environment and persist in each phase of its enzootic cycle, B. burgdorferi wields three systems to regulate the expression of genes: the RpoN-RpoS alternative sigma factor cascade, the Hk1/Rrp1 two-component system and its product c-di-GMP, and the stringent response mediated by Rel_Bbu and DksA. These regulatory systems respond to enzootic phase-specific signals and are controlled or fine- tuned by transcription factors, including BosR and BadR, as well as small RNAs, including DsrABb and Bb6S RNA. In addition, several other DNA-binding and RNA-binding proteins have been identified, although their functions have not all been defined. Global changes in gene expression revealed by high-throughput transcriptomic studies have elucidated various regulons, albeit technical obstacles have mostly limited this experimental approach to cultivated spirochetes. Regardless, we know that the spirochete, which carries a relatively small genome, regulates the expression of a considerable number of genes required for the transitions between the tick vector and the vertebrate host as well as the adaptation to each.

Elucidation of a Novel Role of YebC in Surface Polysaccharides Regulation of Escherichia coli bipA-Deletion

The BipA (BPI-inducible protein A) protein is ubiquitously conserved in various bacterial species and belongs to the translational GTPase family. Interestingly, the function of Escherichia coli BipA is not essential for cell growth under normal growth conditions. However, cultivation of bipA-deleted cells at 20°C leads to cold-sensitive growth defect and several phenotypic changes in ribosome assembly, capsule production, and motility, suggesting its global regulatory roles. Previously, our genomic library screening revealed that the overexpressed ribosomal protein (r-protein) L20 partially suppressed cold-sensitive growth defect by resolving the ribosomal abnormality in bipA-deleted cells at low temperature. Here, we explored another genomic library clone containing yebC, which encodes a predicted transcriptional factor that is not directly associated with ribosome biogenesis. Interestingly, overexpression of yebC in bipA-deleted cells diminished capsule synthesis and partially restored lipopolysaccharide (LPS) core maturation at a low temperature without resolving defects in ribosome assembly or motility, indicating that YebC may be specifically involved in the regulation of exopolysaccharide and LPS core synthesis. In this study, we collectively investigated the impacts of bipA-deletion on E. coli capsule, LPS, biofilm formation, and motility and revealed novel roles of YebC in extracellular polysaccharide production and LPS core synthesis at low temperature using this mutant strain. Furthermore, our findings suggest that ribosomal defects as well as increased capsule synthesis, and changes in LPS composition may contribute independently to the cold-sensitivity of bipA-deleted cells, implying multiple regulatory roles of BipA.

YebC regulates variable surface antigen VlsE expression and is required for host immune evasion in Borrelia burgdorferi

Borrelia burgdorferi, the Lyme disease pathogen causes persistent infection by evading the host immune response. Differential expression of the surface-exposed lipoprotein VlsE that undergoes antigenic variation is a key immune evasion strategy employed by B. burgdorferi. Most studies focused on the mechanism of VlsE antigen variation, but little is known about VlsE regulation and factor(s) that regulates differential vlsE expression. In this study, we investigated BB0025, a putative YebC family transcriptional regulator (and hence designated BB0025 as YebC of B. burgdorferi herein). We constructed yebC mutant and complemented strain in an infectious strain of B. burgdorferi. The yebC mutant could infect immunocompromised SCID mice but not immunocompetent mice, suggesting that YebC plays an important role in evading host adaptive immunity. RNA-seq analyses identified vlsE as one of the genes whose expression was most affected by YebC. Quantitative RT-PCR and Western blot analyses confirmed that vlsE expression was dependent on YebC. In vitro, YebC and VlsE were co-regulated in response to growth temperature. In mice, both yebC and vlsE were inversely expressed with ospC in response to the host adaptive immune response. Furthermore, EMSA proved that YebC directly binds to the vlsE promoter, suggesting a direct transcriptional control. These data demonstrate that YebC is a new regulator that modulates expression of vlsE and other genes important for spirochetal infection and immune evasion in the mammalian host.

Loss of the RNA-binding protein TACO1 causes late-onset mitochondrial dysfunction in mice

The recognition and translation of mammalian mitochondrial mRNAs are poorly understood. To gain further insights into these processes in vivo, we characterized mice with a missense mutation that causes loss of the translational activator of cytochrome oxidase subunit I (TACO1). We report that TACO1 is not required for embryonic survival, although the mutant mice have substantially reduced COXI protein, causing an isolated complex IV deficiency. We show that TACO1 specifically binds the mt-Co1 mRNA and is required for translation of COXI through its association with the mitochondrial ribosome. We determined the atomic structure of TACO1, revealing three domains in the shape of a hook with a tunnel between domains 1 and 3. Mutations in the positively charged domain 1 reduce RNA binding by TACO1. The Taco1 mutant mice develop a late-onset visual impairment, motor dysfunction and cardiac hypertrophy and thus provide a useful model for future treatment trials for mitochondrial disease.

Mutations in the translational activator of cytochrome c oxidase subunit I (TACO1) causes cytochrome c oxidase deficiency and Leigh Syndrome in patients. Here, the authors characterize mice with a mutation that causes lack of TACO1 expression, identifying a mouse model that could be useful for preclinical trials.

In silico identification of a multi-functional regulatory protein involved in Holliday junction resolution in bacteria

BACKGROUND

Homologous recombination is a fundamental cellular process that is most widely used by cells to rearrange genes and accurately repair DNA double-strand breaks. It may result in the formation of a critical intermediate named Holliday junction, which is a four-way DNA junction and needs to be resolved to allow chromosome segregation. Different Holliday junction resolution systems and enzymes have been characterized from all three domains of life. In bacteria, the RuvABC complex is the most important resolution system.

RESULTS

In this study, we conducted comparative genomics studies to identify a novel DNA-binding protein, YebC, which may serve as a key transcriptional regulator that mainly regulates the gene expression of RuvABC resolvasome in bacteria. On the other hand, the presence of YebC orthologs in some organisms lacking RuvC implied that it might participate in other biological processes. Further phylogenetic analysis of YebC protein sequences revealed two functionally different subtypes: YebC_I and YebC_II. Distribution of YebC_I is much wider than YebC_II. Only YebC_I proteins may play an important role in regulating RuvABC gene expression in bacteria. Investigation of YebC-like proteins in eukaryotes suggested that they may have originated from YebC_II proteins and evolved a new function as a specific translational activator in mitochondria. Finally, additional phylum-specific genes associated with Holliday junction resolution were predicted.

CONCLUSIONS

Overall, our data provide new insights into the basic mechanism of Holliday junction resolution and homologous recombination in bacteria.

From complete genome sequence to 'complete' understanding?

The rapidly accumulating genome sequence data allow researchers to address fundamental biological questions that were not even asked just a few years ago. A major problem in genomics is the widening gap between the rapid progress in genome sequencing and the comparatively slow progress in the functional characterization of sequenced genomes. Here we discuss two key questions of genome biology: whether we need more genomes, and how deep is our understanding of biology based on genomic analysis. We argue that overly specific annotations of gene functions are often less useful than the more generic, but also more robust, functional assignments based on protein family classification. We also discuss problems in understanding the functions of the remaining 'conserved hypothetical' genes.

Rapid model building of alpha-helices in electron-density maps

A method for the identification of alpha-helices in electron-density maps at low resolution followed by interpretation at moderate to high resolution is presented. Rapid identification is achieved at low resolution, where alpha-helices appear as tubes of density. The positioning and direction of the alpha-helices is obtained at moderate to high resolution, where the positions of side chains can be seen. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. An average of 63% of the alpha-helical residues in these proteins were built and an average of 76% of the residues built matched helical residues in the refined models of the proteins. The overall average r.m.s.d. between main-chain atoms in the modeled alpha-helices and the nearest atom with the same name in the refined models of the proteins was 1.3 A.

Rapid model building of beta-sheets in electron-density maps

A method for rapidly building beta-sheets into electron-density maps is presented. beta-Strands are identified as tubes of high density adjacent to and nearly parallel to other tubes of density. The alignment and direction of each strand are identified from the pattern of high density corresponding to carbonyl and C(beta) atoms along the strand averaged over all repeats present in the strand. The beta-strands obtained are then assembled into a single atomic model of the beta-sheet regions. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. The beta-sheet regions were nearly completely built in all but two cases, the exceptions being one structure at 2.5 A resolution in which a third of the residues in beta-sheets were built and a structure at 3.8 A in which under 10% were built. The overall average r.m.s.d. of main-chain atoms in the residues built using this method compared with refined models of the structures was 1.5 A.

Rapid chain tracing of polypeptide backbones in electron-density maps

A method for rapid chain tracing of polypeptide backbones at moderate resolution is presented.

A method for the rapid tracing of polypeptide backbones has been developed. The method creates an approximate chain tracing that is useful for visual evaluation of whether a structure has been solved and for use in scoring the quality of electron-density maps. The essence of the method is to (i) sample candidate C^α positions at spacings of approximately 0.6 Å along ridgelines of high electron density, (ii) list all possible nonapeptides that satisfy simple geometric and density criteria using these candidate C^α positions, (iii) score the nonapeptides and choose the highest scoring ones, and (iv) find the longest chains that can be made by connecting nonamers. An indexing and storage scheme that allows a single calculation of most distances and density values is used to speed up the process. The method was applied to 42 density-modified electron-density maps at resolutions from 1.5 to 3.8 Å. A total of 21 428 residues in these maps were traced in 24 CPU min with an overall r.m.s.d. of 1.61 Å for C^α atoms compared with the known refined structures. The method appears to be suitable for rapid evaluation of electron-density map quality.

Decision-making in structure solution using Bayesian estimates of map quality: the PHENIX AutoSol wizard

Ten measures of experimental electron-density-map quality are examined and the skewness of electron density is found to be the best indicator of actual map quality. A Bayesian approach to estimating map quality is developed and used in the PHENIX AutoSol wizard to make decisions during automated structure solution.

Estimates of the quality of experimental maps are important in many stages of structure determination of macromolecules. Map quality is defined here as the correlation between a map and the corresponding map obtained using phases from the final refined model. Here, ten different measures of experimental map quality were examined using a set of 1359 maps calculated by re-analysis of 246 solved MAD, SAD and MIR data sets. A simple Bayesian approach to estimation of map quality from one or more measures is presented. It was found that a Bayesian estimator based on the skewness of the density values in an electron-density map is the most accurate of the ten individual Bayesian estimators of map quality examined, with a correlation between estimated and actual map quality of 0.90. A combination of the skewness of electron density with the local correlation of r.m.s. density gives a further improvement in estimating map quality, with an overall correlation coefficient of 0.92. The PHENIX AutoSol wizard carries out automated structure solution based on any combination of SAD, MAD, SIR or MIR data sets. The wizard is based on tools from the PHENIX package and uses the Bayesian estimates of map quality described here to choose the highest quality solutions after experimental phasing.

The YebC family protein PA0964 negatively regulates the Pseudomonas aeruginosa quinolone signal system and pyocyanin production

Bacterial pathogenicity is often manifested by the expression of various cell-associated and secreted virulence factors, such as exoenzymes, protease, and toxins. In Pseudomonas aeruginosa, the expression of virulence genes is coordinately controlled by the global regulatory quorum-sensing systems, which includes the las and rhl systems as well as the Pseudomonas quinolone signal (PQS) system. Phenazine compounds are among the virulence factors under the control of both the rhl and PQS systems. In this study, regulation of the phzA1B1C1D1E1 (phzA1) operon, which is involved in phenazine synthesis, was investigated. In an initial study of inducing conditions, we observed that phzA1 was induced by subinhibitory concentrations of tetracycline. Screening of 13,000 mutants revealed 32 genes that altered phzA1 expression in the presence of subinhibitory tetracycline concentrations. Among them, the gene PA0964, designated pmpR (pqsR-mediated PQS regulator), has been identified as a novel regulator of the PQS system. It belongs to a large group of widespread conserved hypothetical proteins with unknown function, the YebC protein family (Pfam family DUF28). It negatively regulates the quorum-sensing response regulator pqsR of the PQS system by binding at its promoter region. Alongside phzA1 expression and phenazine and pyocyanin production, a set of virulence factors genes controlled by both rhl and the PQS were shown to be modulated by PmpR. Swarming motility and biofilm formation were also significantly affected. The results added another layer of regulation in the rather complex quorum-sensing systems in P. aeruginosa and demonstrated a clear functional clue for the YebC family proteins.

Preliminary X-ray analysis of XC5848, a hypothetical ORFan protein with an Sm-like motif from Xanthomonas campestris

XC5848, a hypothetical protein from the pathogenic bacterium Xanthomonas campestris that causes black rot, has been chosen as a potential target for the discovery of novel folds. It is unique to the Xanthomonas genus and has significant sequence identity mainly to corresponding proteins from the Xanthomonas genus. In this paper, the cloning, overexpression, purification and crystallization of the XC5848 protein are reported. The XC5848 crystals diffracted to a resolution of at least 1.68 A. They belong to the orthorhombic space group P2(1)2(1)2(1), with unit-cell parameters a = 48.13, b = 51.62, c = 82.32 A. Two molecules were found in each asymmetric unit. Preliminary structural studies nevertheless indicate that XC5848 belongs to the highly conserved Sm-like alpha-beta-beta-beta-beta fold. However, significant differences in sequence and structure were observed. It therefore represents a novel variant of the crucial Sm-like motif that is heavily involved in mRNA splicing and degradation.

The crystal structure of Aq_328 from the hyperthermophilic bacteria Aquifex aeolicus shows an ancestral histone fold

The structure of Aq_328, an uncharacterized protein from hyperthermophilic bacteria Aquifex aeolicus, has been determined to 1.9 A by using multi-wavelength anomalous diffraction (MAD) phasing. Although the amino acid sequence analysis shows that Aq_328 has no significant similarity to proteins with a known structure and function, the structure comparison by using the Dali server reveals that it: (1) assumes a histone-like fold, and (2) is similar to an ancestral nuclear histone protein (PDB code 1F1E) with z-score 8.1 and RMSD 3.6 A over 124 residues. A sedimentation equilibrium experiment indicates that Aq_328 is a monomer in solution, with an average sedimentation coefficient of 2.4 and an apparent molecular weight of about 20 kDa. The overall architecture of Aq_328 consists of two noncanonical histone domains in tandem repeat within a single chain, and is similar to eukaryotic heterodimer (H2A/H2B and H3/H4) and an archaeal histone heterodimer (HMfA/HMfB). The sequence comparisons between the two histone domains of Aq_328 and six eukaryotic/archaeal histones demonstrate that most of the conserved residues that underlie the Aq_328 architecture are used to build and stabilize the two cross-shaped antiparallel histone domains. The high percentage of salt bridges in the structure could be a factor in the protein's thermostability. The structural similarities to other histone-like proteins, molecular properties, and potential function of Aq_328 are discussed in this paper.

Preparation, crystallization and preliminary X-ray characterization of a conserved hypothetical protein XC1692 from Xanthomonas campestris

Xanthomonas campestris pv. campestris strain 17 is a Gram-negative yellow-pigmented pathogenic bacterium that causes black rot, one of the major worldwide diseases of cruciferous crops. Its genome contains approximately 4500 genes, one third of which have no known structure and/or function yet are highly conserved among several different bacterial genuses. One of these gene products is XC1692 protein, containing 141 amino acids. It was overexpressed in Escherichia coli, purified and crystallized in a variety of forms using the hanging-drop vapour-diffusion method. The crystals diffract to at least 1.45 A resolution. They are hexagonal and belong to space group P6(3), with unit-cell parameters a = b = 56.9, c = 71.0 A. They contain one molecule per asymmetric unit.

Cloning, purification, crystallization and preliminary X-ray analysis of XC229, a conserved hypothetical protein from Xanthomonas campestris

Xanthomonas campestris pv. campestris is a Gram-negative yellow-pigmented pathogenic bacterium that causes black rot, one of the major worldwide diseases of cruciferous crops. Its genome contains approximately 4500 genes, roughly one third of which have no known structure and/or function. However, some of these unknown genes are highly conserved among several different bacterial genuses. XC229 is one such protein containing 134 amino acids. It was overexpressed in Escherichia coli, purified and crystallized using the hanging-drop vapour-diffusion method. The crystal diffracted to a resolution of at least 1.80 A. It is cubic and belongs to space group I2(x)3, with unit-cell parameters a = b = c = 106.8 A. It contains one or two molecules per asymmetric unit.

Functional clues for hypothetical proteins based on genomic context analysis in prokaryotes

Three integrated genomic context methods were used to annotate uncharacterized proteins in 102 bacterial genomes. Of 7853 orthologous groups with unknown function containing 45,110 proteins, 1738 groups could be linked to functionally associated partners. In many cases, those partners are uncharacterized themselves (hinting at newly identified modules) or have been described in general terms only. However, we were able to assign pathways, cellular processes or physical complexes for 273 groups (encompassing 3624 previously functionally uncharacterized proteins).

'Conserved hypothetical' proteins: prioritization of targets for experimental study

Comparative genomics shows that a substantial fraction of the genes in sequenced genomes encodes 'conserved hypothetical' proteins, i.e. those that are found in organisms from several phylogenetic lineages but have not been functionally characterized. Here, we briefly discuss recent progress in functional characterization of prokaryotic 'conserved hypothetical' proteins and the possible criteria for prioritizing targets for experimental study. Based on these criteria, the chief one being wide phyletic spread, we offer two 'top 10' lists of highly attractive targets. The first list consists of proteins for which biochemical activity could be predicted with reasonable confidence but the biological function was predicted only in general terms, if at all ('known unknowns'). The second list includes proteins for which there is no prediction of biochemical activity, even if, for some, general biological clues exist ('unknown unknowns'). The experimental characterization of these and other 'conserved hypothetical' proteins is expected to reveal new, crucial aspects of microbial biology and could also lead to better functional prediction for medically relevant human homologs.

Structure-based functional inference in structural genomics

The dramatically increasing number of new protein sequences arising from genomics and proteomics requires the need for methods to rapidly and reliably infer the molecular and cellular functions of these proteins. One such approach, structural genomics, aims to delineate the total repertoire of protein folds in nature, thereby providing three-dimensional folding patterns for all proteins and to infer molecular functions of the proteins based on the combined information of structures and sequences. The goal of obtaining protein structures on a genomic scale has motivated the development of high throughput technologies and protocols for macromolecular structure determination that have begun to produce structures at a greater rate than previously possible. These new structures have revealed many unexpected functional inferences and evolutionary relationships that were hidden at the sequence level. Here, we present samples of structures determined at Berkeley Structural Genomics Center and collaborators' laboratories to illustrate how structural information provides and complements sequence information to deduce the functional inferences of proteins with unknown molecular functions. Two of the major premises of structural genomics are to discover a complete repertoire of protein folds in nature and to find molecular functions of the proteins whose functions are not predicted from sequence comparison alone. To achieve these objectives on a genomic scale, new methods, protocols, and technologies need to be developed by multi-institutional collaborations worldwide. As part of this effort, the Protein Structure Initiative has been launched in the United States (PSI; www.nigms.nih.gov/funding/psi.html). Although infrastructure building and technology development are still the main focus of structural genomics programs, a considerable number of protein structures have already been produced, some of them coming directly out of semiautomated structure determination pipelines. The Berkeley Structural Genomics Center (BSGC) has focused on the proteins of Mycoplasma or their homologues from other organisms as its structural genomics targets because of the minimal genome size of the Mycoplasmas as well as their relevance to human and animal pathogenicity (http://www.strgen.org). Here we present several protein examples encompassing a spectrum of functional inferences obtainable from their three-dimensional structures in five situations, where the inferences are new and testable, and are not predictable from protein sequence information alone.

Identification and functional analysis of 'hypothetical' genes expressed in Haemophilus influenzae

The progress in genome sequencing has led to a rapid accumulation in GenBank submissions of uncharacterized 'hypothetical' genes. These genes, which have not been experimentally characterized and whose functions cannot be deduced from simple sequence comparisons alone, now comprise a significant fraction of the public databases. Expression analyses of Haemophilus influenzae cells using a combination of transcriptomic and proteomic approaches resulted in confident identification of 54 'hypothetical' genes that were expressed in cells under normal growth conditions. In an attempt to understand the functions of these proteins, we used a variety of publicly available analysis tools. Close homologs in other species were detected for each of the 54 'hypothetical' genes. For 16 of them, exact functional assignments could be found in one or more public databases. Additionally, we were able to suggest general functional characterization for 27 more genes (comprising approximately 80% total). Findings from this analysis include the identification of a pyruvate-formate lyase-like operon, likely to be expressed not only in H.influenzae but also in several other bacteria. Further, we also observed three genes that are likely to participate in the transport and/or metabolism of sialic acid, an important component of the H.influenzae lipo-oligosaccharide. Accurate functional annotation of uncharacterized genes calls for an integrative approach, combining expression studies with extensive computational analysis and curation, followed by eventual experimental verification of the computational predictions.

What we have learned about prokaryotes from structural genomics

Five years ago systematic determination and theoretical analysis of all protein structures encoded in model prokaryotic organisms was proposed as a powerful way to obtain new insights into protein function and the variety of protein folds. What has been the pay-off from studying structures in genomic context? Have we learned anything new about protein structure? Can we now predict protein function better? In this contribution, I summarize the status of large-scale structure determination projects on prokaryotes and provide an overview of the main results obtained from experimental and theoretical studies in this dynamic research field.

Protein informatics towards function identification

The study of structural genomics and structural proteomics has determined the tertiary structures of many hypothetical proteins, whose molecular functions could not be understood using conventional methods. In order to infer the geometrical location of the functional site, the biochemical function and the biological function of the hypothetical protein, much effort has been made in protein informatics. The importance of heterogeneous databases and various descriptors of amino acid sequences, tertiary structures and pathways on the proteome scale has been emphasised.

Prediction of catalytic residues in enzymes based on known tertiary structure, stability profile, and sequence conservation

The catalytic or functionally important residues of a protein are known to exist in evolutionarily constrained regions. However, the patterns of residue conservation alone are sometimes not very informative, depending on the homologous sequences available for a given query protein. Here, we present an integrated method to locate the catalytic residues in an enzyme from its sequence and structure. Mutations of functional residues usually decrease the activity, but concurrently often increase stability. Also, catalytic residues tend to occupy partially buried sites in holes or clefts on the molecular surface. After confirming these general tendencies by carrying out statistical analyses on 49 representative enzymes, these data together with amino acid conservation were evaluated. This novel method exhibited better sensitivity in the prediction accuracy than traditional methods that consider only the residue conservation. We applied it to some so-called "hypothetical" proteins, with known structures but undefined functions. The relationships among the catalytic, conserved, and destabilizing residues in enzymatic proteins are discussed.