Searching protein structure databases with DaliLite v.3

Advances and pitfalls of protein structural alignment

Structure comparison opens a window into the distant past of protein evolution, which has been unreachable by sequence comparison alone. With 55,000 entries in the Protein Data Bank and about 500 new structures added each week, automated processing, comparison, and classification are necessary. A variety of methods use different representations, scoring functions, and optimization algorithms, and they generate contradictory results even for moderately distant structures. Sequence mutations, insertions, and deletions are accommodated by plastic deformations of the common core, retaining the precise geometry of the active site, and peripheral regions may refold completely. Therefore structure comparison methods that allow for flexibility and plasticity generate the most biologically meaningful alignments. Active research directions include both the search for fold invariant features and the modeling of structural transitions in evolution. Advances have been made in algorithmic robustness, multiple alignment, and speeding up database searches.

The SALAMI protein structure search server

Protein structures often show similarities to another which would not be seen at the sequence level. Given the coordinates of a protein chain, the SALAMI server atwww.zbh.uni-hamburg.de/salami will search the protein data bank and return a set of similar structures without using sequence information. The results page lists the related proteins, details of the sequence and structure similarity and implied sequence alignments. Via a simple structure viewer, one can view superpositions of query and library structures and finally download superimposed coordinates. The alignment method is very tolerant of large gaps and insertions, and tends to produce slightly longer alignments than other similar programs.

COPS--a novel workbench for explorations in fold space

The COPS (Classification Of Protein Structures) web server provides access to the complete repertoire of known protein structures and protein structural domains. The COPS classification encodes pairwise structural similarities as quantified metric relationships. The resulting metrical structure is mapped to a hierarchical tree, which is largely equivalent to the structure of a file browser. Exploiting this relationship we implemented the Fold Space Navigator, a tool that makes navigation in fold space as convenient as browsing through a file system. Moreover, pairwise structural similarities among the domains can be visualized and inspected instantaneously. COPS is updated weekly and stays concurrent with the PDB repository. The server also exposes the COPS classification pipeline. Newly determined structures uploaded to the server are chopped into domains, the locations of the new domains in the classification tree are determined, and their neighborhood can be immediately explored through the Fold Space Navigator. The COPS web server is accessible at http://cops.services.came.sbg.ac.at/.

Crystal structure and association behaviour of the GluR2 amino-terminal domain

Fast excitatory neurotransmission is mediated largely by ionotropic glutamate receptors (iGluRs), tetrameric, ligand-gated ion channel proteins comprised of three subfamilies, AMPA, kainate and NMDA receptors, with each subfamily sharing a common, modular-domain architecture. For all receptor subfamilies, active channels are exclusively formed by assemblages of subunits within the same subfamily, a molecular process principally encoded by the amino-terminal domain (ATD). However, the molecular basis by which the ATD guides subfamily-specific receptor assembly is not known. Here we show that AMPA receptor GluR1- and GluR2-ATDs form tightly associated dimers and, by the analysis of crystal structures of the GluR2-ATD, propose mechanisms by which the ATD guides subfamily-specific receptor assembly.

Crystal structure of the membrane-bound bifunctional transglycosylase PBP1b from Escherichia coli

Drug-resistant bacteria have caused serious medical problems in recent years, and the need for new antibacterial agents is undisputed. Transglycosylase, a multidomain membrane protein essential for cell wall synthesis, is an excellent target for the development of new antibiotics. Here, we determined the X-ray crystal structure of the bifunctional transglycosylase penicillin-binding protein 1b (PBP1b) from Escherichia coli in complex with its inhibitor moenomycin to 2.16-A resolution. In addition to the transglycosylase and transpeptidase domains, our structure provides a complete visualization of this important antibacterial target, and reveals a domain for protein-protein interaction and a transmembrane helix domain essential for substrate binding, enzymatic activity, and membrane orientation.

Crystal structure of the EutL shell protein of the ethanolamine ammonia lyase microcompartment

