Faculty Bibliography

We report the discovery of a new potent allosteric effector of sickle cell hemoglobin, GBT440 (36), that increases the affinity of hemoglobin for oxygen and consequently inhibits its polymerization when subjected to hypoxic conditions. Unlike earlier allosteric activators that bind covalently to hemoglobin in a 2:1 stoichiometry, 36 binds with a 1:1 stoichiometry. Compound 36 is orally bioavailable and partitions highly and favorably into the red blood cell with a RBC/plasma ratio of ∼150. This partitioning onto the target protein is anticipated to allow therapeutic concentrations to be achieved in the red blood cell at low plasma concentrations. GBT440 (36) is in Phase 3 clinical trials for the treatment of sickle cell disease (NCT03036813).

A major driver of the pathophysiology of sickle cell disease (SCD) is polymerization of deoxygenated haemoglobin S (HbS), which leads to sickling and destruction of red blood cells (RBCs) and end-organ damage. Pharmacologically increasing the proportion of oxygenated HbS in RBCs may inhibit polymerization, prevent sickling and provide long term disease modification. We report that GBT440, a small molecule which binds to the N-terminal α chain of Hb, increases HbS affinity for oxygen, delays in vitro HbS polymerization and prevents sickling of RBCs. Moreover, in a murine model of SCD, GBT440 extends the half-life of RBCs, reduces reticulocyte counts and prevents ex vivo RBC sickling. Importantly, oral dosing of GBT440 in animals demonstrates suitability for once daily dosing in humans and a highly selective partitioning into RBCs, which is a key therapeutic safety attribute. Thus, GBT440 has the potential for clinical use as a disease-modifying agent in sickle cell patients.

5-Carboxyvanillate decarboxylase (LigW) catalyzes the conversion of 5-carboxyvanillate to vanillate in the biochemical pathway for the degradation of lignin. This enzyme was shown to require Mn(2+) for catalytic activity and the kinetic constants for the decarboxylation of 5-carboxyvanillate by the enzymes from Sphingomonas paucimobilis SYK-6 (kcat = 2.2 s(-1) and kcat/Km = 4.0 × 10(4) M(-1) s(-1)) and Novosphingobium aromaticivorans (kcat = 27 s(-1) and kcat/Km = 1.1 × 10(5) M(-1) s(-1)) were determined. The three-dimensional structures of both enzymes were determined in the presence and absence of ligands bound in the active site. The structure of LigW from N. aromaticivorans, bound with the substrate analogue, 5-nitrovanillate (Kd = 5.0 nM), was determined to a resolution of 1.07 Å. The structure of this complex shows a remarkable enzyme-induced distortion of the nitro-substituent out of the plane of the phenyl ring by approximately 23°. A chemical reaction mechanism for the decarboxylation of 5-carboxyvanillate by LigW was proposed on the basis of the high resolution X-ray structures determined in the presence ligands bound in the active site, mutation of active site residues, and the magnitude of the product isotope effect determined in a mixture of H2O and D2O. In the proposed reaction mechanism the enzyme facilitates the transfer of a proton to C5 of the substrate prior to the decarboxylation step.

We describe a general integrated bioinformatic and experimental strategy to discover the in vitro enzymatic activities and in vivo functions (metabolic pathways) of uncharacterized enzymes discovered in microbial genome projects using the ligand specificities of the solute binding proteins (SBPs) for ABC transporters. Using differential scanning fluorimetry, we determined that the SBP for an ABC transporter encoded by the genome of Mycobacterium smegmatis is stabilized by d-threitol. Using sequence similarity networks and genome neighborhood networks to guide selection of target proteins for pathway enzymes, we applied both in vitro and in vivo experimental approaches to discover novel pathways for catabolism of d-threitol, l-threitol, and erythritol.

Many immune ligands and receptors are potential drug targets, which delicately manipulate a wide range of immune responses. We describe here the successful application of an efficient method to dramatically improve the heterologous expression levels in Drosophila Schneider 2 cells, which enables the high-throughput production of several important immune ligands/receptors for raising antibodies, and for the structural and functional analyses. As an example, we purified the protein and characterized the structure of the immune receptor herpesvirus entry mediator (HVEM, TNFRSF14). HVEM is a member of tumor necrosis factor receptor superfamily, which is recognized by herpes simplex virus glycoprotein D (gD) and facilitates viral entry. HVEM participates in a range of interactions with other cell surface molecules, including LIGHT, BTLA, and CD160 to modulate a wide range of immune processes in CD4(+) and CD8(+) T cells, as well as NK cells. Due to the involvement of HVEM in these diverse signaling interactions, crystal structures of HVEM in complex with gD or BTLA have been previously reported. Here, we report the structure of HVEM in the absence of any ligands.

