Peter M. U. Ung

Automated genome mining for natural products

BackgroundDiscovery of new medicinal agents from natural sources has largely been an adventitious process based on screening of plant and microbial extracts combined with bioassay-guided identification and natural product structure elucidation. Increasingly rapid and more cost-effective genome sequencing technologies coupled with advanced computational power have converged to transform this trend toward a more rational and predictive pursuit.ResultsWe have developed a rapid method of scanning genome sequences for multiple polyketide, nonribosomal peptide, and mixed combination natural products with output in a text format that can be readily converted to two and three dimensional structures using conventional software. Our open-source and web-based program can assemble various small molecules composed of twenty standard amino acids and twenty two other chain-elongation intermediates used in nonribosomal peptide systems, and four acyl-CoA extender units incorporated into polyketides by reading a hidden Markov model of DNA. This process evaluates and selects the substrate specificities along the assembly line of nonribosomal synthetases and modular polyketide synthases.ConclusionUsing this approach we have predicted the structures of natural products from a diverse range of bacteria based on a limited number of signature sequences. In accelerating direct DNA to metabolomic analysis, this method bridges the interface between chemists and biologists and enables rapid scanning for compounds with potential therapeutic value.

CSAR Benchmark Exercise of 2010: Combined Evaluation Across All Submitted Scoring Functions

As part of the Community Structure-Activity Resource (CSAR) center, a set of 343 high-quality, protein–ligand crystal structures were assembled with experimentally determined Kd or Ki information from the literature. We encouraged the community to score the crystallographic poses of the complexes by any method of their choice. The goal of the exercise was to (1) evaluate the current ability of the field to predict activity from structure and (2) investigate the properties of the complexes and methods that appear to hinder scoring. A total of 19 different methods were submitted with numerous parameter variations for a total of 64 sets of scores from 16 participating groups. Linear regression and nonparametric tests were used to correlate scores to the experimental values. Correlation to experiment for the various methods ranged R2 = 0.58–0.12, Spearman ρ = 0.74–0.37, Kendall τ = 0.55–0.25, and median unsigned error = 1.00–1.68 pKd units. All types of scoring functions—force field based, knowledge based, and empirical—had examples with high and low correlation, showing no bias/advantage for any particular approach. The data across all the participants were combined to identify 63 complexes that were poorly scored across the majority of the scoring methods and 123 complexes that were scored well across the majority. The two sets were compared using a Wilcoxon rank-sum test to assess any significant difference in the distributions of >400 physicochemical properties of the ligands and the proteins. Poorly scored complexes were found to have ligands that were the same size as those in well-scored complexes, but hydrogen bonding and torsional strain were significantly different. These comparisons point to a need for CSAR to develop data sets of congeneric series with a range of hydrogen-bonding and hydrophobic characteristics and a range of rotatable bonds.

Binding of a Small Molecule at a Protein–Protein Interface Regulates the Chaperone Activity of Hsp70–Hsp40

Heat shock protein 70 (Hsp70) is a highly conserved molecular chaperone that plays multiple roles in protein homeostasis. In these various tasks, the activity of Hsp70 is shaped by interactions with co-chaperones, such as Hsp40. The Hsp40 family of co-chaperones binds to Hsp70 through a conserved J-domain, and these factors stimulate ATPase and protein-folding activity. Using chemical screens, we identified a compound, 115-7c, which acts as an artificial co-chaperone for Hsp70. Specifically, the activities of 115-7c mirrored those of a Hsp40; the compound stimulated the ATPase and protein-folding activities of a prokaryotic Hsp70 (DnaK) and partially compensated for a Hsp40 loss-of-function mutation in yeast. Consistent with these observations, NMR and mutagenesis studies indicate that the binding site for 115-7c is adjacent to a region on DnaK that is required for J-domain-mediated stimulation. Interestingly, we found that 115-7c and the Hsp40 do not compete for binding but act in concert. Using this information, we introduced additional steric bulk to 115-7c and converted it into an inhibitor. Thus, these chemical probes either promote or inhibit chaperone functions by regulating Hsp70-Hsp40 complex assembly at a native protein-protein interface. This unexpected mechanism may provide new avenues for exploring how chaperones and co-chaperones cooperate to shape protein homeostasis.

