Journal of structural and functional genomics最新文献

We present a new method for predicting protein-ligand-binding sites based on protein three-dimensional structure and amino acid conservation. This method involves calculation of the van der Waals interaction energy between a protein and many probes placed on the protein surface and subsequent clustering of the probes with low interaction energies to identify the most energetically favorable locus. In addition, it uses amino acid conservation among homologous proteins. Ligand-binding sites were predicted by combining the interaction energy and the amino acid conservation score. The performance of our prediction method was evaluated using a non-redundant dataset of 348 ligand-bound and ligand-unbound protein structure pairs, constructed by filtering entries in a ligand-binding site structure database, LigASite. Ligand-bound structure prediction (bound prediction) indicated that 74.0 % of predicted ligand-binding sites overlapped with real ligand-binding sites by over 25 % of their volume. Ligand-unbound structure prediction (unbound prediction) indicated that 73.9 % of predicted ligand-binding residues overlapped with real ligand-binding residues. The amino acid conservation score improved the average prediction accuracy by 17.0 and 17.6 points for the bound and unbound predictions, respectively. These results demonstrate the effectiveness of the combined use of the interaction energy and amino acid conservation in the ligand-binding site prediction.

ProMOL, a plugin for the PyMOL molecular graphics system, is a structure-based protein function prediction tool. ProMOL includes a set of routines for building motif templates that are used for screening query structures for enzyme active sites. Previously, each motif template was generated manually and required supervision in the optimization of parameters for sensitivity and selectivity. We developed an algorithm and workflow for the automation of motif building and testing routines in ProMOL. The algorithm uses a set of empirically derived parameters for optimization and requires little user intervention. The automated motif generation algorithm was first tested in a performance comparison with a set of manually generated motifs based on identical active sites from the same 112 PDB entries. The two sets of motifs were equally effective in identifying alignments with homologs and in rejecting alignments with unrelated structures. A second set of 296 active site motifs were generated automatically, based on Catalytic Site Atlas entries with literature citations, as an expansion of the library of existing manually generated motif templates. The new motif templates exhibited comparable performance to the existing ones in terms of hit rates against native structures, homologs with the same EC and Pfam designations, and randomly selected unrelated structures with a different EC designation at the first EC digit, as well as in terms of RMSD values obtained from local structural alignments of motifs and query structures. This research is supported by NIH grant GM078077.

Multiprotein complexes play essential roles in all cells and X-ray crystallography can provide unparalleled insight into their structure and function. Many of these complexes are believed to be sufficiently stable for structural biology studies, but the production of protein-protein complexes using recombinant technologies is still labor-intensive. We have explored several strategies for the identification and cloning of heterodimers and heterotrimers that are compatible with the high-throughput (HTP) structural biology pipeline developed for single proteins. Two approaches are presented and compared which resulted in co-expression of paired genes from a single expression vector. Native operons encoding predicted interacting proteins were selected from a repertoire of genomes, and cloned directly to expression vector. In an alternative approach, Helicobacter pylori proteins predicted to interact strongly were cloned, each associated with translational control elements, then linked into an artificial operon. Proteins were then expressed and purified by standard HTP protocols, resulting to date in the structure determination of two H. pylori complexes.

The MazG family proteins, which are highly conserved in bacteria, are nucleoside triphosphate pyrophosphohydrolases that hydrolyze all canonical nucleoside triphosphates, and are also involved in removing noncanonical nucleoside triphosphates to prevent their incorporation into DNA or RNA. The primary structure of TM0360 from Thermotoga maritima MSB8 suggested that TM0360 is a MazG-related nucleoside triphosphate pyrophosphohydrolase. The crystal structure of the TM0360 protein was determined by the MAD technique at 2.0 Å resolution. The asymmetric unit contains an intact dimer molecule. The overall structure of TM0360 is similar to the known structures of the dimeric MazG protein and dUTPases. The putative NTP binding pocket in TM0360, identified by considering the probable NTP-interacting residues and structural features, suggested that TM0360 resembles the C-terminal domain of Escherichia coli MazG, although TM0360 may be a truncated paralog of the N-terminal domain of T. maritima MazG (TM0913), according to its primary structure. The putative function of TM0360 is discussed, based on structural homology.

