Design of activated serine-containing catalytic triads with atomic-level accuracy.

A challenge in the computational design of enzymes is that multiple properties, including substrate binding, transition state stabilization and product release, must be simultaneously optimized, and this has limited the absolute activity of successful designs. Here, we focus on a single critical property of many enzymes: the nucleophilicity of an active site residue that initiates catalysis. We design proteins with idealized serine-containing catalytic triads and assess their nucleophilicity directly in native biological systems using activity-based organophosphate probes. Crystal structures of the most successful designs show unprecedented agreement with computational models, including extensive hydrogen bonding networks between the catalytic triad (or quartet) residues, and mutagenesis experiments demonstrate that these networks are critical for serine activation and organophosphate reactivity. Following optimization by yeast display, the designs react with organophosphate probes at rates comparable to natural serine hydrolases. Co-crystal structures with diisopropyl fluorophosphate bound to the serine nucleophile suggest that the designs could provide the basis for a new class of organophosphate capture agents.

Structural and biochemical characterization of the bilin lyase CpcS from Thermosynechococcus elongatus.

Cyanobacterial phycobiliproteins have evolved to capture light energy over most of the visible spectrum due to their bilin chromophores, which are linear tetrapyrroles that have been covalently attached by enzymes called bilin lyases. We report here the crystal structure of a bilin lyase of the CpcS family from Thermosynechococcus elongatus (TeCpcS-III). TeCpcS-III is a 10-stranded ß barrel with two alpha helices and belongs to the lipocalin structural family. TeCpcS-III catalyzes both cognate as well as noncognate bilin attachment to a variety of phycobiliprotein subunits. TeCpcS-III ligates phycocyanobilin, phycoerythrobilin, and phytochromobilin to the alpha and beta subunits of allophycocyanin and to the beta subunit of phycocyanin at the Cys82-equivalent position in all cases. The active form of TeCpcS-III is a dimer, which is consistent with the structure observed in the crystal. With the use of the UnaG protein and its association with bilirubin as a guide, a model for the association between the native substrate, phycocyanobilin, and TeCpcS was produced.

The structure of vimentin linker 1 and rod 1B domains characterized by site-directed spin-labeling electron paramagnetic resonance (SDSL-EPR) and X-ray crystallography.

Despite the passage of ∼30 years since the complete primary sequence of the intermediate filament (IF) protein vimentin was reported, the structure remains unknown for both an individual protomer and the assembled filament. In this report, we present data describing the structure of vimentin linker 1 (L1) and rod 1B. Electron paramagnetic resonance spectra collected from samples bearing site-directed spin labels demonstrate that L1 is not a flexible segment between coiled-coils (CCs) but instead forms a rigid, tightly packed structure. An x-ray crystal structure of a construct containing L1 and rod 1B shows that it forms a tetramer comprising two equivalent parallel CC dimers that interact with one another in the form of a symmetrical anti-parallel dimer. Remarkably, the parallel CC dimers are themselves asymmetrical, which enables them to tetramerize rather than undergoing higher order oligomerization. This functionally vital asymmetry in the CC structure, encoded in the primary sequence of rod 1B, provides a striking example of evolutionary exploitation of the structural plasticity of proteins. EPR and crystallographic data consistently suggest that a very short region within L1 represents a minor local distortion in what is likely to be a continuous CC from the end of rod 1A through the entirety of rod 1B. The concordance of this structural model with previously published cross-linking and spectral data supports the conclusion that the crystallographic oligomer represents a native biological structure.

Computational design of catalytic dyads and oxyanion holes for ester hydrolysis.

