Functional implications from crystal structures of the conserved Bacillus subtilis protein Maf with and without dUTP

Identification of a new nucleotide binding site by structural alignment and site directed mutagenesis

Nucleotide binding proteins are involved in many important cellular processes and form one of the largest protein families. Traditionally, the identification of nucleotide binding motif, such as the ATP binding P-loop, has relied on the comparison of protein sequences, consideration of the function of each of the proteins and the identification of signature motifs within the sequence. Sometimes, it is difficult to identify nucleotide binding proteins based on sequence alignment because of increased evolutionary distances. In such cases, structural alignments can provide a better guide for comparing specific features of sequences because the overall structures of these motifs are conserved despite low sequence identity. In the present study, on the basis of bioinformatics and structural comparison of three representative protein structures of Ham1 superfamily, YjjX, YggV, and YhdE, previously identified as nucleotide binding proteins, we have identified a novel nucleotide binding motif (T/SXXXXK/R). The importance of this signature motif in binding of nucleotides was validated using site directed mutagenesis. Mutations of conserved residues of the loop either decreased or completely abolished the nucleotide binding activity of the protein. We used the conserved motif identified in the study to search for other proteins having a similar motif. Two proteins, GTP cyclohydrolase II and dephospho-CoA pyrophosphorylase showed presence of the loop, suggesting that this nucleotide binding motif is not unique in the Ham1 superfamily, but represents a novel NTP recognition motif.

A new method to investigate the catalytic mechanism of YhdE pyrophosphatase by using a pyrophosphate fluorescence probe

YhdE is a Maf (multicopy associated filamentation) proteins from Escherichia coli which exhibits pyrophosphatase activity towards selected nucleotides, although its catalytic mechanism remains unclear. Herein we used a novel fluorescence probe (4-isoACBA–Zn(II) complex) to characterize the enzymatic properties of YhdE and its mutant, establishing a new method for assaying pyrophosphatase catalytic function. Our results reveal for the first time that the new fluorescence sensor confers high sensitivity and specificity and pyrophosphate (PPi) is the direct catalytic product of YhdE. Crystal structures of a mutant in the active-site loop (YhdE_E33A) show conformational flexibility implicated in the catalytic mechanism of YhdE. ITC experiments and computational docking further reveal that Asp70 and substrate dTTP coordinate Mn²⁺. Quantum mechanics calculations indicate that YhdE hydrolysis appears to follow a stepwise pathway in which a water molecule first attacks the α-phosphorus atom in the substrate, followed by the release of PPi from the pentavalent intermediate.

Molecular Dynamics Simulation Studies of dTTP Binding and Catalysis Mediated by YhdE Dimerization

YhdE is a Maf-like (multicopy associated filamentation) protein that primarily acts as dTTPase to hydrolyze dTTP into dTMP and two phosphate molecules in cell metabolism pathway. Two crystal structures of YhdE have been previously determined, representing the open and closed active site conformations, respectively. Based on the structures, we have carried out molecular dynamics simulations and free energy calculations to investigate dTTP binding to and hydrolysis by YhdE. Our results suggest that YhdE closed state is structurally more compact than its open state at room temperature. YhdE open state is a favorable conformation for dTTP binding and closed state is a structurally favorable conformation for catalytic reaction. This observation is supported by the structure of YhdE homolog in complex with a nucleotide analog. Free energy calculations reveal that YhdE dimerization occurs preferentially in dTTP binding and is favorable for successive cooperative reaction. The key residues R11, R12 and K80, are found to contribute to the substrate stabilization. Further, YhdE dimerization and binding of dTTP induce the cooperative effect through a direct allosteric communication network in YhdE from the dTTP binding sites in the catalytic center to the intermolecular β-strand in YhdE dimer.

Crystal structure of the MazG-related nucleoside triphosphate pyrophosphohydrolase from Thermotoga maritima MSB8

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.

Insights into the cellular function of YhdE, a nucleotide pyrophosphatase from Escherichia coli

