Small-molecule inhibitors of the pseudaminic acid biosynthetic pathway: targeting motility as a key bacterial virulence factor.

Helicobacter pylori is motile by means of polar flagella, and this motility has been shown to play a critical role in pathogenicity. The major structural flagellin proteins have been shown to be glycosylated with the nonulosonate sugar, pseudaminic acid (Pse). This glycan is unique to microorganisms, and the process of flagellin glycosylation is required for H. pylori flagellar assembly and consequent motility. As such, the Pse biosynthetic pathway offers considerable potential as an antivirulence drug target, especially since motility is required for H. pylori colonization and persistence in the host. This report describes screening the five Pse biosynthetic enzymes for small-molecule inhibitors using both high-throughput screening (HTS) and in silico (virtual screening [VS]) approaches. Using a 100,000-compound library, 1,773 hits that exhibited a 40% threshold inhibition at a 10 µM concentration were identified by HTS. In addition, VS efforts using a 1.6-million compound library directed at two pathway enzymes identified 80 hits, 4 of which exhibited reasonable inhibition at a 10 µM concentration in vitro. Further secondary screening which identified 320 unique molecular structures or validated hits was performed. Following kinetic studies and structure-activity relationship (SAR) analysis of selected inhibitors from our refined list of 320 compounds, we demonstrated that three inhibitors with 50% inhibitory concentrations (IC50s) of approximately 14 µM, which belonged to a distinct chemical cluster, were able to penetrate the Gram-negative cell membrane and prevent formation of flagella.

Structural and mechanistic insight into covalent substrate binding by Escherichia coli dihydroxyacetone kinase.

The Escherichia coli dihydroxyacetone (Dha) kinase is an unusual kinase because (i) it uses the phosphoenolpyruvate carbohydrate: phosphotransferase system (PTS) as the source of high-energy phosphate, (ii) the active site is formed by two subunits, and (iii) the substrate is covalently bound to His218(K)* of the DhaK subunit. The PTS transfers phosphate to DhaM, which in turn phosphorylates the permanently bound ADP coenzyme of DhaL. This phosphoryl group is subsequently transferred to the Dha substrate bound to DhaK. Here we report the crystal structure of the E. coli Dha kinase complex, DhaK-DhaL. The structure of the complex reveals that DhaK undergoes significant conformational changes to accommodate binding of DhaL. Combined mutagenesis and enzymatic activity studies of kinase mutants allow us to propose a catalytic mechanism for covalent Dha binding, phosphorylation, and release of the Dha-phosphate product. Our results show that His56(K) is involved in formation of the covalent hemiaminal bond with Dha. The structure of H56N(K) with noncovalently bound substrate reveals a somewhat different positioning of Dha in the binding pocket as compared to covalently bound Dha, showing that the covalent attachment to His218(K) orients the substrate optimally for phosphoryl transfer. Asp109(K) is critical for activity, likely acting as a general base activating the γ-OH of Dha. Our results provide a comprehensive picture of the roles of the highly conserved active site residues of dihydroxyacetone kinases.

Contribution of active site residues to substrate hydrolysis by USP2: insights into catalysis by ubiquitin specific proteases.

The ubiquitin-specific protease (USP) structural class represents the largest and most diverse family of deubiquitinating enzymes (DUBs). Many USPs assume important biological roles and emerge as potential targets for therapeutic intervention. A clear understanding of USP catalytic mechanism requires a functional evaluation of the proposed key active site residues. Crystallographic data of ubiquitin aldehyde adducts of USP catalytic cores provided structural details on the catalytic triad residues, namely the conserved Cys and His, and a variable putative third residue, and inferred indirect structural roles for two other conserved residues (Asn and Asp), in stabilizing via a bridging water molecule the oxyanion of the tetrahedral intermediate (TI). We have expressed the catalytic domain of USP2 and probed by site-directed mutagenesis the role of these active site residues in the hydrolysis of peptide and isopeptide substrates, including a synthetic K48-linked diubiquitin substrate for which a label-free, mass spectrometry based assay has been developed to monitor cleavage. Hydrolysis of ubiquitin-AMC, a model substrate, was not affected by the mutations. Molecular dynamics simulations of USP2, free and complexed with the TI of a bona fide isopeptide substrate, were carried out. We found that Asn271 is structurally poised to directly stabilize the oxyanion developed in the acylation step, while being structurally supported by the adjacent absolutely conserved Asp575. Mutagenesis data functionally confirmed this structural role independent of the nature (isopeptide vs peptide) of the bond being cleaved. We also found that Asn574, structurally located as the third member of the catalytic triad, does not fulfill this role functionally. A dual supporting role is inferred from double-point mutation and structural data for the absolutely conserved residue Asp575, in oxyanion hole formation, and in maintaining the correct alignment and protonation of His557 for catalytic competency.