CSAR Benchmark Exercise of 2010: Selection of the Protein–Ligand Complexes

A major goal in drug design is the improvement of computational methods for docking and scoring. The Community Structure Activity Resource (CSAR) aims to collect available data from industry and academia which may be used for this purpose (www.csardock.org). Also, CSAR is charged with organizing community-wide exercises based on the collected data. The first of these exercises was aimed to gauge the overall state of docking and scoring, using a large and diverse data set of protein–ligand complexes. Participants were asked to calculate the affinity of the complexes as provided and then recalculate with changes which may improve their specific method. This first data set was selected from existing PDB entries which had binding data (Kd or Ki) in Binding MOAD, augmented with entries from PDBbind. The final data set contains 343 diverse protein–ligand complexes and spans 14 pKd. Sixteen proteins have three or more complexes in the data set, from which a user could start an inspection of congeneric series. Inherent experimental error limits the possible correlation between scores and measured affinity; R2 is limited to ∼0.9 when fitting to the data set without over parametrizing. R2 is limited to ∼0.8 when scoring the data set with a method trained on outside data. The details of how the data set was initially selected, and the process by which it matured to better fit the needs of the community are presented. Many groups generously participated in improving the data set, and this underscores the value of a supportive, collaborative effort in moving our field forward.

Mutagenesis Reveals the Complex Relationships between ATPase Rate and the Chaperone Activities of Escherichia coli Heat Shock Protein 70 (Hsp70/DnaK)

The Escherichia coli 70-kDa heat shock protein, DnaK, is a molecular chaperone that engages in a variety of cellular activities, including the folding of proteins. During this process, DnaK binds its substrates in coordination with a catalytic ATPase cycle. Both the ATPase and protein folding activities of DnaK are stimulated by its co-chaperones, DnaJ and GrpE. However, it is not yet clear how changes in the stimulated ATPase rate of DnaK impact the folding process. In this study, we performed mutagenesis throughout the nucleotide-binding domain of DnaK to generate a collection of mutants in which the stimulated ATPase rates varied from 0.7 to 13.6 pmol/μg/min−1. We found that this range was largely established by differences in the ability of the mutants to be stimulated by one or both of the co-chaperones. Next, we explored how changes in ATPase rate might impact refolding of denatured luciferase in vitro and found that the two activities were poorly correlated. Unexpectedly, we found several mutants that refold luciferase normally in the absence of significant ATP turnover, presumably by increasing the flexibility of DnaK. Finally, we tested whether DnaK mutants could complement growth of ΔdnaK E. coli cells under heat shock and found that the ability to refold luciferase was more predictive of in vivo activity than ATPase rate. This study provides insights into how flexibility and co-chaperone interactions affect DnaK-mediated ATP turnover and protein folding.

A poke in the eye: Inhibiting HIV-1 protease through its flap-recognition pocket

A novel mechanism of inhibiting HIV-1 protease (HIVp) is presented. Using computational solvent mapping to identify complementary interactions and the Multiple Protein Structure method to incorporate protein flexibility, we generated a receptor-based pharmacophore model of the flexible flap region of the semiopen, apo state of HIVp. Complementary interactions were consistently observed at the base of the flap, only within a cleft with a specific structural role. In the closed, bound state of HIVp, each flap tip docks against the opposite monomer, occupying this cleft. This flap-recognition site is filled by the protein and cannot be identified using traditional approaches based on bound, closed structures. Virtual screening and dynamics simulations show how small molecules can be identified to complement this cleft. Subsequent experimental testing confirms inhibitory activity of this new class of inhibitor. This may be the first new inhibitor class for HIVp since dimerization inhibitors were introduced 17 years ago.

Chemical Probes That Selectively Recognize the Earliest Aβ Oligomers in Complex Mixtures

Alzheimers disease (AD) is characterized by the self-assembly of amyloid beta (Aβ) peptides. Recent models implicate some of the earliest Aβ oligomers, such as trimers and tetramers, in disease. However, the roles of these structures remain uncertain, in part, because selective probes of their formation are not available. Toward that goal, we generated bivalent versions of the known Aβ ligand, the pentapeptide KLVFF. We found that compounds containing sufficiently long linkers (∼19 to 24 Å) recognized primarily Aβ trimers and tetramers, with little binding to either monomer or higher order structures. These compounds might be useful probes for early Aβ oligomers.