YhdE, a Maf-like protein in Escherichia coli, exhibits nucleotide pyrophosphatase (PPase) activity, yet its cellular function remains unknown. Here, we characterized the PPase activity of YhdE on dTTP, UTP and TTP and determined two crystal structures of YhdE, revealing ‘closed’ and ‘open’ conformations of an adaptive active site. Our functional studies demonstrated that YhdE retards cell growth by prolonging the lag and log phases, particularly under stress conditions. Morphology studies showed that yhdE-knockout cells transformed the normal rod shape of wild-type cells to a more spherical form, and the cell wall appeared to become more flexible. In contrast, YhdE overexpression resulted in filamentous cells. This study reveals the previously unknown involvement of YhdE in cell growth inhibition under stress conditions, cell-division arrest and cell-shape maintenance, highlighting YhdE’s important role in E. coli cell-cycle checkpoints.

Preventive DNA repair by sanitizing the cellular (deoxy)nucleoside triphosphate pool

The occurrence of modified bases in DNA is attributed to some major factors: incorporation of altered nucleotide building blocks and chemical reactions or radiation effects on bases within the DNA structure. Several enzyme families are involved in preventing the incorporation of noncanonical bases playing a 'sanitizing' role. The catalytic mechanism of action of these enzymes has been revealed for a number of representatives in clear structural and kinetic detail. In this review, we focus in detail on those examples where clear evidence has been produced using high-resolution structural studies. Comparing the protein fold and architecture of the enzyme active sites, two main classes of sanitizing deoxyribonucleoside triphosphate pyrophosphatases can be assigned that are distinguished by the site of nucleophilic attack. In enzymes associated with attack at the α-phosphorus, it is shown that coordination of the γ-phosphate group is also ensured by multiple interactions. By contrast, enzymes catalyzing attack at the β-phosphorus atom mainly coordinate the α- and the β-phosphate only. Characteristic differences are also observed with respect to the role of the metal ion cofactor (Mg(2+) ) and the coordination of nucleophilic water. Using different catalytic mechanisms embedded in different protein folds, these enzymes present a clear example of convergent evolution.

Biochemical and structural studies of conserved Maf proteins revealed nucleotide pyrophosphatases with a preference for modified nucleotides

Maf (for multicopy associated filamentation) proteins represent a large family of conserved proteins implicated in cell division arrest but whose biochemical activity remains unknown. Here, we show that the prokaryotic and eukaryotic Maf proteins exhibit nucleotide pyrophosphatase activity against 5-methyl-UTP, pseudo-UTP, 5-methyl-CTP, and 7-methyl-GTP, which represent the most abundant modified bases in all organisms, as well as against canonical nucleotides dTTP, UTP, and CTP. Overexpression of the Maf protein YhdE in E. coli cells increased intracellular levels of dTMP and UMP, confirming that dTTP and UTP are the in vivo substrates of this protein. Crystal structures and site-directed mutagenesis of Maf proteins revealed the determinants of their activity and substrate specificity. Thus, pyrophosphatase activity of Maf proteins toward canonical and modified nucleotides might provide the molecular mechanism for a dual role of these proteins in cell division arrest and house cleaning.

•

Maf proteins represent a family of nucleoside triphosphate pyrophosphatases

•

Maf proteins hydrolyze the canonical nucleotides dTTP, UTP, and CTP

•

Maf proteins are also active against m⁵UTP, m⁵CTP, pseudo-UTP, and m⁷GTP

•

Maf structures reveal the molecular mechanisms of their substrate selectivity

Structural and functional characterization of a noncanonical nucleoside triphosphate pyrophosphatase from Thermotoga maritima

The hyperthermophilic bacterium Thermotoga maritima has a noncanonical nucleoside triphosphatase that catalyzes the conversion of inosine triphosphate (ITP), deoxyinosine triphosphate (dITP) and xanthosine triphosphate (XTP) into inosine monophosphate (IMP), deoxyinosine monophosphate (IMP) and xanthosine monophosphate (XMP), respectively. The k(cat)/K(m) values determined at 323 and 353 K fall between 1.31 × 10(4) and 7.80 × 10(4) M(-1) s(-1). ITP and dITP are slightly preferred over XTP. Activity towards canonical nucleoside triphosphates (ATP and GTP) was not detected. The enzyme has an absolute requirement for Mg(2+) as a cofactor and has a preference for alkaline conditions. A protein X-ray structure of the enzyme with bound IMP was obtained at 2.15 Å resolution. The active site houses a well conserved network of residues that are critical for substrate recognition and catalysis. The crystal structure shows a tetramer with two possible dimer interfaces. One of these interfaces strongly resembles the dimer interface that is found in the structures of other noncanonical nucleoside pyrophosphatases from human (human ITPase) and archaea (Mj0226 and PhNTPase).