Catalytic mechanism of heparinase II investigated by site-directed mutagenesis and the crystal structure with its substrate.

Heparinase II (HepII) is an 85-kDa dimeric enzyme that depolymerizes both heparin and heparan sulfate glycosaminoglycans through a beta-elimination mechanism. Recently, we determined the crystal structure of HepII from Pedobacter heparinus (previously known as Flavobacterium heparinum) in complex with a heparin disaccharide product, and identified the location of its active site. Here we present the structure of HepII complexed with a heparan sulfate disaccharide product, proving that the same binding/active site is responsible for the degradation of both uronic acid epimers containing substrates. The key enzymatic step involves removal of a proton from the C5 carbon (a chiral center) of the uronic acid, posing a topological challenge to abstract the proton from either side of the ring in a single active site. We have identified three potential active site residues equidistant from C5 and located on both sides of the uronate product and determined their role in catalysis using a set of defined tetrasaccharide substrates. HepII H202A/Y257A mutant lost activity for both substrates and we determined its crystal structure complexed with a heparan sulfate-derived tetrasaccharide. Based on kinetic characterization of various mutants and the structure of the enzyme-substrate complex we propose residues participating in catalysis and their specific roles.

Molecular dynamics study of small molecule inhibitors of the Bcl-2 family.

We carried out docking and molecular dynamics simulations on ABT-737 and obatoclax, which are inhibitors of the Bcl-2 family of proteins. We modeled the binding mode of ABT-737 with Bcl-x(L) , Bcl-2, and Mcl-1 and examined their dynamical behavior. We found that the binding of the chlorobiphenyl end of ABT-737 was quite stable across all three proteins. However, the phenylpiperazine linker group was dramatically more mobile in Mcl-1 compared to either Bcl-x(L) or Bcl-2. The S-phenyl group at the p4 binding site was well-anchored in Bcl-x(L) and Bcl-2 but was somewhat more mobile in Mcl-1 although the phenyl ring itself on average stayed close to the p4 binding site in Mcl-1. This greater mobility is likely due to the greater openness of the p3 and p4 binding sites on Mcl-1. The calculated binding free energies were consistent with the much weaker binding affinity of ABT-737 for Mcl-1. Obatoclax was predicted to bind at the p1 and p2 binding sites of Mcl-1 and the binding mode was quite stable during the molecular dynamics simulation with Mcl-1 wrapping around the molecule. The modeled binding mode suggests that obatoclax is able to inhibit all three proteins because it makes use of the p1 and p2 binding sites alone, which is a fairly narrow groove in all three proteins unlike the p4 binding site, which is much broader in Mcl-1.

Solvated interaction energy (SIE) for scoring protein-ligand binding affinities. 2. Benchmark in the CSAR-2010 scoring exercise.

Solvated interaction energy (SIE) is an end-point physics-based scoring function for predicting binding affinities from force-field nonbonded interaction terms, continuum solvation, and configurational entropy linear compensation. We tested the SIE function in the Community Structure-Activity Resource (CSAR) scoring challenge consisting of high-resolution cocrystal structures for 343 protein-ligand complexes with high-quality binding affinity data and high diversity with respect to protein targets. Particular emphasis was placed on the sensitivity of SIE predictions to the assignment of protonation and tautomeric states in the complex and the treatment of metal ions near the protein-ligand interface. These were manually curated from an originally distributed CSAR-HiQ data set version, leading to the currently distributed CSAR-NRC-HiQ version. We found that this manual curation was a critical step for accurately testing the performance of the SIE function. The standard SIE parametrization, previously calibrated on an independent data set, predicted absolute binding affinities with a mean-unsigned-error (MUE) of 2.41 kcal/mol for the CSAR-HiQ version, which improved to 1.98 kcal/mol for the upgraded CSAR-NRC-HiQ version. Half-half retraining-testing of SIE parameters on two predefined subsets of CSAR-NRC-HiQ led to only marginal further improvements to an MUE of 1.83 kcal/mol. Hence, we do not recommend altering the current default parameters of SIE at this time. For a sample of SIE outliers, additional calculations by molecular dynamics-based SIE averaging with or without incorporation of ligand strain, by MM-PB(GB)/SA methods with or without entropic estimates, or even by the linear interaction energy (LIE) formalism with an explicit solvent model, did not further improve predictions.

