Faculty Bibliography

Proteins of unknown function belonging to cog1816 and cog0402 were characterized. Sav2595 from Steptomyces avermitilis MA-4680, Acel0264 from Acidothermus cellulolyticus 11B, Nis0429 from Nitratiruptor sp. SB155-2 and Dr0824 from Deinococcus radiodurans R1 were cloned, purified, and their substrate profiles determined. These enzymes were previously incorrectly annotated as adenosine deaminases or chlorohydrolases. It was shown here that these enzymes actually deaminate 6-aminodeoxyfutalosine. The deamination of 6-aminodeoxyfutalosine is part of an alternative menaquinone biosynthetic pathway that involves the formation of futalosine. 6-Aminodeoxyfutalosine is deaminated by these enzymes with catalytic efficiencies greater than 10(5) M(-1) s(-1), Km values of 0.9-6.0 μM, and kcat values of 1.2-8.6 s(-1). Adenosine, 2'-deoxyadenosine, thiomethyladenosine, and S-adenosylhomocysteine are deaminated at least an order of magnitude slower than 6-aminodeoxyfutalosine. The crystal structure of Nis0429 was determined and the substrate, 6-aminodeoxyfutalosine, was positioned in the active site on the basis of the presence of adventitiously bound benzoic acid. In this model, Ser-145 interacts with the carboxylate moiety of the substrate. The structure of Dr0824 was also determined, but a collapsed active site pocket prevented docking of substrates. A computational model of Sav2595 was built on the basis of the crystal structure of adenosine deaminase and substrates were docked. The model predicted a conserved arginine after β-strand 1 to be partially responsible for the substrate specificity of Sav2595.

The haloacid dehalogenase enzyme superfamily (HADSF) is largely composed of phosphatases that have been particularly successful at adaptating to novel biological functions relative to members of other phosphatase families. Herein, we examine the structural basis for the divergence of function in two bacterial homologues: 2-keto-3-deoxy-d-manno-octulosonate 8-phosphate phosphohydrolase (KDO8P phosphatase, KDO8PP) and 2-keto-3-deoxy-9-O-phosphonononic acid phosphohydrolase (KDN9P phosphatase, KDN9PP). KDO8PP and KDN9PP catalyze the final step in KDO and KDN synthesis, respectively, prior to transfer to CMP to form the activated sugar nucleotide. KDO8PP and KDN9PP orthologs derived from an evolutionarily diverse collection of bacterial species were subjected to steady-state kinetic analysis to determine their specificities toward catalyzed KDO8P and KDN9P hydrolysis. Although each enzyme was more active with its biological substrate, the degree of selectivity (as defined by the ratio of kcat/Km for KDO8P vs KDN9P) varied significantly. High-resolution X-ray structure determination of Haemophilus influenzae KDO8PP bound to KDO/VO3(-) and Bacteriodes thetaiotaomicron KDN9PP bound to KDN/VO3(-) revealed the substrate-binding residues. The structures of the KDO8PP and KDN9PP orthologs were also determined to reveal the differences in their active-site structures that underlie the variation in substrate preference. Bioinformatic analysis was carried out to define the sequence divergence among KDN9PP and KDO8PP orthologs. The KDN9PP orthologs were found to exist as single-domain proteins or fused with the pathway nucleotidyl transferases; the fusion of KDO8PP with the transferase is rare. The KDO8PP and KDN9PP orthologs share a stringently conserved Arg residue that forms a salt bridge with the substrate carboxylate group. The split of the KDN9PP lineage from the KDO8PP orthologs is easily tracked by the acquisition of a Glu/Lys pair that supports KDN9P binding. Moreover, independently evolved lineages of KDO8PP orthologs exist, and are separated by diffuse active-site sequence boundaries. We infer a high tolerance of the KDO8PP catalytic platform to amino acid replacements that in turn influence substrate specificity changes and thereby facilitate the divergence in biological function.

The identification of novel metabolites and the characterization of their biological functions are major challenges in biology. X-ray crystallography can reveal unanticipated ligands that persist through purification and crystallization. These adventitious protein-ligand complexes provide insights into new activities, pathways and regulatory mechanisms. We describe a new metabolite, carboxy-S-adenosyl-l-methionine (Cx-SAM), its biosynthetic pathway and its role in transfer RNA modification. The structure of CmoA, a member of the SAM-dependent methyltransferase superfamily, revealed a ligand consistent with Cx-SAM in the catalytic site. Mechanistic analyses showed an unprecedented role for prephenate as the carboxyl donor and the involvement of a unique ylide intermediate as the carboxyl acceptor in the CmoA-mediated conversion of SAM to Cx-SAM. A second member of the SAM-dependent methyltransferase superfamily, CmoB, recognizes Cx-SAM and acts as a carboxymethyltransferase to convert 5-hydroxyuridine into 5-oxyacetyl uridine at the wobble position of multiple tRNAs in Gram-negative bacteria, resulting in expanded codon-recognition properties. CmoA and CmoB represent the first documented synthase and transferase for Cx-SAM. These findings reveal new functional diversity in the SAM-dependent methyltransferase superfamily and expand the metabolic and biological contributions of SAM-based biochemistry. These discoveries highlight the value of structural genomics approaches in identifying ligands within the context of their physiologically relevant macromolecular binding partners, and in revealing their functions.