ZFAT is a transcriptional regulator, containing eighteen C2H2-type zinc-fingers and one AT-hook, involved in autoimmune thyroid disease, apoptosis, and immune-related cell survival. We determined the solution structures of the thirteen individual ZFAT zinc-fingers (ZF) and the tandemly arrayed zinc-fingers in the regions from ZF2 to ZF5, by NMR spectroscopy. ZFAT has eight uncommon bulged-out helix-containing zinc-fingers, and six of their structures (ZF4, ZF5, ZF6, ZF10, ZF11, and ZF13) were determined. The distribution patterns of the putative DNA-binding surface residues are different among the ZFAT zinc-fingers, suggesting the distinct DNA sequence preferences of the N-terminal and C-terminal zinc-fingers. Since ZFAT has three to five consecutive tandem zinc-fingers, which may cooperatively function as a unit, we also determined two tandemly arrayed zinc-finger structures, between ZF2 to ZF4 and ZF3 to ZF5. Our NMR spectroscopic analysis detected the interaction between ZF4 and ZF5, which are connected by an uncommon linker sequence, KKIK. The ZF4-ZF5 linker restrained the relative structural space between the two zinc-fingers in solution, unlike the other linker regions with determined structures, suggesting the involvement of the ZF4-ZF5 interfinger linker in the regulation of ZFAT function.

Vectors designed for protein production in Escherichia coli and by wheat germ cell-free translation were tested using 21 well-characterized eukaryotic proteins chosen to serve as controls within the context of a structural genomics pipeline. The controls were carried through cloning, small-scale expression trials, large-scale growth or synthesis, and purification. Successfully purified proteins were also subjected to either crystallization trials or (1)H-(15)N HSQC NMR analyses. Experiments evaluated: (1) the relative efficacy of restriction/ligation and recombinational cloning systems; (2) the value of maltose-binding protein (MBP) as a solubility enhancement tag; (3) the consequences of in vivo proteolysis of the MBP fusion as an alternative to post-purification proteolysis; (4) the effect of the level of LacI repressor on the yields of protein obtained from E. coli using autoinduction; (5) the consequences of removing the His tag from proteins produced by the cell-free system; and (6) the comparative performance of E. coli cells or wheat germ cell-free translation. Optimal promoter/repressor and fusion tag configurations for each expression system are discussed.

The methylmalonyl Co-A mutase-associated GTPase MeaB from Methylobacterium extorquens is involved in glyoxylate regulation and required for growth. In humans, mutations in the homolog methylmalonic aciduria associated protein (MMAA) cause methylmalonic aciduria, which is often fatal. The central role of MeaB from bacteria to humans suggests that MeaB is also important in other, pathogenic bacteria such as Mycobacterium tuberculosis. However, the identity of the mycobacterial MeaB homolog is presently unclear. Here, we identify the M. tuberculosis protein Rv1496 and its homologs in M. smegmatis and M. thermoresistibile as MeaB. The crystal structures of all three homologs are highly similar to MeaB and MMAA structures and reveal a characteristic three-domain homodimer with GDP bound in the G domain active site. A structure of Rv1496 obtained from a crystal grown in the presence of GTP exhibited electron density for GDP, suggesting GTPase activity. These structures identify the mycobacterial MeaB and provide a structural framework for therapeutic targeting of M. tuberculosis MeaB.

Working with a combination of ProMOL (a plugin for PyMOL that searches a library of enzymatic motifs for local structural homologs), BLAST and Pfam (servers that identify global sequence homologs), and Dali (a server that identifies global structural homologs), we have begun the process of assigning functional annotations to the approximately 3,500 structures in the Protein Data Bank that are currently classified as having "unknown function". Using a limited template library of 388 motifs, over 500 promising in silico matches have been identified by ProMOL, among which 65 exceptionally good matches have been identified. The characteristics of the exceptionally good matches are discussed.