With many genomes sequenced, a pressing challenge in biology is predicting the function of the proteins that the genes encode. When proteins are unrelated to others of known activity, bioinformatics inference for function becomes problematic. It would thus be useful to interrogate protein structures for function directly. Here, we predict the function of an enzyme of unknown activity, Tm0936 from Thermotoga maritima, by docking high-energy intermediate forms of thousands of candidate metabolites. The docking hit list was dominated by adenine analogues, which appeared to undergo C6-deamination. Four of these, including 5-methylthioadenosine and Sadenosylhomocysteine (SAH), were tested as substrates, and three had substantial catalytic rate constants (10 5 M −1 s −1 ). The X-ray crystal structure of the complex between Tm0936 and the product resulting from the deamination of SAH, S-inosylhomocysteine, was determined, and it corresponded closely to the predicted structure. The deaminated products can be further metabolized by T. maritima in a previously uncharacterized SAH degradation pathway. Structure-based docking with high-energy forms of potential substrates may be a useful tool to annotate enzymes for function.For enzymes of unknown function, substrate prediction based on structural complementarity is, in principle, an alternative to bio-informatics inference of function 1,2 . Structure-based prediction becomes attractive when the target enzyme has little relationship to orthologues of known activity, making inference unreliable 3,4 . Whereas structure-based prediction has been used with some successes for inhibitor design, substrate prediction has proven difficult 5-8 . In addition to the well-known problems of sampling and scoring in docking, substrate prediction confronts several additional challenges. These include the many possible substrates to consider and the many reactions that an enzyme might catalyse 9-11 . Furthermore, enzymes preferably recognize transition states over the ground state structures that are usually represented in docking databases 12-14 .Correspondence and requests for materials related to docking should be addressed to B.K.S. (shoichet@cgl.ucsf.edu). J.C.H designed the docking database, performed the docking runs, and analysed the docking results. F.M.R. and R.M.-A. performed the enzymatic characterization of Tm0936 and Tm0172, including cloning and purification of the proteins. S.C.A., E.F. and A.A.F. determined the X-ray structure of Tm0936 with S-inosyl-homocysteine. J.C.H. and B.K.S. largely wrote the paper. All authors discussed the results and commented on the manuscript. Author InformationThe complex structure of Tm0936 with SIH has been deposited in the PDB (accession code 2PLM). Reprints and permissions information is available at www.nature.com/reprints. The authors declare no competing financial interests. These retrospective results encouraged us to prospectively predict the substrates of Tm0936 from T. maritima. The X-ray structure of the enzyme had been deter...

Isoaspartyl dipeptidase (IAD) is a member of the amidohydrolase superfamily and catalyzes the hydrolytic cleavage of beta-aspartyl dipeptides. Structural studies of the wild-type enzyme have demonstrated that the active site consists of a binuclear metal center positioned at the C-terminal end of a (beta/alpha)(8)-barrel domain. Steady-state kinetic parameters for the hydrolysis of beta-aspartyl dipeptides were obtained at pH 8.1. The pH-rate profiles for the hydrolysis of beta-Asp-Leu were obtained for the Zn/Zn-, Co/Co-, Ni/Ni-, and Cd/Cd-substituted forms of IAD. Bell-shaped profiles were observed for k(cat) and k(cat)/K(m) as a function of pH for all four metal-substituted forms. The pK(a) of the group that must be unprotonated for catalytic activity varied according to the specific metal ion bound in the active site, whereas the pK(a) of the group that must be protonated for catalytic activity was relatively independent of the specific metal ion present. The identity of the group that must be unprotonated for catalytic activity was consistent with the hydroxide that bridges the two divalent cations of the binuclear metal center. The identity of the group that must be protonated for activity was consistent with the free alpha-amino group of the dipeptide substrate. Kinetic constants were obtained for the mutant enzymes at conserved residues Glu77, Tyr137, Arg169, Arg233, Asp285, and Ser289. The catalytic properties of the wild-type and mutant enzymes, coupled with the X-ray crystal structure of the D285N mutant complexed with beta-Asp-His, are consistent with a chemical reaction mechanism for the hydrolysis of dipeptides that is initiated by the polarization of the amide bond via complexation to the beta-metal ion of the binuclear metal center. Nucleophilic attack by the bridging hydroxide is facilitated by abstraction of its proton by the side chain carboxylate of Asp285. Collapse of the tetrahedral intermediate and cleavage of the carbon-nitrogen bond occur with donation of a proton from the protonated form of Asp285.