Identifying binding hot spots on protein surfaces by mixed‐solvent molecular dynamics: HIV‐1 protease as a test case

Mixed‐solvent molecular dynamics (MixMD) simulations use full protein flexibility and competition between water and small organic probes to achieve accurate hot‐spot mapping on protein surfaces. In this study, we improved MixMD using human immunodeficiency virus type‐1 protease (HIVp) as the test case. We used three probe–water solutions (acetonitrile–water, isopropanol–water, and pyrimidine–water), first at 50% w/w concentration and later at 5% v/v. Paradoxically, better mapping was achieved by using fewer probes; 5% simulations gave a superior signal‐to‐noise ratio and far fewer spurious hot spots than 50% MixMD. Furthermore, very intense and well‐defined probe occupancies were observed in the catalytic site and potential allosteric sites that have been confirmed experimentally. The Eye site, an allosteric site underneath the flap of HIVp, has been confirmed by the presence of a 5‐nitroindole fragment in a crystal structure. MixMD also mapped two additional hot spots: the Exo site (between the Gly16‐Gly17 and Cys67‐Gly68 loops) and the Face site (between Glu21‐Ala22 and Val84‐Ile85 loops). The Exo site was observed to overlap with crystallographic additives such as acetate and dimethyl sulfoxide that are present in different crystal forms of the protein. Analysis of crystal structures of HIVp in different symmetry groups has shown that some surface sites are common interfaces for crystal contacts, which means that they are surfaces that are relatively easy to desolvate and complement with organic molecules. MixMD should identify these sites; in fact, their occupancy values help establish a solid cut‐off where “druggable” sites are required to have higher occupancies than the crystal‐packing faces.

An allosteric modulator of HIV-1 protease shows equipotent inhibition of wild-type and drug-resistant proteases.

NMR and MD simulations have demonstrated that the flaps of HIV-1 protease (HIV-1p) adopt a range of conformations that are coupled with its enzymatic activity. Previously, a model was created for an allosteric site located between the flap and the core of HIV-1p, called the Eye site (Biopolymers2008, 89, 643−65218381626). Here, results from our first study were combined with a ligand-based, lead-hopping method to identify a novel compound (NIT). NIT inhibits HIV-1p, independent of the presence of an active-site inhibitor such as pepstatin A. Assays showed that NIT acts on an allosteric site other than the dimerization interface. MD simulations of the ligand–protein complex show that NIT stably binds in the Eye site and restricts the flaps. That bound state of NIT is consistent with a crystal structure of similar fragments bound in the Eye site (Chem. Biol. Drug Des.2010, 75, 257−26820659109). Most importantly, NIT is equally potent against wild-type and a multidrug-resistant mutant of HIV-1p, which highlights the promise of allosteric inhibitors circumventing existing clinical resistance.

Identification of Key Hinge Residues Important for Nucleotide-Dependent Allostery in E. coli Hsp70/DnaK

DnaK is a molecular chaperone that has important roles in protein folding. The hydrolysis of ATP is essential to this activity, and the effects of nucleotides on the structure and function of DnaK have been extensively studied. However, the key residues that govern the conformational motions that define the apo, ATP-bound, and ADP-bound states are not entirely clear. Here, we used molecular dynamics simulations, mutagenesis, and enzymatic assays to explore the molecular basis of this process. Simulations of DnaKs nucleotide-binding domain (NBD) in the apo, ATP-bound, and ADP/Pi-bound states suggested that each state has a distinct conformation, consistent with available biochemical and structural information. The simulations further suggested that large shearing motions between subdomains I-A and II-A dominated the conversion between these conformations. We found that several evolutionally conserved residues, especially G228 and G229, appeared to function as a hinge for these motions, because they predominantly populated two distinct states depending on whether ATP or ADP/Pi was bound. Consistent with the importance of these “hinge” residues, alanine point mutations caused DnaK to have reduced chaperone activities in vitro and in vivo. Together, these results clarify how sub-domain motions communicate allostery in DnaK.