Bacterial microcompartments (BMCs) are specialized organelles that use proteinaceous membranes to confine chemical reaction spaces. The ethanolamine ammonialyase microcompartment of Escherichia coli represents such a class of cytosolic organelles that enables bacteria to survive on small organic molecules such as ethanolamine as the sole source for carbon and nitrogen. We present here the crystal structure of the shell protein EutL at 2.2-A resolution. With 219 residues, it is the largest representative of this BMC's shell proteins. In the crystal, EutL forms a trimer that exhibits a hexagonally shaped tile structure. The tiles arrange into a tightly packed 2D array that is likely to resemble the proteinaceous membrane of the intact BMC. In contrast to other BMC shell proteins, which have only 1 pore per tile, EutL exhibits 3 pores per tile, thereby significantly increasing the overall porosity of this protein membrane. Each of the individual pores is lined with negatively charged residues and aromatic residues that are proposed to facilitate passive transport of specific solutes. The characteristic shape of the hexagonal tile, which is also found in the microcompartments of carbon-fixating bacteria, may present an inherent and fundamental building unit that may provide a general explanation for the formation of differently sized microcompartments.

Crystal structures of Cif from bacterial pathogens Photorhabdus luminescens and Burkholderia pseudomallei

A pre-requisite for bacterial pathogenesis is the successful interaction of a pathogen with a host. One mechanism used by a broad range of Gram negative bacterial pathogens is to deliver effector proteins directly into host cells through a dedicated type III secretion system where they modulate host cell function. The cycle inhibiting factor (Cif) family of effector proteins, identified in a growing number of pathogens that harbour functional type III secretion systems and have a wide host range, arrest the eukaryotic cell cycle. Here, the crystal structures of Cifs from the insect pathogen/nematode symbiont Photorhabdus luminescens (a gamma-proteobacterium) and human pathogen Burkholderia pseudomallei (a beta-proteobacterium) are presented. Both of these proteins adopt an overall fold similar to the papain sub-family of cysteine proteases, as originally identified in the structure of a truncated form of Cif from Enteropathogenic E. coli (EPEC), despite sharing only limited sequence identity. The structure of an N-terminal region, referred to here as the 'tail-domain' (absent in the EPEC Cif structure), suggests a surface likely to be involved in host-cell substrate recognition. The conformation of the Cys-His-Gln catalytic triad is retained, and the essential cysteine is exposed to solvent and addressable by small molecule reagents. These structures and biochemical work contribute to the rapidly expanding literature on Cifs, and direct further studies to better understand the molecular details of the activity of these proteins.

Crystallographic insights into the autocatalytic assembly mechanism of a bacteriophage tail spike

The tailed bacteriophage phi29 has 12 "appendages" (gene product 12, gp12) attached to its neck region that participate in host cell recognition and entry. In the cell, monomeric gp12 undergoes proteolytic processing that releases the C-terminal domain during assembly into trimers. We report here crystal structures of the protein before and after catalytic processing and show that the C-terminal domain of gp12 is an "autochaperone" that aids trimerization. We also show that autocleavage of the C-terminal domain is a posttrimerization event that is followed by a unique ATP-dependent release. The posttranslationally modified N-terminal part has three domains that function to attach the appendages to the phage, digest the cell wall teichoic acids, and bind irreversibly to the host, respectively. Structural and sequence comparisons suggest that some eukaryotic and bacterial viruses as well as bacterial adhesins might have a similar maturation mechanism as is performed by phi29 gp12 for Bacillus subtilis.

ProSMoS server: a pattern-based search using interaction matrix representation of protein structures

Assessing structural similarity and defining common regions through comparison of protein spatial structures is an important task in functional and evolutionary studies of proteins. There are many servers that compare structures and define sub-structures in common between proteins through superposition and closeness of either coordinates or contacts. However, a natural way to analyze a structure for experts working on structure classification is to look for specific three-dimensional (3D) motifs and patterns instead of finding common features in two proteins. Such motifs can be described by the architecture and topology of major secondary structural elements (SSEs) without consideration of subtle differences in 3D coordinates. Despite the importance of motif-based structure searches, currently there is a shortage of servers to perform this task. Widely known TOPS does not fully address this problem, as it finds only topological match but does not take into account other important spatial properties, such as interactions and chirality. Here, we implemented our approach to protein structure pattern search (ProSMoS) as a web-server. ProSMoS converts 3D structure into an interaction matrix representation including the SSE types, handednesses of connections between SSEs, coordinates of SSE starts and ends, types of interactions between SSEs and beta-sheet definitions. For a user-defined structure pattern, ProSMoS lists all structures from a database that contain this pattern. ProSMoS server will be of interest to structural biologists who would like to analyze very general and distant structural similarities. The ProSMoS web server is available at: http://prodata.swmed.edu/ProSMoS/.