Maf acts downstream of ComGA to arrest cell division in competent cells of B. subtilis

Transformable (competent) cells of Bacillus subtilis are blocked in cell division because the traffic ATPase ComGA prevents the formation of FtsZ rings. Although ComGA-deficient cells elongate and form FtsZ rings, cell division remains blocked at a later stage and the cells become mildly filamented. Here we show that the highly conserved protein Maf is synthesized predominantly in competent cells under the direct control of the transcription factor ComK and is solely responsible for the later block in cell division. In vivo and in vitro data show that Maf binds to both ComGA and DivIVA. A point mutation in maf that interferes with Maf binding to DivIVA also interferes with the ability of Maf to inhibit cell division. Based on these findings, we propose that Maf and ComGA mediate mechanisms for the inhibition of cell division in competent cells with Maf acting downstream of ComGA. We further suggest that Maf must interact with DivIVA to inhibit cell division.

Cell division blockage: but this time by a surprisingly conserved protein

Activation of DNA repair proteins is often accompanied by an arrest in cell division. Several proteins have been identified that regulate the division blockage associated with the SOS response. When Bacillus subtilis cells become genetically competent they also activate DNA repair proteins and stop dividing. In this issue of Molecular Microbiology, Briley et al., 2011 describe a new protein involved in this process. This protein, Maf, does not prevent FtsZ polymerization, but it inhibits synthesis of the division septum. What is fascinating about Maf is that it is conserved and can be found in all kingdoms of life.

Structure of a putative NTP pyrophosphohydrolase: YP_001813558.1 from Exiguobacterium sibiricum 255-15