Nucleophilic catalysis is a general strategy for accelerating ester and amide hydrolysis. In natural active sites, nucleophilic elements such as catalytic dyads and triads are usually paired with oxyanion holes for substrate activation, but it is difficult to parse out the independent contributions of these elements or to understand how they emerged in the course of evolution. Here we explore the minimal requirements for esterase activity by computationally designing artificial catalysts using catalytic dyads and oxyanion holes. We found much higher success rates using designed oxyanion holes formed by backbone NH groups rather than by side chains or bridging water molecules and obtained four active designs in different scaffolds by combining this motif with a Cys-His dyad. Following active site optimization, the most active of the variants exhibited a catalytic efficiency (k(cat)/K(M)) of 400 M(-1) s(-1) for the cleavage of a p-nitrophenyl ester. Kinetic experiments indicate that the active site cysteines are rapidly acylated as programmed by design, but the subsequent slow hydrolysis of the acyl-enzyme intermediate limits overall catalytic efficiency. Moreover, the Cys-His dyads are not properly formed in crystal structures of the designed enzymes. These results highlight the challenges that computational design must overcome to achieve high levels of activity.

Quantum chemical 13C(alpha) chemical shift calculations for protein NMR structure determination, refinement, and validation.

A recently determined set of 20 NMR-derived conformations of a 48-residue all-alpha-helical protein, (PDB ID code 2JVD), is validated here by comparing the observed (13)C(alpha) chemical shifts with those computed at the density functional level of theory. In addition, a recently introduced physics-based method, aimed at determining protein structures by using NOE-derived distance constraints together with observed and computed (13)C(alpha) chemical shifts, was applied to determine a new set of 10 conformations, (Set-bt), as a blind test for the same protein. A cross-validation of these two sets of conformations in terms of the agreement between computed and observed (13)C(alpha) chemical shifts, several stereochemical quality factors, and some NMR quality assessment scores reveals the good quality of both sets of structures. We also carried out an analysis of the agreement between the observed and computed (13)C(alpha) chemical shifts for a slightly longer construct of the protein solved by x-ray crystallography at 2.0-A resolution (PDB ID code 3BHP) with an identical amino acid residue sequence to the 2JVD structure for the first 46 residues. Our results reveal that both of the NMR-derived sets, namely 2JVD and Set-bt, are somewhat better representations of the observed (13)C(alpha) chemical shifts in solution than the 3BHP crystal structure. In addition, the (13)C(alpha)-based validation analysis appears to be more sensitive to subtle structural differences across the three sets of structures than any other NMR quality-assessment scores used here, and, although it is computationally intensive, this analysis has potential value as a standard procedure to determine, refine, and validate protein structures.

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study.

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.

Crystal structures of complexes between the R61 DD-peptidase and peptidoglycan-mimetic beta-lactams: a non-covalent complex with a "perfect penicillin".

The bacterial D-alanyl-D-alanine transpeptidases (DD-peptidases) are the killing targets of beta-lactams, the most important clinical defense against bacterial infections. However, due to the constant development of antibiotic-resistance mechanisms by bacteria, there is an ever-present need for new, more effective antimicrobial drugs. While enormous numbers of beta-lactam compounds have been tested for antibiotic activity in over 50 years of research, the success of a beta-lactam structure in terms of antibiotic activity remains unpredictable. Tipper and Strominger suggested long ago that beta-lactams inhibit DD-peptidases because they mimic the D-alanyl-D-alanine motif of the peptidoglycan substrate of these enzymes. They also predicted that beta-lactams having a peptidoglycan-mimetic side-chain might be better antibiotics than their non-specific counterparts, but decades of research have not provided any evidence for this. We have recently described two such novel beta-lactams. The first is a penicillin having the glycyl-L-alpha-amino-epsilon-pimelyl side-chain of Streptomyces strain R61 peptidoglycan, making it the "perfect penicillin" for this organism. The other is a cephalosporin with the same side-chain. Here, we describe the X-ray crystal structures of the perfect penicillin in non-covalent and covalent complexes with the Streptomyces R61 DD-peptidase. The structure of the non-covalent enzyme-inhibitor complex is the first such complex to be trapped crystallographically with a DD-peptidase. In addition, the covalent complex of the peptidyl-cephalosporin with the R61 DD-peptidase is described. Finally, two covalent complexes with the traditional beta-lactams benzylpenicillin and cephalosporin C were determined for comparison with the peptidyl beta-lactams. These structures, together with relevant kinetics data, support Tipper and Strominger's assertion that peptidoglycan-mimetic side-chains should improve beta-lactams as inhibitors of DD-peptidases.

A community resource of experimental data for NMR / X-ray crystal structure pairs.