HorA web server to infer homology between proteins using sequence and structural similarity

The biological properties of proteins are often gleaned through comparative analysis of evolutionary relatives. Although protein structure similarity search methods detect more distant homologs than purely sequence-based methods, structural resemblance can result from either homology (common ancestry) or analogy (similarity without common ancestry). While many existing web servers detect structural neighbors, they do not explicitly address the question of homology versus analogy. Here, we present a web server named HorA (Homology or Analogy) that identifies likely homologs for a query protein structure. Unlike other servers, HorA combines sequence information from state-of-the-art profile methods with structure information from spatial similarity measures using an advanced computational technique. HorA aims to identify biologically meaningful connections rather than purely 3D-geometric similarities. The HorA method finds approximately 90% of remote homologs defined in the manually curated database SCOP. HorA will be especially useful for finding remote homologs that might be overlooked by other sequence or structural similarity search servers. The HorA server is available at http://prodata.swmed.edu/horaserver.

Atomic structure of a folate/FAD-dependent tRNA T54 methyltransferase

tRNAs from all 3 phylogenetic domains have a 5-methyluridine at position 54 (T54) in the T-loop. The methyl group is transferred from S-adenosylmethionine by TrmA methyltransferase in most Gram-negative bacteria and some archaea and eukaryotes, whereas it is transferred from 5,10-methylenetetrahydrofolate (MTHF) by TrmFO, a folate/FAD-dependent methyltransferase, in most Gram-positive bacteria and some Gram-negative bacteria. However, the catalytic mechanism remains unclear, because the crystal structure of TrmFO has not been solved. Here, we report the crystal structures of Thermus thermophilus TrmFO in its free form, tetrahydrofolate (THF)-bound form, and glutathione-bound form at 2.1-, 1.6-, and 1.05-A resolutions, respectively. TrmFO consists of an FAD-binding domain and an insertion domain, which both share structural similarity with those of GidA, an enzyme involved in the 5-carboxymethylaminomethylation of U34 of some tRNAs. However, the overall structures of TrmFO and GidA are basically different because of their distinct domain orientations, which are consistent with their respective functional specificities. In the THF complex, the pteridin ring of THF is sandwiched between the flavin ring of FAD and the imidazole ring of a His residue. This structure provides a snapshot of the folate/FAD-dependent methyl transfer, suggesting that the transferring methylene group of MTHF is located close to the redox-active N5 atom of FAD. Furthermore, we established an in vitro system to measure the methylation activity. Our TrmFO-tRNA docking model, in combination with mutational analyses, suggests a catalytic mechanism, in which the methylene of MTHF is directly transferred onto U54, and then the exocyclic methylene of U54 is reduced by FADH(2).

The solution structure of the first PHD finger of autoimmune regulator in complex with non-modified histone H3 tail reveals the antagonistic role of H3R2 methylation

Plant homeodomain (PHD) fingers are often present in chromatin-binding proteins and have been shown to bind histone H3 N-terminal tails. Mutations in the autoimmune regulator (AIRE) protein, which harbours two PHD fingers, cause a rare monogenic disease, autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED). AIRE activates the expression of tissue-specific antigens by directly binding through its first PHD finger (AIRE-PHD1) to histone H3 tails non-methylated at K4 (H3K4me0). Here, we present the solution structure of AIRE-PHD1 in complex with H3K4me0 peptide and show that AIRE-PHD1 is a highly specialized non-modified histone H3 tail reader, as post-translational modifications of the first 10 histone H3 residues reduce binding affinity. In particular, H3R2 dimethylation abrogates AIRE-PHD1 binding in vitro and reduces the in vivo activation of AIRE target genes in HEK293 cells. The observed antagonism by R2 methylation on AIRE-PHD1 binding is unique among the H3K4me0 histone readers and represents the first case of epigenetic negative cross-talk between non-methylated H3K4 and methylated H3R2. Collectively, our results point to a very specific histone code responsible for non-modified H3 tail recognition by AIRE-PHD1 and describe at atomic level one crucial step in the molecular mechanism responsible for antigen expression in the thymus.