An enzyme from Pseudomonas aeruginosa, Pa0142 (gi|9945972) has been identified for the first time that is able to catalyze the deamination of 8-oxoguanine (8-oxoG) to uric acid. 8-Oxoguanine is formed by the oxidation of guanine residues within DNA by reactive oxygen species and this lesion results in the G:C to T:A transversions. The value of kcat/Km for the deamination of 8-oxoG by Pa0142 at pH 8.0 and 30 °C is 2.0 × 104 M−1 s−1. This enzyme can also catalyze the deamination of isocystosine and guanine at rates that are approximately an order of magnitude slower. The three-dimensional structure of a homologous enzyme (gi|44264246) from the Sargasso Sea has been determined by x-ray diffraction methods to a resolution of 2.2Å (PDB code: 3h4u). The enzyme folds as a (β/α)8− barrel and it is a member of the amidohydrolase superfamily with a single zinc in the active site. This enzyme catalyzes the deamination of 8-oxoG with a value of kcat/Km of 2.7 × 105 M−1 s−1. Computational docking of potential high energy intermediates for the deamination reaction to the x-ray crystal structure suggests that the active site binding of 8-oxoG is facilitated by hydrogen bond interactions from a conserved glutamine that follows β-strand 1 with O6, a conserved tyrosine that follows β-strand 2 with N7, and a conserved cysteine residue that follows β-strand 4 with O8. A bioinformatic analysis of available protein sequences suggest that approximately 200 other bacteria possess an enzyme capable of catalyzing the deamination of 8-oxoG.

The functional assignment of enzymes that catalyze unknown chemical transformations is a difficult problem. The protein Pa5106 from Pseudomonas aeruginosa has been identified as a member of the amidohydrolase superfamily by a comprehensive amino acid sequence comparison with structurally authenticated members of this superfamily. The function of Pa5106 has been annotated as a probablechlorohydrolase or cytosine deaminase. A close examination of the genomic content of P. aeruginosa reveals that the gene for this protein is in close proximity to genes included in the histidine degradation pathway. The first three steps for the degradation of histidine include the action of HutH, HutU, and HutI to convert L-histidine to N-formimino-L-glutamate. The degradation of N-formimino-L-glutamate to L-glutamate can occur by three different pathways. Three proteins in P. aeruginosa have been identified that catalyze two of the three possible pathways for the degradation of N-formimino-L-glutamate. The protein Pa5106 was shown to catalyze the deimination of N-formimino-L-glutamate to ammonia and N-formyl-L-glutamate, while Pa5091 catalyzed the hydrolysis of N-formyl-L-glutamate to formate and L-glutamate. The protein Pa3175 is dislocated from the hut operon and was shown to catalyze the hydrolysis of N-formimino-L-glutamate to formamide and L-glutamate. The reason for the coexistence of two alternative pathways for the degradation of N-formimino-L-glutamate in P. aeruginosa is unknown.

The functional role of the alphaM3 transmembrane domain of the Torpedo nicotinic acetylcholine receptor (AChR) was characterized by performing tryptophan-scanning mutagenesis at 13 positions within alphaM3, from residue M278 through I290. The expression of the mutants in Xenopus oocytes was measured by [(125)I]-alpha-bungarotoxin binding, and ACh receptor function was evaluated by using a two-electrode voltage clamp. Six mutants (L279W, F280W, I283W, V285W, S288W, and I289W) were expressed at lower levels than the wild type. Most of these residues have been proposed to face the interior of the protein. The I286W mutant was expressed at 2.4-fold higher levels than the wild type, and the two lipid-exposed mutations, F284W and S287W, were expressed at similar levels as wild type. Binding assays indicated that the alphaM3 domain can accommodate bulky groups in almost all positions. Three mutations, M282W, V285W, and I289W, caused a loss of receptor function, suggesting that the tryptophan side chains alter the conformational changes required for channel assembly or ion channel function. This loss of function suggests that these positions may be involved in helix-helix contacts that are critical for channel gating. The lipid-exposed mutation F284W enhances the receptor macroscopic response at low ACh concentrations and decreases the EC(50). Taken together, our results suggest that alphaM3 contributes to the gating machinery of the nicotinic ACh receptor and that alphaM3 is comprised of a mixture of two types of helical structures.

