Alexandre P. Kuzin

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.

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.

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. Proteins 2009.

Quantum chemical 13Cα 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-α-helical protein, (PDB ID code 2JVD), is validated here by comparing the observed 13Cα 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 13Cα 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 13Cα 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 13Cα chemical shifts for a slightly longer construct of the protein solved by x-ray crystallography at 2.0-Å 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 13Cα chemical shifts in solution than the 3BHP crystal structure. In addition, the 13Cα-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.

Enzymes of vancomycin resistance: the structure of D-alanine-D-lactate ligase of naturally resistant Leuconostoc mesenteroides.

BACKGROUND The bacterial cell wall and the enzymes that synthesize it are targets of glycopeptide antibiotics (vancomycins and teicoplanins) and beta-lactams (penicillins and cephalosporins). Biosynthesis of cell wall peptidoglycan requires a crosslinking of peptidyl moieties on adjacent glycan strands. The D-alanine-D-alanine transpeptidase, which catalyzes this crosslinking, is the target of beta-lactam antibiotics. Glycopeptides, in contrast, do not inhibit an enzyme, but bind directly to D-alanine-D-alanine and prevent subsequent crosslinking by the transpeptidase. Clinical resistance to vancomycin in enterococcal pathogens has been traced to altered ligases producing D-alanine-D-lactate rather than D-alanine-D-alanine. RESULTS The structure of a D-alanine-D-lactate ligase has been determined by multiple anomalous dispersion (MAD) phasing to 2.4 A resolution. Co-crystallization of the Leuconostoc mesenteroides LmDdl2 ligase with ATP and a di-D-methylphosphinate produced ADP and a phosphinophosphate analog of the reaction intermediate of cell wall peptidoglycan biosynthesis. Comparison of this D-alanine-D-lactate ligase with the known structure of DdlB D-alanine-D-alanine ligase, a wild-type enzyme that does not provide vancomycin resistance, reveals alterations in the size and hydrophobicity of the site for D-lactate binding (subsite 2). A decrease was noted in the ability of the ligase to hydrogen bond a substrate molecule entering subsite 2. CONCLUSIONS Structural differences at subsite 2 of the D-alanine-D-lactate ligase help explain a substrate specificity shift (D-alanine to D-lactate) leading to remodeled cell wall peptidoglycan and vancomycin resistance in Gram-positive pathogens.

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.

Background: The complete structure is not known for any intermediate filament (IF) protein. Results: Linker 1 and rod 1B in human vimentin were characterized using electron paramagnetic resonance spectroscopy and x-ray crystallography. Conclusion: The rod 1B adopts two functional conformations that mediate formation of an anti-parallel “A11” tetramer. Significance: Understanding vimentin structure provides insight into all IFs and the related human pathologies. 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.

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

A challenge in the computational design of enzymes is that multiple properties must be simultaneously optimized -- substrate-binding, transition state stabilization, and product release -- 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 the designs could provide the basis for a new class of organophosphate captures agents.

Molecular insights into the biosynthesis of the f420 coenzyme.

Coenzyme F420, 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 F420-0 (F420 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 F420-producing organisms, and weak sequence homologs are also found in non-F420-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 Pi 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 F420 coenzyme. Large structural differences in the active site region of the non-F420-producing CofD homologs suggest that they catalyze a different biochemical reaction.

X-RAY STUDIES OF ENZYMES THAT INTERACT WITH PENICILLINS

Abstract. The technique of X-ray diffraction has been successfully applied to enzymes associated with peptidoglycan biosynthesis. The technique has taught us a great deal about the structures and catalytic mechanisms of penicillin-binding proteins and β-lactamases. An insight into the structural basis for antibiotic resistance is given.

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 Mn2+-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.