Structural basis for m7G-cap hypermethylation of small nuclear, small nucleolar and telomerase RNA by the dimethyltransferase TGS1

The 5′-cap of spliceosomal small nuclear RNAs, some small nucleolar RNAs and of telomerase RNA was found to be hypermethylated in vivo. The Trimethylguanosine Synthase 1 (TGS1) mediates this conversion of the 7-methylguanosine-cap to the 2,2,7-trimethylguanosine (m₃G)-cap during maturation of the RNPs. For mammalian UsnRNAs the generated m^2,2,7G-cap is one part of a bipartite import signal mediating the transport of the UsnRNP-core complex into the nucleus. In order to understand the structural organization of human TGS1 as well as substrate binding and recognition we solved the crystal structure of the active TGS1 methyltransferase domain containing both, the minimal substrate m⁷GTP and the reaction product S-adenosyl-l-homocysteine (AdoHcy). The methyltransferase of human TGS1 harbors the canonical class 1 methyltransferase fold as well as an unique N-terminal, α-helical domain of 40 amino acids, which is essential for m⁷G-cap binding and catalysis. The crystal structure of the substrate bound methyltransferase domain as well as mutagenesis studies provide insight into the catalytic mechanism of TGS1.

Characterization of 3-ketosteroid 9{alpha}-hydroxylase, a Rieske oxygenase in the cholesterol degradation pathway of Mycobacterium tuberculosis

KshAB (3-Ketosteroid 9alpha-hydroxylase) is a two-component Rieske oxygenase (RO) in the cholesterol catabolic pathway of Mycobacterium tuberculosis. Although the enzyme has been implicated in pathogenesis, it has largely been characterized by bioinformatics and molecular genetics. Purified KshB, the reductase component, was a monomeric protein containing a plant-type [2Fe-2S] cluster and FAD. KshA, the oxygenase, was a homotrimer containing a Rieske [2Fe-2S] cluster and mononuclear ferrous iron. Of two potential substrates, reconstituted KshAB had twice the specificity for 1,4-androstadiene-3,17-dione as for 4-androstene-3,17-dione. The transformation of both substrates was well coupled to the consumption of O(2). Nevertheless, the reactivity of KshAB with O(2) was low in the presence of 1,4-androstadiene-3,17-dione, with a k(cat)/K(m)(O(2)) of 2450 +/- 80 m(-1) s(-1). The crystallographic structure of KshA, determined to 2.3A(,) revealed an overall fold and a head-to-tail subunit arrangement typical of ROs. The central fold of the catalytic domain lacks all insertions found in characterized ROs, consistent with a minimal and perhaps archetypical RO catalytic domain. The structure of KshA is further distinguished by a C-terminal helix, which stabilizes subunit interactions in the functional trimer. Finally, the substrate-binding pocket extends farther into KshA than in other ROs, consistent with the large steroid substrate, and the funnel accessing the active site is differently orientated. This study provides a solid basis for further studies of a key steroid-transforming enzyme of biotechnological and medical importance.

A structural model for the damage-sensing complex in bacterial nucleotide excision repair

Nucleotide excision repair is distinguished from other DNA repair pathways by its ability to process a wide range of structurally unrelated DNA lesions. In bacteria, damage recognition is achieved by the UvrA.UvrB ensemble. Here, we report the structure of the complex between the interaction domains of UvrA and UvrB. These domains are necessary and sufficient for full-length UvrA and UvrB to associate and thereby form the DNA damage-sensing complex of bacterial nucleotide excision repair. The crystal structure and accompanying biochemical analyses suggest a model for the complete damage-sensing complex.

Protein structure prediction center in CASP8