Structural and functional analysis of Campylobacter jejuni PseG: a udp-sugar hydrolase from the pseudaminic acid biosynthetic pathway.

Flagella of the bacteria Helicobacter pylori and Campylobacter jejuni are important virulence determinants, whose proper assembly and function are dependent upon glycosylation at multiple positions by sialic acid-like sugars, such as 5,7-diacetamido-3,5,7,9-tetradeoxy-l-glycero-l-manno-nonulosonic acid (pseudaminic acid (Pse)). The fourth enzymatic step in the pseudaminic acid pathway, the hydrolysis of UDP-2,4-diacetamido-2,4,6-trideoxy-beta-l-altropyranose to generate 2,4-diacetamido-2,4,6-trideoxy-l-altropyranose, is performed by the nucleotide sugar hydrolase PseG. To better understand the molecular basis of the PseG catalytic reaction, we have determined the crystal structures of C. jejuni PseG in apo-form and as a complex with its UDP product at 1.8 and 1.85 A resolution, respectively. In addition, molecular modeling was utilized to provide insight into the structure of the PseG-substrate complex. This modeling identifies a His(17)-coordinated water molecule as the putative nucleophile and suggests the UDP-sugar substrate adopts a twist-boat conformation upon binding to PseG, enhancing the exposure of the anomeric bond cleaved and favoring inversion at C-1. Furthermore, based on these structures a series of amino acid substitution derivatives were constructed, altering residues within the active site, and each was kinetically characterized to examine its contribution to PseG catalysis. In conjunction with structural comparisons, the almost complete inactivation of the PseG H17F and H17L derivatives suggests that His(17) functions as an active site base, thereby activating the nucleophilic water molecule for attack of the anomeric C-O bond of the UDP-sugar. As the PseG structure reveals similarity to those of glycosyltransferase family-28 members, in particular that of Escherichia coli MurG, these findings may also be of relevance for the mechanistic understanding of this important enzyme family.

Trapping open and closed forms of FitE: a group III periplasmic binding protein.

Periplasmic binding proteins (PBPs) are essential components of bacterial transport systems, necessary for bacterial growth and survival. The two-domain structures of PBPs are topologically classified into three groups based on the number of crossovers or hinges between the globular domains: group I PBPs have three connections, group II have two, and group III have only one. Although a large number of structures for group I or II PBPs are known, fewer group III PBPs have been structurally characterized. Group I and II PBPs exhibit significant domain motions during transition from the unbound to ligand-bound form, however, no large conformational changes have been observed to date in group III PBPs. We have solved the crystal structure of a periplasmic binding protein FitE, part of an iron transport system, fit, recently identified in a clinical E. coli isolate. The structure, determined at 1.8 A resolution, shows that FitE is a group III PBP containing a single alpha-helix bridging the two domains. Among the individual FitE molecules present in two crystal forms we observed three different conformations (open, closed, intermediate). Our crystallographic and molecular dynamics results strongly support the notion that group III PBPs also adopt the same Venus flytrap mechanism as do groups I and II PBPs. Unlike other group III PBPs, FitE forms dimers both in solution and in the crystals. The putative siderophore binding pocket is lined with arginine residues, suggesting an anionic nature of the iron-containing siderophore.

Structure of L-xylulose-5-Phosphate 3-epimerase (UlaE) from the anaerobic L-ascorbate utilization pathway of Escherichia coli: identification of a novel phosphate binding motif within a TIM barrel fold.

Three catabolic enzymes, UlaD, UlaE, and UlaF, are involved in a pathway leading to fermentation of l-ascorbate under anaerobic conditions. UlaD catalyzes a beta-keto acid decarboxylation reaction to produce L-xylulose-5-phosphate, which undergoes successive epimerization reactions with UlaE (L-xylulose-5-phosphate 3-epimerase) and UlaF (L-ribulose-5-phosphate 4-epimerase), yielding D-xylulose-5-phosphate, an intermediate in the pentose phosphate pathway. We describe here crystallographic studies of UlaE from Escherichia coli O157:H7 that complete the structural characterization of this pathway. UlaE has a triosephosphate isomerase (TIM) barrel fold and forms dimers. The active site is located at the C-terminal ends of the parallel beta-strands. The enzyme binds Zn(2+), which is coordinated by Glu155, Asp185, His211, and Glu251. We identified a phosphate-binding site formed by residues from the beta1/alpha1 loop and alpha3' helix in the N-terminal region. This site differs from the well-characterized phosphate-binding motif found in several TIM barrel superfamilies that is located at strands beta7 and beta8. The intrinsic flexibility of the active site region is reflected by two different conformations of loops forming part of the substrate-binding site. Based on computational docking of the L-xylulose 5-phosphate substrate to UlaE and structural similarities of the active site of this enzyme to the active sites of other epimerases, a metal-dependent epimerization mechanism for UlaE is proposed, and Glu155 and Glu251 are implicated as catalytic residues. Mutation and activity measurements for structurally equivalent residues in related epimerases supported this mechanistic proposal.