The number of available protein sequences has increased exponentially with the advent of high-throughput genomic sequencing, creating a significant challenge for functional annotation. Here, we describe a large-scale study on assigning function to unknown members of the trans-polyprenyl transferase (E-PTS) subgroup in the isoprenoid synthase superfamily, which provides substrates for the biosynthesis of the more than 55,000 isoprenoid metabolites. Although the mechanism for determining the product chain length for these enzymes is known, there is no simple relationship between function and primary sequence, so that assigning function is challenging. We addressed this challenge through large-scale bioinformatics analysis of >5,000 putative polyprenyl transferases; experimental characterization of the chain-length specificity of 79 diverse members of this group; determination of 27 structures of 19 of these enzymes, including seven cocrystallized with substrate analogs or products; and the development and successful application of a computational approach to predict function that leverages available structural data through homology modeling and docking of possible products into the active site. The crystallographic structures and computational structural models of the enzyme-ligand complexes elucidate the structural basis of specificity. As a result of this study, the percentage of E-PTS sequences similar to functionally annotated ones (BLAST e-value ≤ 1e(-70)) increased from 40.6 to 68.8%, and the percentage of sequences similar to available crystal structures increased from 28.9 to 47.4%. The high accuracy of our blind prediction of newly characterized enzymes indicates the potential to predict function to the complete polyprenyl transferase subgroup of the isoprenoid synthase superfamily computationally.

The rapid advance in genome sequencing presents substantial challenges for protein functional assignment, with half or more of new protein sequences inferred from these genomes having uncertain assignments. The assignment of enzyme function in functionally diverse superfamilies represents a particular challenge, which we address through a combination of computational predictions, enzymology, and structural biology. Here we describe the results of a focused investigation of a group of enzymes in the enolase superfamily that are involved in epimerizing dipeptides. The first members of this group to be functionally characterized were Ala-Glu epimerases in Eschericiha coli and Bacillus subtilis, based on the operon context and enzymological studies; these enzymes are presumed to be involved in peptidoglycan recycling. We have subsequently studied more than 65 related enzymes by computational methods, including homology modeling and metabolite docking, which suggested that many would have divergent specificities;, i.e., they are likely to have different (unknown) biological roles. In addition to the Ala-Phe epimerase specificity reported previously, we describe the prediction and experimental verification of: (i) a new group of presumed Ala-Glu epimerases; (ii) several enzymes with specificity for hydrophobic dipeptides, including one from Cytophaga hutchinsonii that epimerizes D-Ala-D-Ala; and (iii) a small group of enzymes that epimerize cationic dipeptides. Crystal structures for certain of these enzymes further elucidate the structural basis of the specificities. The results highlight the potential of computational methods to guide experimental characterization of enzymes in an automated, large-scale fashion.

The explosion of protein sequence information requires that current strategies for function assignment evolve to complement experimental approaches with computationally based function prediction. This necessitates the development of strategies based on the identification of sequence markers in the form of specificity determinants and a more informed definition of orthologues. Herein, we have undertaken the function assignment of the unknown haloalkanoate dehalogenase superfamily member BT2127 (Uniprot accession code Q8A5 V9) from Bacteroides thetaiotaomicron using an integrated bioinformatics-structure-mechanism approach. The substrate specificity profile and steady-state rate constants of BT2127 (with a k(cat)/K(m) value for pyrophosphate of ~1 × 10(5) M(-1) s(-1)), together with the gene context, support the assigned in vivo function as an inorganic pyrophosphatase. The X-ray structural analysis of wild-type BT2127 and several variants generated by site-directed mutagenesis shows that substrate discrimination is based, in part, on active site space restrictions imposed by the cap domain (specifically by residues Tyr76 and Glu47). Structure-guided site-directed mutagenesis coupled with kinetic analysis of the mutant enzymes identified the residues required for catalysis, substrate binding, and domain-domain association. On the basis of this structure-function analysis, the catalytic residues Asp11, Asp13, Thr113, and Lys147 as well the metal binding residues Asp171, Asn172, and Glu47 were used as markers to confirm BT2127 orthologues identified via sequence searches. This bioinformatic analysis demonstrated that the biological range of BT2127 orthologue is restricted to the phylum Bacteroidetes/Chlorobi. The key structural determinants in the divergence of BT2127 and its closest homologue, β-phosphoglucomutase, control the leaving group size (phosphate vs glucose phosphate) and the position of the Asp acid/base in the open versus closed conformations. HADSF pyrophosphatases represent a third mechanistic and fold type for bacterial pyrophosphatases.