We present an outline of the Critical Assessment of Protein Structure Prediction (CASP) infrastructure implemented at the University of California, Davis, Protein Structure Prediction Center. The infrastructure supports selection and validation of prediction targets, collection of predictions, standard evaluation of submitted predictions, and presentation of results. The Center also supports information exchange relating to CASP experiments and structure prediction in general. Technical aspects of conducting the CASP8 experiment and relevant statistics are also provided.

The other 90% of the protein: assessment beyond the Calphas for CASP8 template-based and high-accuracy models

For template-based modeling in the CASP8 Critical Assessment of Techniques for Protein Structure Prediction, this work develops and applies six new full-model metrics. They are designed to complement and add value to the traditional template-based assessment by the global distance test (GDT) and related scores (based on multiple superpositions of Calpha atoms between target structure and predictions labeled "Model 1"). The new metrics evaluate each predictor group on each target, using all atoms of their best model with above-average GDT. Two metrics evaluate how "protein-like" the predicted model is: the MolProbity score used for validating experimental structures, and a mainchain reality score using all-atom steric clashes, bond length and angle outliers, and backbone dihedrals. Four other new metrics evaluate match of model to target for mainchain and sidechain hydrogen bonds, sidechain end positioning, and sidechain rotamers. Group-average Z-score across the six full-model measures is averaged with group-average GDT Z-score to produce the overall ranking for full-model, high-accuracy performance. Separate assessments are reported for specific aspects of predictor-group performance, such as robustness of approximately correct template or fold identification, and self-scoring ability at identifying the best of their models. Fold identification is distinct from but correlated with group-average GDT Z-score if target difficulty is taken into account, whereas self-scoring is done best by servers and is uncorrelated with GDT performance. Outstanding individual models on specific targets are identified and discussed. Predictor groups excelled at different aspects, highlighting the diversity of current methodologies. However, good full-model scores correlate robustly with high Calpha accuracy.

Structural analysis of the DNA-binding domain of the Helicobacter pylori response regulator ArsR

The Helicobacter pylori ArsS-ArsR two-component signal transduction system, comprised of a sensor histidine kinase (ArsS) and a response regulator (ArsR), allows the bacteria to regulate gene expression in response to acidic pH. We expressed and purified the full-length ArsR protein and the DNA-binding domain of ArsR (ArsR-DBD), and we analyzed the tertiary structure of the ArsR-DBD using solution nuclear magnetic resonance (NMR) methods. Both the full-length ArsR and the ArsR-DBD behaved as monomers in size exclusion chromatography experiments. The structure of ArsR-DBD consists of an N-terminal four-stranded beta-sheet, a helical core, and a C-terminal beta-hairpin. The overall tertiary fold of the ArsR-DBD is most closely related to DBD structures of the OmpR/PhoB subfamily of bacterial response regulators. However, the orientation of the N-terminal beta-sheet with respect to the rest of the DNA-binding domain is substantially different in ArsR compared with the orientation in related response regulators. Molecular modeling of an ArsR-DBD-DNA complex permits identification of protein elements that are predicted to bind target DNA sequences and thereby regulate gene transcription in H. pylori.

Crystal structure of the N-terminal domain of anaphase-promoting complex subunit 7

Anaphase-promoting complex or cyclosome (APC/C) is an unusual E3 ubiquitin ligase and an essential protein that controls mitotic progression. APC/C includes at least 13 subunits, but no structure has been determined for any tetratricopeptide repeat (TPR)-containing subunit (Apc3 and -6-8) in the TPR subcomplex of APC/C. Apc7 is a TPR-containing subunit that exists only in vertebrate APC/C. Here we report the crystal structure of quad mutant of nApc7 (N-terminal fragment, residues 1-147) of human Apc7 at a resolution of 2.5 A. The structure of nApc7 adopts a TPR-like motif and has a unique dimerization interface, although the protein does not contain the conserved TPR sequence. Based on the structure of nApc7, in addition to previous experimental findings, we proposed a putative homodimeric structure for full-length Apc7. This model suggests that TPR-containing subunits self-associate and bind to adaptors and substrates via an IR peptide in TPR-containing subunits of APC/C.