Structure and function of the glycopeptide N-methyltransferase MtfA, a tool for the biosynthesis of modified glycopeptide antibiotics.

There is a considerable interest in the modification of existing antibiotics to generate new antimicrobials. Glycopeptide antibiotics (GPAs) are effective against serious Gram-positive bacterial pathogens including methicillin-resistant Staphylococcus aureus. However, resistance to these antibiotics is becoming a serious problem requiring new strategies. We show that the Amycolatopsis orientalis (S)-adenosyl-L-methionine-dependent methyltransferase MtfA, from the vancomycin-class GPA chloroeremomycin biosynthetic pathway, catalyzes in vivo and in vitro methyl transfer to generate methylated GPA derivatives of the teicoplanin class. The crystal structure of MtfA complexed with (S)-adenosyl-L-methionine, (S)-adenosylhomocysteine, or sinefungin inhibitor, coupled with mutagenesis, identified His228 as a likely general base required for methyl transfer to the N terminus of the glycopeptide. Computational docking and molecular dynamics simulations were used to model binding of demethyl-vancomycin aglycone to MtfA. These results demonstrate its utility as a tool for engineering methylated analogs of GPAs.

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex.

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.

High-resolution field-cycling NMR studies of a DNA octamer as a probe of phosphodiester dynamics and comparison with computer simulation.

Phosphorus-spin longitudinal relaxation rates of the DNA duplex octamer [d(GGAATTCC)](2) have been measured from 0.1 to 17.6 T by means of conventional and new field-cycling NMR methods. The high-resolution field-cycling method is identical to a conventional relaxation experiment, except that after preparation the sample is moved pneumatically from its usual position at the center of the high-resolution magnet upward to a lower field above its normal position and then returned to the center for readout after it has relaxed for the programmed relaxation delay at the low field. This is the first measurement of all longitudinal relaxation rates R(1) of a nuclear species in a macromolecule over virtually the entire accessible magnetic field range. For detailed analysis, three magnetic field regions can be delineated: (i) dipolar relaxation dominates at fields below 2 T, (ii) chemical shift anisotropy (CSA) relaxation is roughly constant from 2 to 6 T, and (iii) a square-law increasing dependence is seen at fields higher than approximately 6 T due to internal motion CSA relaxation. The analysis provides a rotational correlation time (tau(r) = 4.1 +/- 0.3 ns) for the duplex at both 1.5 and 0.25 mM concentrations (of duplex) at 22 degrees C. For comparison, extraction of tau(r) in the conventional way from the ratio of T(1)/T(2) at 14 T yields 3.2 ns. The tau(r) discrepancy disappears when we exclude the contribution of internal motion from the R(1) in the ratio. The low-field dipolar relaxation provides a weighted inverse sixth power sum of the distances from the phosphorus to the protons responsible for relaxation. This average is similar for all phosphates in the octamer and similar to that in previous B-DNA structures (its inverse sixth root is about 2.40 A for two different concentrations of octamer). The CSA relaxation at intermediate field provides an estimate of the order parameter squared, S(c)(2), for each phosphorus. S(c)(2) is about 0.7-1, clearly different for different phosphate linkages in the octamer duplex. The increasing R(1) at high fields reflects CSA relaxation due to internal motions, for which a correlation time, tau(hf), can be approximately extracted with the aid of additional measurements at 14.0 and 17.6 T. We conclude that tau(hf) values are relatively large, in the range of about 150 ps. Insight into the motions leading to this correlation time was gained by a 28 ns molecular dynamics simulation of the molecule. S(2) and tau(s) (corresponding to tau(hf)) predicted by this simulation were in good agreement with the experimental values from the field-cycling data. Both the effect of Mg(2+) on the dynamic parameters extracted from (31)P relaxation rates and the field dependence of relaxation rates for several protons of the octamer were measured. High-resolution field cycling opens up the possibility of monitoring residue-specific dipolar interactions and dynamics for the phosphorus nuclei of diverse oligonucleotides.