We have developed an online NMR / X-ray Structure Pair Data Repository. The NIGMS Protein Structure Initiative (PSI) has provided many valuable reagents, 3D structures, and technologies for structural biology. The Northeast Structural Genomics Consortium was one of several PSI centers. NESG used both X-ray crystallography and NMR spectroscopy for protein structure determination. A key goal of the PSI was to provide experimental structures for at least one representative of each of hundreds of targeted protein domain families. In some cases, structures for identical (or nearly identical) constructs were determined by both NMR and X-ray crystallography. NMR spectroscopy and X-ray diffraction data for 41 of these "NMR / X-ray" structure pairs determined using conventional triple-resonance NMR methods with extensive sidechain resonance assignments have been organized in an online NMR / X-ray Structure Pair Data Repository. In addition, several NMR data sets for perdeuterated, methyl-protonated protein samples are included in this repository. As an example of the utility of this repository, these data were used to revisit questions about the precision and accuracy of protein NMR structures first outlined by Levy and coworkers several years ago (Andrec et al., Proteins 2007;69:449-465). These results demonstrate that the agreement between NMR and X-ray crystal structures is improved using modern methods of protein NMR spectroscopy. The NMR / X-ray Structure Pair Data Repository will provide a valuable resource for new computational NMR methods development.

The SufE sulfur-acceptor protein contains a conserved core structure that mediates interdomain interactions in a variety of redox protein complexes.

The isc and suf operons in Escherichia coli represent alternative genetic systems optimized to mediate the essential metabolic process of iron-sulfur cluster (Fe-S) assembly under basal or oxidative-stress conditions, respectively. Some of the proteins in these two operons share strong sequence homology, e.g. the cysteine desulfurases IscS and SufS, and presumably play the same role in the oxygen-sensitive assembly process. However, other proteins in these operons share no significant homology and occur in a mutually exclusive manner in Fe-S assembly operons in other organisms (e.g. IscU and SufE). These latter proteins presumably play distinct roles adapted to the different assembly mechanisms used by the two systems. IscU has three invariant cysteine residues that function as a template for Fe-S assembly while accepting a sulfur atom from IscS. SufE, in contrast, does not function as an Fe-S assembly template but has been suggested to function as a shuttle protein that uses a persulfide linkage to a single invariant cysteine residue to transfer a sulfur atom from SufS to an alternative Fe-S assembly template. Here, we present and analyze the 2.0A crystal structure of E.coli SufE. The structure shows that the persulfide-forming cysteine occurs at the tip of a loop with elevated B-factors, where its side-chain is buried from solvent exposure in a hydrophobic cavity located beneath a highly conserved surface. Despite the lack of sequence homology, the core of SufE shows strong structural similarity to IscU, and the sulfur-acceptor site in SufE coincides with the location of the cysteine residues mediating Fe-S cluster assembly in IscU. Thus, a conserved core structure is implicated in mediating the interactions of both SufE and IscU with the mutually homologous cysteine desulfurase enzymes present in their respective operons. A similar core structure is observed in a domain found in a variety of Fe-S cluster containing flavoenzymes including xanthine dehydrogenase, where it also mediates interdomain interactions. Therefore, the core fold of SufE/IscU has been adapted to mediate interdomain interactions in diverse redox protein systems in the course of evolution.

Precise assembly of complex beta sheet topologies from de novo designed building blocks.

Design of complex alpha-beta protein topologies poses a challenge because of the large number of alternative packing arrangements. A similar challenge presumably limited the emergence of large and complex protein topologies in evolution. Here, we demonstrate that protein topologies with six and seven-stranded beta sheets can be designed by insertion of one de novo designed beta sheet containing protein into another such that the two beta sheets are merged to form a single extended sheet, followed by amino acid sequence optimization at the newly formed strand-strand, strand-helix, and helix-helix interfaces. Crystal structures of two such designs closely match the computational design models. Searches for similar structures in the SCOP protein domain database yield only weak matches with different beta sheet connectivities. A similar beta sheet fusion mechanism may have contributed to the emergence of complex beta sheets during natural protein evolution.