Pmi1525, an enzyme of unknown function from Proteus mirabilis HI4320 and the amidohydrolase superfamily, was cloned, purified to homogeneity, and functionally characterized. The three-dimensional structure of Pmi1525 was determined with zinc and cacodylate bound in the active site (PDB id: 3RHG ). The structure was also determined with manganese and butyrate in the active site (PDB id: 4QSF ). Pmi1525 folds as a distorted (β/α)8-barrel that is typical for members of the amidohydrolase superfamily and cog1735. The substrate profile for Pmi1525 was determined via a strategy that marshaled the utilization of bioinformatics, structural characterization, and focused library screening. The protein was found to efficiently catalyze the hydrolysis of organophosphonate and carboxylate esters. The best substrates identified for Pmi1525 are ethyl 4-nitrophenylmethyl phosphonate (kcat and kcat/Km values of 580 s(-1) and 1.2 × 10(5) M(-1) s(-1), respectively) and 4-nitrophenyl butyrate (kcat and kcat/Km values of 140 s(-1) and 1.4 × 10(5) M(-1) s(-1), respectively). Pmi1525 is stereoselective for the hydrolysis of chiral methylphosphonate esters. The enzyme hydrolyzes the (SP)-enantiomer of isobutyl 4-nitrophenyl methylphosphonate 14 times faster than the corresponding (RP)-enantiomer. The catalytic properties of this enzyme make it an attractive template for the evolution of novel enzymes for the detection, destruction, and detoxification of organophosphonate nerve agents.

The rate of protein evolution is determined by a combination of selective pressure on protein function and biophysical constraints on protein folding and structure. Determining the relative contributions of these properties is an unsolved problem in molecular evolution with broad implications for protein engineering and function prediction. As a case study, we examined the structural divergence of the rapidly evolving o-succinylbenzoate synthase (OSBS) family, which catalyzes a step in menaquinone synthesis in diverse microorganisms and plants. On average, the OSBS family is much more divergent than other protein families from the same set of species, with the most divergent family members sharing <15% sequence identity. Comparing 11 representative structures revealed that loss of quaternary structure and large deletions or insertions are associated with the family's rapid evolution. Neither of these properties has been investigated in previous studies to identify factors that affect the rate of protein evolution. Intriguingly, one subfamily retained a multimeric quaternary structure and has small insertions and deletions compared with related enzymes that catalyze diverse reactions. Many proteins in this subfamily catalyze both OSBS and N-succinylamino acid racemization (NSAR). Retention of ancestral structural characteristics in the NSAR/OSBS subfamily suggests that the rate of protein evolution is not proportional to the capacity to evolve new protein functions. Instead, structural features that are conserved among proteins with diverse functions might contribute to the evolution of new functions.

The continued increase in the size of the protein sequence databases as a result of advances in genome sequencing technology is overwhelming the ability to perform experimental characterization of function. Consequently, functions are assigned to the vast majority of proteins via automated, homology-based methods, with the result that as many as 50% are incorrectly annotated or unannotated ( Schnoes et al. PLoS Comput. Biol. 2009 , 5 ( 12 ), e1000605 ). This manuscript describes a study of the D-mannonate dehydratase (ManD) subgroup of the enolase superfamily (ENS) to investigate how function diverges as sequence diverges. Previously, one member of the subgroup had been experimentally characterized as ManD [dehydration of D-mannonate to 2-keto-3-deoxy-D-mannonate (equivalently, 2-keto-3-deoxy-D-gluconate)]. In this study, 42 additional members were characterized to sample sequence-function space in the ManD subgroup. These were found to differ in both catalytic efficiency and substrate specificity: (1) high efficiency (kcat/KM = 10(3) to 10(4) M(-1) s(-1)) for dehydration of D-mannonate, (2) low efficiency (kcat/KM = 10(1) to 10(2) M(-1) s(-1)) for dehydration of d-mannonate and/or D-gluconate, and 3) no-activity with either D-mannonate or D-gluconate (or any other acid sugar tested). Thus, the ManD subgroup is not isofunctional and includes D-gluconate dehydratases (GlcDs) that are divergent from the GlcDs that have been characterized in the mandelate racemase subgroup of the ENS (Lamble et al. FEBS Lett. 2004 , 576 , 133 - 136 ) (Ahmed et al. Biochem. J. 2005 , 390 , 529 - 540 ). These observations signal caution for functional assignment based on sequence homology and lay the foundation for the studies of the physiological functions of the GlcDs and the promiscuous ManDs/GlcDs.

The cytosolic glutathione transferase (cytGST) superfamily comprises more than 13,000 nonredundant sequences found throughout the biosphere. Their key roles in metabolism and defense against oxidative damage have led to thousands of studies over several decades. Despite this attention, little is known about the physiological reactions they catalyze and most of the substrates used to assay cytGSTs are synthetic compounds. A deeper understanding of relationships across the superfamily could provide new clues about their functions. To establish a foundation for expanded classification of cytGSTs, we generated similarity-based subgroupings for the entire superfamily. Using the resulting sequence similarity networks, we chose targets that broadly covered unknown functions and report here experimental results confirming GST-like activity for 82 of them, along with 37 new 3D structures determined for 27 targets. These new data, along with experimentally known GST reactions and structures reported in the literature, were painted onto the networks to generate a global view of their sequence-structure-function relationships. The results show how proteins of both known and unknown function relate to each other across the entire superfamily and reveal that the great majority of cytGSTs have not been experimentally characterized or annotated by canonical class. A mapping of taxonomic classes across the superfamily indicates that many taxa are represented in each subgroup and highlights challenges for classification of superfamily sequences into functionally relevant classes. Experimental determination of disulfide bond reductase activity in many diverse subgroups illustrate a theme common for many reaction types. Finally, sequence comparison between an enzyme that catalyzes a reductive dechlorination reaction relevant to bioremediation efforts with some of its closest homologs reveals differences among them likely to be associated with evolution of this unusual reaction. Interactive versions of the networks, associated with functional and other types of information, can be downloaded from the Structure-Function Linkage Database (SFLD; http://sfld.rbvi.ucsf.edu).