Isoaspartyl dipeptidase from Escherichia coli functions in protein degradation by catalyzing the hydrolysis of beta-L-isoaspartyl linkages in dipeptides. The best substrate for the enzyme reported thus far is iso-Asp-Leu. Here we report the X-ray analysis of the enzyme in its resting state and complexed with aspartate to 1.65 and 2.1 A resolution, respectively. The quaternary structure of the enzyme is octameric and can be aptly described as a tetramer of dimers. Each subunit folds into two distinct domains: the N-terminal region containing eight strands of mixed beta-sheet and the C-terminal motif that is dominated by a (beta,alpha)(8)-barrel. A binuclear zinc center is located in each subunit at the C-terminal end of the (beta,alpha)(8)-barrel. Ligands to the binuclear metal center include His 68, His 70, His 201, His 230, and Asp 285. The two zincs are bridged by a carboxylated lysine residue (Lys 162) and a solvent molecule, most likely a hydroxide ion. The product of the reaction, aspartate, binds to the enzyme by displacing the bridging solvent with its side chain functional group. From this investigation it is proposed that the reaction mechanism of the enzyme proceeds through a tetrahedral intermediate and that the bridging solvent attacks the re face of the carbonyl carbon of the scissile peptide bond. This structural analysis confirms the placement of isoaspartyl dipeptidase into the urease-related amidohydrolase superfamily.

Mycolactone, the amphiphilic macrolide toxin secreted by Mycobacterium ulcerans, plays a significant role in the pathology and manifestations of Buruli ulcer (BU). Consequently, it follows that the toxin is a suitable target for the development of diagnostics and therapeutics for this disease. Yet, several challenges have deterred such development. For one, the lipophilic nature of the toxin makes it difficult to handle and store and contributes to variability associated with laboratory experimentation and purification yields. In this manuscript, we have attempted to incorporate our understanding of the lipophilicity of mycolactone in order to define the optimal methods for the storage, handling, and purification of this toxin. We present a systematic correlation of variability associated with measurement techniques (thin-layer chromatography (TLC), mass spectrometry (MS), and UV-Vis spectrometry), storage conditions, choice of solvents, as well as the impact of each of these on toxin function as assessed by cellular cytotoxicity. We also compared natural mycolactone extracted from bacterial culture with synthesized toxins in laboratory (solvents, buffers) and physiologically relevant (serum) matrices. Our results point to the greater stability of mycolactone in organic, as well as detergent-containing, solvents, regardless of the container material (plastic, glass, or silanized tubes). They also highlight the presence of toxin in samples that may be undetectable by any one technique, suggesting that each detection approach captures different configurations of the molecule with varying specificity and sensitivity. Most importantly, our results demonstrate for the very first time that amphiphilic mycolactone associates with host lipoproteins in serum, and that this association will likely impact our ability to study, diagnose, and treat Buruli ulcers in patients.

Until recently, engineering strategies for altering gene expression have focused on transcription control using strong inducible promoters or one of several methods to knock down wasteful genes. Recently, synthetic riboregulators have been developed for translational regulation of gene expression. Here, we report a new modular synthetic riboregulator class that has the potential to finely tune protein expression and independently control the concentration of each enzyme in an engineered metabolic pathway. This development is important because the most straightforward approach to altering the flux through a particular metabolic step is to increase or decrease the concentration of the enzyme. Our design includes a cis-repressor at the 5' end of the mRNA that forms a stem-loop helix, occluding the ribosomal binding sequence and blocking translation. A trans-expressed activating-RNA frees the ribosomal-binding sequence, which turns on translation. The overall architecture of the riboregulators is designed using Watson-Crick base-pairing stability. We describe here a cis-repressor that can completely shut off translation of antibiotic-resistance reporters and a trans-activator that restores translation. We have established that it is possible to use these riboregulators to achieve translational control of gene expression over a wide dynamic range. We have also found that a targeting sequence can be modified to develop riboregulators that can, in principle, independently regulate translation of many genes. In a selection experiment, we demonstrated that by subtly altering the sequence of the trans-activator it is possible to alter the ratio of the repressed and activated states and to achieve intermediate translational control.