Structure determination and biochemical characterization of a putative HNH endonuclease from Geobacter metallireducens GS-15.

The crystal structure of a putative HNH endonuclease, Gmet_0936 protein from Geobacter metallireducens GS-15, has been determined at 2.6 Å resolution using single-wavelength anomalous dispersion method. The structure contains a two-stranded anti-parallel ß-sheet that are surrounded by two helices on each face, and reveals a Zn ion bound in each monomer, coordinated by residues Cys38, Cys41, Cys73, and Cys76, which likely plays an important structural role in stabilizing the overall conformation. Structural homologs of Gmet_0936 include Hpy99I endonuclease, phage T4 endonuclease VII, and other HNH endonucleases, with these enzymes sharing 15-20% amino acid sequence identity. An overlay of Gmet_0936 and Hpy99I structures shows that most of the secondary structure elements, catalytic residues as well as the zinc binding site (zinc ribbon) are conserved. However, Gmet_0936 lacks the N-terminal domain of Hpy99I, which mediates DNA binding as well as dimerization. Purified Gmet_0936 forms dimers in solution and a dimer of the protein is observed in the crystal, but with a different mode of dimerization as compared to Hpy99I. Gmet_0936 and its N77H variant show a weak DNA binding activity in a DNA mobility shift assay and a weak Mn²âº-dependent nicking activity on supercoiled plasmids in low pH buffers. The preferred substrate appears to be acid and heat-treated DNA with AP sites, suggesting Gmet_0936 may be a DNA repair enzyme.

Crystal structures of Phd-Doc, HigA, and YeeU establish multiple evolutionary links between microbial growth-regulating toxin-antitoxin systems.

Bacterial toxin-antitoxin (TA) systems serve a variety of physiological functions including regulation of cell growth and maintenance of foreign genetic elements. Sequence analyses suggest that TA families are linked by complex evolutionary relationships reflecting likely swapping of functional domains between different TA families. Our crystal structures of Phd-Doc from bacteriophage P1, the HigA antitoxin from Escherichia coli CFT073, and YeeU of the YeeUWV systems from E. coli K12 and Shigella flexneri confirm this inference and reveal additional, unanticipated structural relationships. The growth-regulating Doc toxin exhibits structural similarity to secreted virulence factors that are toxic for eukaryotic target cells. The Phd antitoxin possesses the same fold as both the YefM and NE2111 antitoxins that inhibit structurally unrelated toxins. YeeU, which has an antitoxin-like activity that represses toxin expression, is structurally similar to the ribosome-interacting toxins YoeB and RelE. These observations suggest extensive functional exchanges have occurred between TA systems during bacterial evolution.

Molecular insights into the biosynthesis of the F420 coenzyme.