The crystal structure of a putative NTPase, YP_001813558.1 from Exiguobacterium sibiricum 255-15 (PF09934, DUF2166) was determined to 1.78 Å resolution. YP_001813558.1 and its homologs (dimeric dUTPases, MazG proteins and HisE-encoded phosphoribosyl ATP pyrophosphohydrolases) form a superfamily of all-α-helical NTP pyrophosphatases. In dimeric dUTPase-like proteins, a central four-helix bundle forms the active site. However, in YP_001813558.1, an unexpected intertwined swapping of two of the helices that compose the conserved helix bundle results in a `linked dimer' that has not previously been observed for this family. Interestingly, despite this novel mode of dimerization, the metal-binding site for divalent cations, such as magnesium, that are essential for NTPase activity is still conserved. Furthermore, the active-site residues that are involved in sugar binding of the NTPs are also conserved when compared with other α-helical NTPases, but those that recognize the nucleotide bases are not conserved, suggesting a different substrate specificity.

The role of acyl-coenzyme A carboxylase complex in lipstatin biosynthesis of Streptomyces toxytricini

Streptomyces toxytricini produces lipstatin, a specific inhibitor of pancreatic lipase, which is derived from two fatty acid moieties with eight and 14 carbon atoms. The pccB gene locus in 10.6 kb fragment of S. toxytricini chromosomal DNA contains three genes for acyl-coenzyme A carboxylase (ACCase) complex accA3, pccB, and pccE that are presumed to be involved in secondary metabolism. The pccB gene encoding a β subunit of ACCase [carboxyltransferase (CT)] was identified upstream of pccE gene for a small protein of ε subunit. The accA3 encoding the α subunit of ACCase [biotin carboxylase (BC)] was also identified downstream of pccB gene. When the pccB and pccE genes were inactivated by homologous recombination, the lipstatin production was reduced as much as 80%. In contrast, the accumulation of another compound, tetradeca-5.8-dienoic acid (the major lipstatin precursor), was 4.5-fold increased in disruptant compared with wild-type. It implies that PccB of S. toxytricini is involved in the activation of octanoic acid to hexylmalonic acid for lipstatin biosynthesis.

Discrimination between biological interfaces and crystal-packing contacts

A discrimination method between biologically relevant interfaces and artificial crystal-packing contacts in crystal structures was constructed. The method evaluates protein-protein interfaces in terms of complementarities for hydrophobicity, electrostatic potential and shape on the protein surfaces, and chooses the most probable biological interfaces among all possible contacts in the crystal. The method uses a discriminator named as "COMP", which is a linear combination of the complementarities for the above three surface features and does not correlate with the contact area. The discrimination of homo-dimer interfaces from symmetry-related crystal-packing contacts based on the COMP value achieved the modest success rate. Subsequent detailed review of the discrimination results raised the success rate to about 88.8%. In addition, our discrimination method yielded some clues for understanding the interaction patterns in several examples in the PDB. Thus, the COMP discriminator can also be used as an indicator of the "biological-ness" of protein-protein interfaces.

Molecular basis of the antimutagenic activity of the house-cleaning inosine triphosphate pyrophosphatase RdgB from Escherichia coli

Inosine triphosphate pyrophosphatases, which are ubiquitous house-cleaning enzymes, hydrolyze noncanonical nucleoside triphosphates (inosine triphosphate (ITP) and xanthosine triphosphate (XTP)) and prevent the incorporation of hypoxanthine or xanthine into nascent DNA or RNA. Here we present the 1.5-A-resolution crystal structure of the inosine triphosphate pyrophosphatase RdgB from Escherichia coli in a free state and in complex with a substrate (ITP+Ca(2+)) or a product (inosine monophosphate (IMP)). ITP binding to RdgB induced a large displacement of the alpha1 helix, closing the enzyme active site. This positions the conserved Lys13 close to the bridging oxygen between the alpha- and beta-phosphates of the substrate, weakening the P(alpha)-O bond. On the other side of the substrate, the conserved Asp69 is proposed to act as a base coordinating the catalytic water molecule. Our data provide insight into the molecular mechanisms of the substrate selectivity and catalysis of RdgB and other ITPases.

Functional specificity lies within the properties and evolutionary changes of amino acids

The rapid increase in the amount of protein sequence data has created a need for automated identification of sites that determine functional specificity among related subfamilies of proteins. A significant fraction of subfamily specific sites are only marginally conserved, which makes it extremely challenging to detect those amino acid changes that lead to functional diversification. To address this critical problem we developed a method named SPEER (specificity prediction using amino acids' properties, entropy and evolution rate) to distinguish specificity determining sites from others. SPEER encodes the conservation patterns of amino acid types using their physico-chemical properties and the heterogeneity of evolutionary changes between and within the subfamilies. To test the method, we compiled a test set containing 13 protein families with known specificity determining sites. Extensive benchmarking by comparing the performance of SPEER with other specificity site prediction algorithms has shown that it performs better in predicting several categories of subfamily specific sites.

Piecing together the structure–function puzzle: Experiences in structure-based functional annotation of hypothetical proteins

The combination of genomic sequencing with structural genomics has provided a wealth of new structures for previously uncharacterized ORFs, more commonly referred to as hypothetical proteins. This rapid growth has been the direct result of high-throughput, automated approaches in both the identification of new ORFs and the determination of high-resolution 3-D protein structures. A significant bottleneck is reached, however, at the stage of functional annotation in that the assignment of function is not readily automatable. It is often the case that the initial structural analysis at best indicates a functional family for a given hypothetical protein, but further identification of a relevant ligand or substrate is impeded by the diversity of function in a particular structural classification of proteins family, a highly selective and specific ligand-binding site, or the identification of a novel protein fold. Our approach to the functional annotation of hypothetical proteins relies on the combination of structural information with additional bioinformatics evidence garnered from operon prediction, loose functional information of additional operon members, conservation of catalytic residues, as well as cocrystallization trials and virtual ligand screening. The synthesis of all available information for each protein has permitted the functional annotation of several hypothetical proteins from Escherichia coli and each assignment has been confirmed through generally accepted biochemical methods.

Crystal Structure of RS21-C6, Involved in Nucleoside Triphosphate Pyrophosphohydrolysis

RS21-C6, which is highly expressed in all vertebrate genomes and green plants, is proposed to have nucleoside triphosphate pyrophosphohydrolase activity. Here, we report the crystal structures of the core fragment of RS21-C6, named RSCUT, and the complex with the substrate 5-methyl dCTP. The refined structure of RSCUT consists mainly of alpha-helices and shows formation of a tightly associated tetramer. On the basis of the structure of the RSCUT-m5dCTP complex and the results of pyrophosphatase activity assays, several key residues involved in the substrate binding of RS21-C6 have been identified. Tetramer formation is shown to be required for substrate binding.

Crystal structure of VC0702 at 2.0 A: conserved hypothetical protein from Vibrio cholerae

VC0702, a conserved hypothetical protein of unknown function from Vibrio cholerae, resides in a three-gene operon containing the MbaA gene that encodes for a GGDEF and EAL domain-containing protein which is involved in regulating formation of the extracellular matrix of biofilms in Vibrio cholerae. The VC0702 crystal structure has been determined at 2.0 A and refined to Rwork = 22.8% and Rfree = 26.3%. VC0702 crystallized in an orthorhombic crystal lattice in the C222(1) space group with dimensions of a = 66.61 A, b = 88.118 A, and c = 118.35 A with a homodimer in the asymmetric unit. VC0702, which forms a mixed alpha + beta three-layered alphabetaalpha sandwich, belongs to the Pfam DUF84 and COG1986 families of proteins. Sequence conservation within the DUF84 and COG1986 families was used to identify a conserved patch of surface residues that define a cleft and potential substrate-binding site in VC0702. The three-dimensional structure of VC0702 is similar to that of Mj0226 from Methanococcus janeschii, which has been identified as a novel NTPase that binds NTP in a deep cleft similarly located to the conserved patch of surface residues that define an analogous cleft in VC0702. Collectively, the data suggest that VC0702 may have a biochemical function that involves NTP binding and phosphatase activity of some kind, and is likely involved in regulation of the signaling pathway that controls biofilm formation and maintenance in Vibrio cholerae.

House cleaning, a part of good housekeeping

Cellular metabolism constantly generates by-products that are wasteful or even harmful. Such compounds are excreted from the cell or are removed through hydrolysis to normal cellular metabolites by various 'house-cleaning' enzymes. Some of the most important contaminants are non-canonical nucleoside triphosphates (NTPs) whose incorporation into the nascent DNA leads to increased mutagenesis and DNA damage. Enzymes intercepting abnormal NTPs from incorporation by DNA polymerases work in parallel with DNA repair enzymes that remove lesions produced by modified nucleotides. House-cleaning NTP pyrophosphatases targeting non-canonical NTPs belong to at least four structural superfamilies: MutT-related (Nudix) hydrolases, dUTPase, ITPase (Maf/HAM1) and all-alpha NTP pyrophosphatases (MazG). These enzymes have high affinity (Km's in the micromolar range) for their natural substrates (8-oxo-dGTP, dUTP, dITP, 2-oxo-dATP), which allows them to select these substrates from a mixture containing a approximately 1000-fold excess of canonical NTPs. To date, many house-cleaning NTPases have been identified only on the basis of their side activity towards canonical NTPs and NDP derivatives. Integration of growing structural and biochemical data on these superfamilies suggests that their new family members cleanse the nucleotide pool of the products of oxidative damage and inappropriate methylation. House-cleaning enzymes, such as 6-phosphogluconolactonase, are also part of normal intermediary metabolism. Genomic data suggest that house-cleaning systems are more abundant than previously thought and include numerous analogous enzymes with overlapping functions. We discuss the structural diversity of these enzymes, their phylogenetic distribution, substrate specificity and the problem of identifying their true substrates.

YjjX: from structure "Tu" function

It has been shown by structural analysis that YjjX, a hypothetical protein in E. coli, is an ITPase/XTPase and suggest that it may play dual roles in prokaryotic translational regulation and oxidative cell stress response.

Identification of an ITPase/XTPase in Escherichia coli by structural and biochemical analysis

Inosine triphosphate (ITP) and xanthosine triphosphate (XTP) are formed upon deamination of ATP and GTP as a result of exposure to chemical mutagens and oxidative damage. Nucleic acid synthesis requires safeguard mechanisms to minimize undesired lethal incorporation of ITP and XTP. Here, we present the crystal structure of YjjX, a protein of hitherto unknown function. The three-dimensional fold of YjjX is similar to those of Mj0226 from Methanococcus janschii, which possesses nucleotidase activity, and of Maf from Bacillus subtilis, which can bind nucleotides. Biochemical analyses of YjjX revealed it to exhibit specific phosphatase activity for inosine and xanthosine triphosphates and have a possible interaction with elongation factor Tu. The enzymatic activity of YjjX as an inosine/xanthosine triphosphatase provides evidence for a plausible protection mechanism by clearing the noncanonical nucleotides from the cell during oxidative stress in E. coli.

The utility of residual dipolar couplings in detecting motion in carbohydrates: application to sucrose

The solution structure and dynamics of sucrose are examined using a combination of NMR residual dipolar coupling and molecular mechanics force fields. It is found that the alignment tensors of the individual rings are different, and that fitting 35 measured residual dipolar couplings to structures with specific phi, psi values indicates the presence of three major conformations: phi, psi=(120 degrees ,270 degrees), (45 degrees, 300 degrees) and (90 degrees ,180 degrees). Furthermore, fitting two structures simultaneously to the 35 residual dipolar couplings results in a substantial improvement in the fits. The existence of multiple conformations having similar stabilities is a strong indication of motion, due to the interconversion among these states. Results from four molecular mechanics force fields are in general agreement with the experimental results. However, there are major disagreements between force fields. Because fits of residual dipolar couplings to structures are dependent on the force field used to calculate the structures, multiple force fields were used to interpret NMR data. It is demonstrated that the pucker of the fructofuranosyl ring affects the calculated potential energy surface, and the fit to the residual dipolar couplings data. Previously published 13C nuclear relaxation results suggesting that sucrose is rigid are not inconsistent with the present results when motional timescales are considered.

Dimeric dUTPases, HisE, and MazG belong to a new superfamily of all-alpha NTP pyrophosphohydrolases with potential "house-cleaning" functions

Structure-guided analysis of the new dimeric dUTPase family revealed its sequence relationship to the phage T4 dCTPase, phosphoribosyl-ATP pyrophosphatase HisE, NTP pyrophosphatase MazG, and several uncharacterized protein families, including the human protein XTP3TPA (RS21-C6), which is overexpressed in embryonic and cancer cells. Comparison with the recently determined structure of a MazG-like protein from Sulfolobus solfataricus supported the unification of these enzymes in one superfamily of all-alpha NTP pyrophosphatases, suggesting that dimeric dUTPases evolved from a tetrameric MazG-like ancestor by gene duplication. Analysis of the structure of the Sulfolobus MazG points to 2-hydroxyadenosine (isoguanosine) triphosphate, a product of oxidative damage of ATP, as the most likely substrate. We predict that uncharacterized members of this superfamily perform "house-cleaning" functions by hydrolyzing abnormal NTPs and are functionally analogous to the structurally unrelated hydrolases of the Nudix superfamily. We outline probable tertiary and quaternary structures of the all-alpha NTP pyrophosphatase superfamily members.

Analysis of protein structures reveals regions of rare backbone conformation at functional sites

Regions of rare conformation were located in 300 protein crystal structures representing seven major protein folds. A distance matrix algorithm was used to search rapidly for 9-residue fragments of rare backbone conformation using a comparison to a relational database of encoded fragments derived from the database of nonredundant structures. Rare fragments were found in 61% of the analyzed protein structures. Detailed analysis was performed for 78 proteins of different folds. The rare fragments were located near functional sites in 72% of the protein structures. The rare fragments often formed parts of ligand-binding sites (59%), protein-protein interfaces (8%), and domain-domain contacts (5%). Of the remaining structures, 5% had a high average B-factor or high local B-factors. Statistical analysis suggests that the association between ligands and rare regions does not occur by chance alone. The present study is likely to underestimate the number of functional sites, because not all analyzed protein structures contained a ligand. The results suggest that rapid searches for regions with rare local backbone conformations can assist in prediction of functional sites in novel proteins.

Annotation in three dimensions. PINTS: Patterns in Non-homologous Tertiary Structures

The detection of local structural patterns in proteins (e.g. active sites) can provide insights into protein function in the absence of sequence or fold similarity. Methods to detect such similarities are key during structural annotation, for example with results from Structural Genomics initiatives. PINTS (Patterns in Non-homologous Tertiary Structures, http://pints.embl.de) performs database searches for such patterns and most importantly provides a measure of statistical significance for any similarity uncovered. To aid functional annotation of proteins, we allow comparisons of pre-defined patterns against databases of complete structures and of entire structures to databases of particular residues likely to be functionally important.

Selenium-assisted nucleic acid crystallography: use of phosphoroselenoates for MAD phasing of a DNA structure

The combination of synchrotron radiation and a variety of atoms or ions (either covalently attached to the biomolecule prior to crystallization or soaked into crystals) that serve as anomalous scatterers constitutes a powerful tool in the X-ray crystallographer's repertoire of structure determination techniques. Phosphoroselenoates in which one of the nonbridging phosphate oxygens in the backbone is replaced by selenium offer a simplified means for introducing an anomalous scatterer into oligonucleotides by conventional solid-phase synthesis. Unlike other methods that are used to derivatize DNA or RNA by covalent attachment of a heavy atom (i.e., bromine at the C5 position of pyrimidines), tedious synthesis of specialized nucleosides is not required. Introduction of selenium is readily accomplished in solid-phase oligonucleotide synthesis by replacing the standard oxidation agent with a solution of potassium selenocyanide. This results in a diastereomeric mixture of phosphoroselenoates that can be separated by strong anion-exchange HPLC. As a test case, all 10 DNA hexamers of the sequence CGCGCG containing a single phosphoroselenoate linkage (PSe) were prepared. Crystals were grown for a subset of them, and the structure of [d(C(PSe)GCGCG)](2) was determined by the multiwavelength anomalous dispersion technique and refined to 1.1 A resolution.

Covalent incorporation of selenium into oligonucleotides for X-ray crystal structure determination via MAD: proof of principle. Multiwavelength anomalous dispersion

Selenium was incorporated into an oligodeoxynucleotide in the form of 2'-methylseleno-uridine (U(Se)). The X-ray crystal structure of the duplex left open bracket d(GCGTA)U(Se)d(ACGC) right open bracket (2) was determined by the multiwavelength anomalous dispersion (MAD) technique and refined to a resolution of 1.3 A, demonstrating that selenium can selectively substitute oxygen in DNA and that the resulting compounds are chemically stable. Since derivatization at the 2'-alpha-position with selenium does not affect the preference of the sugar for the C3'-endo conformation, this strategy is suitable for incorporating selenium into RNA. The availability of selenium-containing nucleic acids for crystallographic phasing offers an attractive alternative to the commonly used halogenated pyrimidines.

Microarray analysis of the Bacillus subtilis K-state: genome-wide expression changes dependent on ComK

In Bacillus subtilis, the competence transcription factor ComK activates its own transcription as well as the transcription of genes that encode DNA transport proteins. ComK is expressed in about 10% of the cells in a culture grown to competence. Using DNA microarrays representing approximately 95% of the protein-coding open reading frames in B. subtilis, we compared the expression profiles of wild-type and comK strains, as well as of a mecA mutant (which produces active ComK in all the cells of the population) and a comK mecA double mutant. In these comparisons, we identified at least 165 genes that are upregulated by ComK and relatively few that are downregulated. The use of reporter fusions has confirmed these results for several genes. Many of the ComK-regulated genes are organized in clusters or operons, and 23 of these clusters are preceded by apparent ComK-box promoter motifs. In addition to those required for DNA uptake, other genes that are upregulated in the presence of ComK are probably involved in DNA repair and in the uptake and utilization of nutritional sources. From this and previous work, we conclude that the ComK regulon defines a growth-arrested state, distinct from sporulation, of which competence for genetic transformation is but one notable feature. We suggest that this is a unique adaptation to stress and that it be termed the 'K-state'.

A tour of structural genomics

Structural genomics projects aim to provide an experimental or computational three-dimensional model structure for all of the tractable macromolecules that are encoded by complete genomes. To this end, pilot centres worldwide are now exploring the feasibility of large-scale structure determination. Their experimental structures and computational models are expected to yield insight into the molecular function and mechanism of thousands of proteins. The pervasiveness of this information is likely to change the use of structure in molecular biology and biochemistry.

Solution structure and function of a conserved protein SP14.3 encoded by an essential Streptococcus pneumoniae gene

Streptococcus pneumoniae is a major human pathogen that causes high mortality and morbidity rates and has developed resistance to many antibiotics. The genome of S. pneumoniae has recently been completely sequenced revealing many genes encoding hypothetical proteins of unknown function. We have found that the gene encoding one such conserved protein, SP14.3, is essential for growth of S. pneumonia. Since it is essential, SP14.3 represents a potential target for drug discovery. Here, we describe the three-dimensional solution structure of SP14.3 as determined by NMR spectroscopy. The structure consists of two domains each with an alpha/beta-fold. The N-terminal domain contains two alpha-helices and a three-stranded beta-sheet, while the C-terminal domain is composed of one alpha-helix and a five-stranded beta-sheet. The N-terminal domain of the protein contains a highly negatively charged surface and resembles the fold of the N-terminal domain of Thermus thermophilus ribosomal protein S3. The C-terminal domain has a protein fold similar to human small nuclear ribonucleoprotein Sm D3 and Haloarcula marismortui ribosomal protein L21E. The two domains of the protein tumble in solution overall as a whole with an overall molecular rotational correlation time (tau(m)) of 12.9 ns at 25 degrees C. The relative orientation of the two domains is not defined by the nuclear Overhauser effect data. Indeed, residual dipolar couplings and the structure calculations indicate that the relative orientation of the two domains is not rigidly oriented with respect to one another in solution.

Structural genomics: opportunities and challenges

Following the complete genome sequencing of an increasing number of organisms, structural biology is engaging in a systematic approach of high-throughput structure determination called structural genomics to create a complete inventory of protein folds/structures that will help predict functions for all proteins. First results show that structural genomics will be highly effective in finding functional annotations for proteins of unknown function.

Determination of protein function, evolution and interactions by structural genomics

The genome sequencing projects and knowledge of the entire protein repertoires of many organisms have prompted new procedures and techniques for the large-scale determination of protein structure, function and interactions. Recently, new work has been carried out on the determination of the function and evolutionary relationships of proteins by experimental structural genomics, and the discovery of protein-protein interactions by computational structural genomics.

The structure of the yrdC gene product from Escherichia coli reveals a new fold and suggests a role in RNA binding

The yrdC family of genes codes for proteins that occur both independently and as a domain in proteins that have been implicated in regulation. An example for the latter case is the sua5 gene from yeast. SuaS was identified as a suppressor of a translation initiation defect in cytochrome c and is required for normal growth in yeast (Na JG, Pinto I, Hampsey M, 1992, Genetics 11:791-801). However, the function of the Sua5 protein remains unknown; Sua5 could act either at the transcriptional or the posttranscriptional levels to compensate for an aberrant translation start codon in the cyc gene. To potentially learn more about the function of YrdC and proteins featuring this domain, the crystal structure of the YrdC protein from Escherichia coli was determined at a resolution of 2.0 A. YrdC adopts a new fold with no obvious similarity to those of other proteins with known three-dimensional (3D) structure. The protein features a large concave surface on one side that exhibits a positive electrostatic potential. The dimensions of this depression, its curvature, and the fact that conserved basic amino acids are located at its floor suggest that YrdC may be a nucleic acid binding protein. An investigation of YrdC's binding affinities for single- and double-stranded RNA and DNA fragments as well as tRNAs demonstrates that YrdC binds preferentially to double-stranded RNA. Our work provides evidence that 3D structures of functionally uncharacterized gene products with unique sequences can yield novel folds and functional insights.