Decoy Receptor 3 (DcR3), a secreted member of the Tumor Necrosis Factor (TNF) receptor superfamily, neutralizes three different TNF ligands: FasL, LIGHT, and TL1A. Each of these ligands engages unique signaling receptors which direct distinct and critical immune responses. We report the crystal structures of the unliganded DcR3 ectodomain and its complex with TL1A, as well as complementary mutagenesis and biochemical studies. These analyses demonstrate that DcR3 interacts with invariant backbone and side-chain atoms in the membrane-proximal half of TL1A which supports recognition of its three distinct TNF ligands. Additional features serve as antideterminants that preclude interaction with other members of the TNF superfamily. This mode of interaction is unique among characterized TNF:TNFR family members and provides a mechanistic basis for the broadened specificity required to support the decoy function of DcR3, as well as for the rational manipulation of specificity and affinity of DcR3 and its ligands.

Plants and microorganisms reduce environmental inorganic nitrogen to ammonium, which then enters various metabolic pathways solely via conversion of 2-oxoglutarate (2OG) to glutamate and glutamine. Cellular 2OG concentrations increase during nitrogen starvation. We recently identified a family of 2OG-sensing proteins--the nitrogen regulatory protein NrpR--that bind DNA and repress transcription of nitrogen assimilation genes. We used X-ray crystallography to determine the structure of NrpR regulatory domain. We identified the NrpR 2OG-binding cleft and show that residues predicted to interact directly with 2OG are conserved among diverse classes of 2OG-binding proteins. We show that high levels of 2OG inhibit NrpRs ability to bind DNA. Electron microscopy analyses document that NrpR adopts different quaternary structures in its inhibited 2OG-bound state compared with its active apo state. Our results indicate that upon 2OG release, NrpR repositions its DNA-binding domains correctly for optimal interaction with DNA thereby enabling gene repression.

Two uncharacterized enzymes from the amidohydrolase superfamily belonging to cog1228 were cloned, expressed, and purified to homogeneity. The two proteins, Sgx9260c ( gi|44242006 ) and Sgx9260b ( gi|44479596 ), were derived from environmental DNA samples originating from the Sargasso Sea. The catalytic function and substrate profiles for Sgx9260c and Sgx9260b were determined using a comprehensive library of dipeptides and N-acyl derivative of l-amino acids. Sgx9260c catalyzes the hydrolysis of Gly-l-Pro, l-Ala-l-Pro, and N-acyl derivatives of l-Pro. The best substrate identified to date is N-acetyl-l-Pro with a value of k(cat)/K(m) of 3 x 10(5) M(-1) s(-1). Sgx9260b catalyzes the hydrolysis of l-hydrophobic l-Pro dipeptides and N-acyl derivatives of l-Pro. The best substrate identified to date is N-propionyl-l-Pro with a value of k(cat)/K(m) of 1 x 10(5) M(-1) s(-1). Three-dimensional structures of both proteins were determined by X-ray diffraction methods (PDB codes 3MKV and 3FEQ ). These proteins fold as distorted (beta/alpha)(8)-barrels with two divalent cations in the active site. The structure of Sgx9260c was also determined as a complex with the N-methylphosphonate derivative of l-Pro (PDB code 3N2C ). In this structure the phosphonate moiety bridges the binuclear metal center, and one oxygen atom interacts with His-140. The alpha-carboxylate of the inhibitor interacts with Tyr-231. The proline side chain occupies a small substrate binding cavity formed by residues contributed from the loop that follows beta-strand 7 within the (beta/alpha)(8)-barrel. A total of 38 other proteins from cog1228 are predicted to have the same substrate profile based on conservation of the substrate binding residues. The structure of an evolutionarily related protein, Cc2672 from Caulobacter crecentus, was determined as a complex with the N-methylphosphonate derivative of l-arginine (PDB code 3MTW ).

ZBP1 (zipcode-binding protein 1) was originally discovered as a trans-acting factor for the "zipcode" in the 3' untranslated region (UTR) of the beta-actin mRNA that is important for its localization and translational regulation. Subsequently, ZBP1 has been found to be a multifunctional regulator of RNA metabolism that controls aspects of localization, stability, and translation for many mRNAs. To reveal how ZBP1 recognizes its RNA targets, we biochemically characterized the interaction between ZBP1 and the beta-actin zipcode. The third and fourth KH (hnRNP K homology) domains of ZBP1 specifically recognize a bipartite RNA element located within the first 28 nucleotides of the zipcode. The spacing between the RNA sequences is consistent with the structure of IMP1 KH34, the human ortholog of ZBP1, that we solved by X-ray crystallography. The tandem KH domains are arranged in an intramolecular anti-parallel pseudodimer conformation with the canonical RNA-binding surfaces at opposite ends of the molecule. This orientation of the KH domains requires that the RNA backbone must undergo an approximately 180 degrees change in direction in order for both KH domains to contact the RNA simultaneously. The RNA looping induced by ZBP1 binding provides a mechanism for specific recognition and may facilitate the assembly of post-transcriptional regulatory complexes by remodeling the bound transcript.