Coenzyme F(420), a hydride carrier, is found in Archaea and some bacteria and has crucial roles in methanogenesis, antibiotic biosynthesis, DNA repair, and activation of antitubercular compounds. CofD, 2-phospho-l-lactate transferase, catalyzes the last step in the biosynthesis of F(420)-0 (F(420) without polyglutamate), by transferring the lactyl phosphate moiety of lactyl(2)diphospho-(5')guanosine to 7,8-didemethyl-8-hydroxy-5-deazariboflavin ribitol (Fo). CofD is highly conserved among F(420)-producing organisms, and weak sequence homologs are also found in non-F(420)-producing organisms. This superfamily does not share any recognizable sequence conservation with other proteins. Here we report the first crystal structures of CofD, the free enzyme and two ternary complexes, with Fo and P(i) or with Fo and GDP, from Methanosarcina mazei. The active site is located at the C-terminal end of a Rossmann fold core, and three large insertions make significant contributions to the active site and dimer formation. The observed binding modes of Fo and GDP can explain known biochemical properties of CofD and are also supported by our binding assays. The structures provide significant molecular insights into the biosynthesis of the F(420) coenzyme. Large structural differences in the active site region of the non-F(420)-producing CofD homologs suggest that they catalyze a different biochemical reaction.

Functional insights from structural genomics.

Structural genomics efforts have produced structural information, either directly or by modeling, for thousands of proteins over the past few years. While many of these proteins have known functions, a large percentage of them have not been characterized at the functional level. The structural information has provided valuable functional insights on some of these proteins, through careful structural analyses, serendipity, and structure-guided functional screening. Some of the success stories based on structures solved at the Northeast Structural Genomics Consortium (NESG) are reported here. These include a novel methyl salicylate esterase with important role in plant innate immunity, a novel RNA methyltransferase (H. influenzae yggJ (HI0303)), a novel spermidine/spermine N-acetyltransferase (B. subtilis PaiA), a novel methyltransferase or AdoMet binding protein (A. fulgidus AF_0241), an ATP:cob(I)alamin adenosyltransferase (B. subtilis YvqK), a novel carboxysome pore (E. coli EutN), a proline racemase homolog with a disrupted active site (B. melitensis BME11586), an FMN-dependent enzyme (S. pneumoniae SP_1951), and a 12-stranded beta-barrel with a novel fold (V. parahaemolyticus VPA1032).

Structure and function of the phenazine biosynthetic protein PhzF from Pseudomonas fluorescens.

Phenazines produced by Pseudomonas and Streptomyces spp. are heterocyclic nitrogen-containing metabolites with antibiotic, antitumor, and antiparasitic activity. The antibiotic properties of pyocyanin, produced by Pseudomonas aeruginosa, were recognized in the 1890s, although this blue phenazine is now known to be a virulence factor in human disease. Despite their biological significance, the biosynthesis of phenazines is not fully understood. Here we present structural and functional studies of PhzF, an enzyme essential for phenazine synthesis in Pseudomonas spp. PhzF shares topology with diaminopimelate epimerase DapF but lacks the same catalytic residues. The structure of PhzF in complex with its substrate, trans-2,3-dihydro-3-hydroxyanthranilic acid, suggests that it is an isomerase using the conserved glutamate E45 to abstract a proton from C3 of the substrate. The proton is returned to C1 of the substrate after rearrangement of the double-bond system, yielding an enol that converts to the corresponding ketone. PhzF is a dimer that may be bifunctional, providing a shielded cavity for ketone dimerization via double Schiff-base formation to produce the phenazine scaffold. Our proposed mechanism is supported by mass and NMR spectroscopy. The results are discussed in the context of related structures and protein sequences of unknown biochemical function.

The 2.3-A crystal structure of the shikimate 5-dehydrogenase orthologue YdiB from Escherichia coli suggests a novel catalytic environment for an NAD-dependent dehydrogenase.

We present here the 2.3-A crystal structure of the Escherichia coli YdiB protein, an orthologue of shikimate 5-dehydrogenase. This enzyme catalyzes the reduction of 3-dehydroshikimate to shikimate as part of the shikimate pathway, which is absent in mammals but required for the de novo synthesis of aromatic amino acids, quinones, and folate in many other organisms. In this context, the shikimate pathway has been promoted as a target for the development of antimicrobial agents. The crystal structure of YdiB shows that the protomer contains two alpha/beta domains connected by two alpha-helices, with the N-terminal domain being novel and the C-terminal domain being a Rossmann fold. The NAD+ cofactor, which co-purified with the enzyme, is bound to the Rossmann domain in an elongated fashion with the nicotinamide ring in the pro-R conformation. Its binding site contains several unusual features, including a cysteine residue in close apposition to the nicotinamide ring and a clamp over the ribose of the adenosine moiety formed by phenylalanine and lysine residues. The structure explains the specificity for NAD versus NADP in different members of the shikimate dehydrogenase family on the basis of variations in the amino acid identity of several other residues in the vicinity of this ribose group. A cavity lined by residues that are 100% conserved among all shikimate dehydrogenases is found between the two domains of YdiB, in close proximity to the hydride acceptor site on the nicotinamide ring. Shikimate was modeled into this site in a geometry such that all of its heteroatoms form high quality hydrogen bonds with these invariant residues. Their strong conservation in all orthologues supports the possibility of developing broad spectrum inhibitors of this enzyme. The nature and disposition of the active site residues suggest a novel reaction mechanism in which an aspartate acts as the general acid/base catalyst during the hydride transfer reaction.