Rafael Toro

Large-Scale Determination of Sequence, Structure, and Function Relationships in Cytosolic Glutathione Transferases across the Biosphere

Global networks of the cytosolic glutathione S-transferases illuminate sequence-structure-function relationships across more than 13,000 members of this superfamily, including experimental confirmation of enzymatic activity for 82 members and new crystal structures for 27.

Homology models guide discovery of diverse enzyme specificities among dipeptide epimerases in the enolase superfamily

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.

Computation-facilitated assignment of the function in the enolase superfamily: a regiochemically distinct galactarate dehydratase from Oceanobacillus iheyensis .

The structure of an uncharacterized member of the enolase superfamily from Oceanobacillus iheyensis (GI 23100298, IMG locus tag Ob2843, PDB entry 2OQY ) was determined by the New York SGX Research Center for Structural Genomics (NYSGXRC). The structure contained two Mg(2+) ions located 10.4 A from one another, with one located in the canonical position in the (beta/alpha)(7)beta-barrel domain (although the ligand at the end of the fifth beta-strand is His, unprecedented in structurally characterized members of the superfamily); the second is located in a novel site within the capping domain. In silico docking of a library of mono- and diacid sugars to the active site predicted a diacid sugar as a likely substrate. Activity screening of a physical library of acid sugars identified galactarate as the substrate (k(cat) = 6.8 s(-1), K(M) = 620 microM, k(cat)/K(M) = 1.1 x 10(4) M(-1) s(-1)), allowing functional assignment of Ob2843 as galactarate dehydratase (GalrD-II). The structure of a complex of the catalytically impaired Y90F mutant with Mg(2+) and galactarate allowed identification of a Tyr 164-Arg 162 dyad as the base that initiates the reaction by abstraction of the alpha-proton and Tyr 90 as the acid that facilitates departure of the beta-OH leaving group. The enzyme product is 2-keto-d-threo-4,5-dihydroxyadipate, the enantiomer of the product obtained in the GalrD reaction catalyzed by a previously characterized bifunctional l-talarate/galactarate dehydratase (TalrD/GalrD). On the basis of the different active site structures and different regiochemistries, we recognize that these functions represent an example of apparent, not actual, convergent evolution of function. The structure of GalrD-II and its active site architecture allow identification of the seventh functionally and structurally characterized subgroup in the enolase superfamily. This study provides an additional example in which an integrated sequence- and structure-based strategy employing computational approaches is a viable approach for directing functional assignment of unknown enzymes discovered in genome projects.

Divergence of Structure and Function in the Haloacid Dehalogenase Enzyme Superfamily: Bacteroides thetaiotaomicron BT2127 is an Inorganic Pyrophosphatase

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.

Biochemical and Structural Characterization of the Human TL1A Ectodomain

TNF-like 1A (TL1A) is a newly described member of the TNF superfamily that is directly implicated in the pathogenesis of autoimmune diseases, including inflammatory bowel disease, atherosclerosis, and rheumatoid arthritis. We report the crystal structure of the human TL1A extracellular domain at a resolution of 2.5 A, which reveals a jelly-roll fold typical of the TNF superfamily. This structural information, in combination with complementary mutagenesis and biochemical characterization, provides insights into the binding interface and the specificity of the interactions between TL1A and the DcR3 and DR3 receptors. These studies suggest that the mode of interaction between TL1A and DcR3 differs from other characterized TNF ligand/receptor complexes. In addition, we have generated functional TL1A mutants with altered disulfide bonding capability that exhibit enhanced solution properties, which will facilitate the production of materials for future cell-based and whole animal studies. In summary, these studies provide insights into the structure and function of TL1A and provide the basis for the rational manipulation of its interactions with cognate receptors.

Structural Characterization of the Glycoprotein GP2 Core Domain from the CAS Virus, a Novel Arenavirus-Like Species

Abstract Fusion of the viral and host cell membranes is a necessary first step for infection by enveloped viruses and is mediated by the envelope glycoprotein. The transmembrane subunits from the structurally defined “class I” glycoproteins adopt an α-helical “trimer-of-hairpins” conformation during the fusion pathway. Here, we present our studies on the envelope glycoprotein transmembrane subunit, GP2, of the CAS virus (CASV). CASV was recently identified from annulated tree boas (Corallus annulatus) with inclusion body disease and is implicated in the disease etiology. We have generated and characterized two protein constructs consisting of the predicted CASV GP2 core domain. The crystal structure of the CASV GP2 post-fusion conformation indicates a trimeric α-helical bundle that is highly similar to those of Ebola virus and Marburg virus GP2 despite CASV genome homology to arenaviruses. Denaturation studies demonstrate that the stability of CASV GP2 is pH dependent with higher stability at lower pH; we propose that this behavior is due to a network of interactions among acidic residues that would destabilize the α-helical bundle under conditions where the side chains are deprotonated. The pH-dependent stability of the post-fusion structure has been observed in Ebola virus and Marburg virus GP2, as well as other viruses that enter via the endosome. Infection experiments with CASV and the related Golden Gate virus support a mechanism of entry that requires endosomal acidification. Our results suggest that, despite being primarily arenavirus like, the transmembrane subunit of CASV is extremely similar to the filoviruses.

Substrate Distortion and the Catalytic Reaction Mechanism of 5-Carboxyvanillate Decarboxylase

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 Mn2+ 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 × 104 M–1 s–1) and Novosphingobium aromaticivorans (kcat = 27 s–1 and kcat/Km = 1.1 × 105 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.

Loss of quaternary structure is associated with rapid sequence divergence in the OSBS family.

Significance The rate at which proteins accumulate amino acid substitutions during evolution depends on the likelihood that mutations will disrupt structure or affect function. Many mutations affect the ability of proteins to fold correctly, and previous studies showed that the burden imposed by misfolded proteins in cells heavily influences evolutionary rates of proteins. However, these studies could not examine the influence of function on evolutionary rates. The work described here examines the relationship between structural and functional divergence in a rapidly evolving protein family. This analysis revealed that family members that evolved a new function retained more ancestral sequence and structural characteristics, suggesting that the rate of protein evolution is not proportional to the capacity to evolve new functions. 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.

Crystal Structure of human Karyopherin β2 bound to the PY-NLS of Saccharomyces cerevisiae Nab2

Import-Karyopherin or Importin proteins bind nuclear localization signals (NLSs) to mediate the import of proteins into the cell nucleus. Karyopherin β2 or Kapβ2, also known as Transportin, is a member of this transporter family responsible for the import of numerous RNA binding proteins. Kapβ2 recognizes a targeting signal termed the PY-NLS that lies within its cargos to target them through the nuclear pore complex. The recognition of PY-NLS by Kapβ2 is conserved throughout eukaryotes. Kap104, the Kapβ2 homolog in Saccharomyces cerevisiae, recognizes PY-NLSs in cargos Nab2, Hrp1, and Tfg2. We have determined the crystal structure of Kapβ2 bound to the PY-NLS of the mRNA processing protein Nab2 at 3.05-Å resolution. A seven-residue segment of the PY-NLS of Nab2 is observed to bind Kapβ2 in an extended conformation and occupies the same PY-NLS binding site observed in other Kapβ2·PY-NLS structures.

Function Discovery and Structural Characterization of a Methylphosphonate Esterase

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.