Iminoguanidines as Allosteric Inhibitors of the Iron-Regulated Heme Oxygenase (HemO) of Pseudomonas aeruginosa

Site Identification by Ligand Competitive Saturation-Biologics Approach for Structure-Based Protein Charge Prediction

Protein-based therapeutics typically require high concentrations of the active protein, which can lead to protein aggregation and high solution viscosity. Such solution behaviors can limit the stability, bioavailability, and manufacturability of protein-based therapeutics and are directly influenced by the charge of a protein. Protein charge is a system property affected by its environment, including the buffer composition, pH, and temperature. Thus, the charge calculated by summing the charges of each residue in a protein, as is commonly done in computational methods, may significantly differ from the effective charge of the protein as these calculations do not account for contributions from bound ions. Here, we present an extension of the structure-based approach termed site identification by ligand competitive saturation-biologics (SILCS-Biologics) to predict the effective charge of proteins. The SILCS-Biologics approach was applied on a range of protein targets in different salt environments for which membrane-confined electrophoresis-determined charges were previously reported. SILCS-Biologics maps the 3D distribution and predicted occupancy of ions, buffer molecules, and excipient molecules bound to the protein surface in a given salt environment. Using this information, the effective charge of the protein is predicted such that the concentrations of the ions and the presence of excipients or buffers are accounted for. Additionally, SILCS-Biologics also produces 3D structures of the binding sites of ions on the proteins, which enable further analyses such as the characterization of protein surface charge distribution and dipole moments in different environments. Notable is the capability of the method to account for competition between salts, excipients, and buffers on the calculated electrostatic properties in different protein formulations. Our study demonstrates the ability of the SILCS-Biologics approach to predict the effective charge of proteins and its applicability in uncovering protein-ion interactions and their contributions to protein solubility and function.

In silico identification of a β2-adrenoceptor allosteric site that selectively augments canonical β2AR-Gs signaling and function

Activation of β₂-adrenoceptors (β₂ARs) causes airway smooth muscle (ASM) relaxation and bronchodilation, and β₂AR agonists (β-agonists) are front-line treatments for asthma and other obstructive lung diseases. However, the therapeutic efficacy of β-agonists is limited by agonist-induced β₂AR desensitization and noncanonical β₂AR signaling involving β-arrestin that is shown to promote asthma pathophysiology. Accordingly, we undertook the identification of an allosteric site on β₂AR that could modulate the activity of β-agonists to overcome these limitations. We employed the site identification by ligand competitive saturation (SILCS) computational method to comprehensively map the entire 3D structure of in silico-generated β₂AR intermediate conformations and identified a putative allosteric binding site. Subsequent database screening using SILCS identified drug-like molecules with the potential to bind to the site. Experimental assays in HEK293 cells (expressing recombinant wild-type human β₂AR) and human ASM cells (expressing endogenous β₂AR) identified positive and negative allosteric modulators (PAMs and NAMs) of β₂AR as assessed by regulation of β-agonist-stimulation of cyclic AMP generation. PAMs/NAMs had no effect on β-agonist-induced recruitment of β-arrestin to β₂AR- or β-agonist-induced loss of cell surface expression in HEK293 cells expressing β₂AR. Mutagenesis analysis of β₂AR confirmed the SILCS identified site based on mutants of amino acids R131, Y219, and F282. Finally, functional studies revealed augmentation of β-agonist-induced relaxation of contracted human ASM cells and bronchodilation of contracted airways. These findings identify a allosteric binding site on the β₂AR, whose activation selectively augments β-agonist-induced Gs signaling, and increases relaxation of ASM cells, the principal therapeutic effect of β-agonists.

Effect of an Eco-Friendly Cuminaldehyde Guanylhydrazone Disinfectant on Shiga Toxin Production and Global Transcription of Escherichia coli

Antimicrobials have been important medicines used to treat various infections. However, some antibiotics increase the expression of Shiga toxin (Stx). Also, the pervasive use of persistent antibiotics has led to ecotoxicity and antibiotic resistance. In this study, a newly developed broad-spectrum and reversible antibiotic (guanylhydrazone disinfectant) was evaluated for its antibiotic activity and effects on Stx production and global transcription of bacteria. No Stx induction was observed in 25 Shiga toxin-producing E. coli (STEC) isolates treated with a sublethal concentration of the guanylhydrazone. A differential gene expression study comparing two guanylhydrazone-treated to non-treated E. coli strains indicated that the expression of a group of stress-responsive genes were enhanced. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that guanylhydrazone treatment significantly downregulated the pathways of ribosome and flagellar assembly in both pathogenic and non-pathogenic strains and differentially regulated some pathways essential for bacteria to maintain cell shape and gain survival advantage in two strains. In addition, upregulation of antibiotic resistant genes related to the multidrug efflux system and virulence genes coding for colibactin, colicin, and adhesin was observed in strains treated with the disinfectant. The knowledge obtained in this study contributes to our understanding of the mode of this disinfectant action and facilitates our effort to better use disinfectants for STEC treatments.

Recombinant Production of Biliverdin IXβ and δ Isomers in the T7 Promoter Compatible Escherichia coli Nissle

The ability to obtain purified biliverdin IX (BVIX) isomers other than the commercially available BVIXα is limited due to the low yields obtained by the chemical coupled oxidation of heme. Chemical oxidation requires toxic chemicals, has very poor BVIX yields (<0.05%), and is not conducive to scalable production. Alternative approaches utilizing recombinant E. coli BL21 expressing a cyanobacterial heme oxygenase have been employed for the production BVIXα, but yields are limited by the rate of endogenous heme biosynthesis. Furthermore, the emerging roles of BVIXβ and BVIXδ in biology and their lack of commercial availability has led to a need for an efficient and scalable method with the flexibility to produce all three physiologically relevant BVIX isomers. Herein, we have taken advantage of an optimized non-pathogenic E. coli Nissle (EcN(T7)) strain that encodes an endogenous heme transporter and an integrated T7 polymerase gene. Protein production of the Pseudomonas aeruginosa BVIXβ and BVIXδ selective heme oxygenase (HemO) or its BVIXα producing mutant (HemOα) in the EcN(T7) strain provides a scalable method to obtain all three isomers, that is not limited by the rate of endogenous heme biosynthesis, due to the natural ability of EcN(T7) to transport extracellular heme. Additionally, we have optimized our previous LC-MS/MS protocol for semi-preparative separation and validation of the BVIX isomers. Utilizing this new methodology for scalable production and separation we have increased the yields of the BVIXβ and -δ isomers >300-fold when compared to the chemical oxidation of heme.

Extracellular haem utilization by the opportunistic pathogen Pseudomonas aeruginosa and its role in virulence and pathogenesis

Iron is an essential micronutrient for all bacteria but presents a significant challenge given its limited bioavailability. Furthermore, iron's toxicity combined with the need to maintain iron levels within a narrow physiological range requires integrated systems to sense, regulate and transport a variety of iron complexes. Most bacteria encode systems to chelate and transport ferric iron (Fe³⁺) via siderophore receptor mediated uptake or via cytoplasmic energy dependent transport systems. Pathogenic bacteria have further lowered the barrier to iron acquisition by employing systems to utilize haem as a source of iron. Haem, a lipophilic and toxic molecule, presents a significant challenge for transport into the cell. As such pathogenic bacteria have evolved sophisticated cell surface signaling (CSS) and transport systems to sense and obtain haem from the host. Once internalized haem is cleaved by both oxidative and non-oxidative mechanisms to release iron. Herein we summarize our current understanding of the mechanism of haem sensing, uptake and utilization in Pseudomonas aeruginosa, its role in pathogenesis and virulence, and the potential of these systems as antimicrobial targets.

Rapid and accurate estimation of protein-ligand relative binding affinities using site-identification by ligand competitive saturation

Predicting relative protein-ligand binding affinities is a central pillar of lead optimization efforts in structure-based drug design. The site identification by ligand competitive saturation (SILCS) methodology is based on functional group affinity patterns in the form of free energy maps that may be used to compute protein-ligand binding poses and affinities. Presented are results obtained from the SILCS methodology for a set of eight target proteins as reported originally in Wang et al. (J. Am. Chem. Soc., 2015, 137, 2695-2703) using free energy perturbation (FEP) methods in conjunction with enhanced sampling and cycle closure corrections. These eight targets have been subsequently studied by many other authors to compare the efficacy of their method while comparing with the outcomes of Wang et al. In this work, we present results for a total of 407 ligands on the eight targets and include specific analysis on the subset of 199 ligands considered previously. Using the SILCS methodology we can achieve an average accuracy of up to 77% and 74% when considering the eight targets with their 199 and 407 ligands, respectively, for rank-ordering ligand affinities as calculated by the percent correct metric. This accuracy increases to 82% and 80%, respectively, when the SILCS atomic free energy contributions are optimized using a Bayesian Markov-chain Monte Carlo approach. We also report other metrics including Pearson's correlation coefficient, Pearlman's predictive index, mean unsigned error, and root mean square error for both sets of ligands. The results obtained for the 199 ligands are compared with the outcomes of Wang et al. and other published works. Overall, the SILCS methodology yields similar or better-quality predictions without a priori need for known ligand orientations in terms of the different metrics when compared to current FEP approaches with significant computational savings while additionally offering quantitative estimates of individual atomic contributions to binding free energies. These results further validate the SILCS methodology as an accurate, computationally efficient tool to support lead optimization and drug discovery.

Repurposing Acitretin as an Antipseudomonal Agent Targeting the Pseudomonas aeruginosa Iron-Regulated Heme Oxygenase

Iron is an essential micronutrient for the survival and virulence of the bacterial pathogen Pseudomonas aeruginosa. To overcome iron withholding and successfully colonize a host, P. aeruginosa uses a variety of mechanisms to acquire iron, including the secretion of high-affinity iron chelators (siderophores) or the uptake and utilization of heme. P. aeruginosa heme oxygenase (HemO) plays pivotal roles in heme sensing, uptake, and utilization and has emerged as a therapeutic target for the development of antipseudomonal agents. Using a high-throughput fluorescence quenching assay combined with minimum inhibitory concentration measurements, we screened the Selleck Bioactive collection of 2100 compounds and identified acitretin, a Food and Drug Administration-approved oral retinoid, as a potent and selective inhibitor of HemO. Acitretin binds to HemO with a KD value of 0.10 ± 0.02 μM and inhibits the growth of P. aeruginosa PAO1 with an IC50 of 70 ± 18 μg/mL. In addition, acitretin showed good selectivity for HemO, which uniquely generates BVIXβ/δ, over human heme oxygenase (hHO1) and other BVIXα-producing homologues such as the heme oxygenases from Neisseria meningitidis (nmHO) and Acinetobacter baumannii (abHO). The binding of acitretin within the HemO active site was confirmed by 1H-15N heteronuclear single-quantum coherence nuclear magnetic resonance, and molecular modeling provided further insight into potential interactions of acitretin with residues specific for orienting heme in the β/δ selective HemO. Moreover, at 20 μM, acitretin inhibited the enzymatic activity of HemO in P. aeruginosa cells by >60% and effectively blocked the ability of P. aeruginosa to sense and acquire heme as demonstrated in the β-galactosidase transcriptional reporter assay.

Small Molecule Inhibitors of the Bacterioferritin (BfrB)-Ferredoxin (Bfd) Complex Kill Biofilm-Embedded Pseudomonas aeruginosa Cells

Bacteria depend on a well-regulated iron homeostasis to survive adverse environments. A key component of the iron homeostasis machinery is the compartmentalization of Fe3+ in bacterioferritin and its subsequent mobilization as Fe2+ to satisfy metabolic requirements. In Pseudomonas aeruginosa Fe3+ is compartmentalized in bacterioferritin (BfrB), and its mobilization to the cytosol requires binding of a ferredoxin (Bfd) to reduce the stored Fe3+ and release the soluble Fe2+. Blocking the BfrB-Bfd complex in P. aeruginosa by deletion of the bfd gene triggers an irreversible accumulation of Fe3+ in BfrB, concomitant cytosolic iron deficiency and significant impairment of biofilm development. Herein we report that small molecules developed to bind BfrB at the Bfd binding site block the BfrB-Bfd complex, inhibit the mobilization of iron from BfrB in P. aeruginosa cells, elicit a bacteriostatic effect on planktonic cells, and are bactericidal to cells embedded in mature biofilms.

Immune Modulation by Inhibitors of the HO System

The heme oxygenase (HO) system involves three isoforms of this enzyme, HO-1, HO-2, and HO-3. The three of them display the same catalytic activity, oxidating the heme group to produce biliverdin, ferrous iron, and carbon monoxide (CO). HO-1 is the isoform most widely studied in proinflammatory diseases because treatments that overexpress this enzyme promote the generation of anti-inflammatory products. However, neonatal jaundice (hyperbilirubinemia) derived from HO overexpression led to the development of inhibitors, such as those based on metaloproto- and meso-porphyrins inhibitors with competitive activity. Further, non-competitive inhibitors have also been identified, such as synthetic and natural imidazole-dioxolane-based, small synthetic molecules, inhibitors of the enzyme regulation pathway, and genetic engineering using iRNA or CRISPR cas9. Despite most of the applications of the HO inhibitors being related to metabolic diseases, the beneficial effects of these molecules in immune-mediated diseases have also emerged. Different medical implications, including cancer, Alzheimer´s disease, and infections, are discussed in this article and as to how the selective inhibition of HO isoforms may contribute to the treatment of these ailments.

Assessing hERG1 Blockade from Bayesian Machine-Learning-Optimized Site Identification by Ligand Competitive Saturation Simulations

Drug-induced cardiotoxicity is a potentially lethal and yet one of the most common side effects with the drugs in clinical use. Most of the drug-induced cardiotoxicity is associated with an off-target pharmacological blockade of K+ currents carried out by the cardiac Human-Ether-a-go-go-Related (hERG1) potassium channel. There is a compulsory preclinical stage safety assessment for the hERG1 blockade for all classes of drugs, which adds substantially to the cost of drug development. The availability of a high-resolution cryogenic electron microscopy (cryo-EM) structure for the channel in its open/depolarized state solved in 2017 enabled the application of molecular modeling for rapid assessment of drug blockade by molecular docking and simulation techniques. More importantly, if successful, in silico methods may allow a path to lead-compound salvaging by mapping out key block determinants. Here, we report the blind application of the site identification by the ligand competitive saturation (SILCS) protocol to map out druggable/regulatory hotspots in the hERG1 channel available for blockers and activators. The SILCS simulations use small solutes representative of common functional groups to sample the chemical space for the entire protein and its environment using all-atom simulations. The resulting chemical maps, FragMaps, explicitly account for receptor flexibility, protein-fragment interactions, and fragment desolvation penalty allowing for rapid ranking of potential ligands as blockers or nonblockers of hERG1. To illustrate the power of the approach, SILCS was applied to a test set of 55 blockers with diverse chemical scaffolds and pIC50 values measured under uniform conditions. The original SILCS model was based on the all-atom modeling of the hERG1 channel in an explicit lipid bilayer and was further augmented with a Bayesian-optimization/machine-learning (BML) stage employing an independent literature-derived training set of 163 molecules. BML approach was used to determine weighting factors for the FragMaps contributions to the scoring function. pIC50 predictions from the combined SILCS/BML approach to the 55 blockers showed a Pearson correlation (PC) coefficient of >0.535 relative to the experimental data. SILCS/BML model was shown to yield substantially improved performance as compared to commonly used rigid and flexible molecular docking methods for a well-established cohort of hERG1 blockers, where no correlation with experimental data was recorded. SILCS/BML results also suggest that a proper weighting of protonation states of common blockers present at physiological pH is essential for accurate predictions of blocker potency. The precalculated and optimized SILCS FragMaps can now be used for the rapid screening of small molecules for their cardiotoxic potential as well as for exploring alternative binding pockets in the hERG1 channel with applications to the rational design of activators.

Metallotherapeutics development in the age of iron-clad bacteria

Drug-resistant infections pose a significant risk to global health as pathogenic bacteria become increasingly difficult to treat. The rapid selection of resistant strains through poor antibiotic stewardship has reduced the number of viable treatments and increased morbidity of infections, especially among the immunocompromised. To circumvent such challenges, new strategies are required to stay ahead of emerging resistance trends, yet research and funding for antibiotic development lags other classes of therapeutics. Though the use of metals in therapeutics has been around for centuries, recent strategies have devoted a great deal of effort into the pathways through which bacteria acquire and utilize iron, which is critical for the establishment of infection. To target iron uptake systems, siderophore-drug conjugates have been developed that hijack siderophore-based iron uptake for delivery of antibiotics. While this strategy has produced several potential leads, the use of siderophores in infection is diminished over time when bacteria adapt to utilize heme as an iron source, leading to a need for the development of porphyrin mimetics as therapeutics. The use of such strategies as well as the inclusion of gallium, a redox-inert iron mimic, are herein reviewed.

Mobilization of Iron Stored in Bacterioferritin Is Required for Metabolic Homeostasis in Pseudomonas aeruginosa

Iron homeostasis offers a significant bacterial vulnerability because pathogens obtain essential iron from their mammalian hosts, but host-defenses maintain vanishingly low levels of free iron. Although pathogens have evolved mechanisms to procure host-iron, these depend on well-regulated iron homeostasis. To disrupt iron homeostasis, our work has targeted iron mobilization from the iron storage protein bacterioferritin (BfrB) by blocking a required interaction with its cognate ferredoxin partner (Bfd). The blockade of the BfrB–Bfd complex by deletion of the bfd gene (Δbfd) causes iron to irreversibly accumulate in BfrB. In this study we used mass spectrometry and NMR spectroscopy to compare the proteomic response and the levels of key intracellular metabolites between wild type (wt) and isogenic ΔbfdP. aeruginosa strains. We find that the irreversible accumulation of unusable iron in BfrB leads to acute intracellular iron limitation, even if the culture media is iron-sufficient. Importantly, the iron limitation and concomitant iron metabolism dysregulation trigger a cascade of events that lead to broader metabolic homeostasis disruption, which includes sulfur limitation, phenazine-mediated oxidative stress, suboptimal amino acid synthesis and altered carbon metabolism.

Identification and characterization of fragment binding sites for allosteric ligand design using the site identification by ligand competitive saturation hotspots approach (SILCS-Hotspots)

BACKGROUND

Fragment-based ligand design is used for the development of novel ligands that target macromolecules, most notably proteins. Central to its success is the identification of fragment binding sites that are spatially adjacent such that fragments occupying those sites may be linked to create drug-like ligands. Current experimental and computational approaches that address this problem typically identify only a limited number of sites as well as use a limited number of fragment types.

METHODS

The site-identification by ligand competitive saturation (SILCS) approach is extended to the identification of fragment bindings sites, with the method termed SILCS-Hotspots. The approach involves precomputation of the SILCS FragMaps following which the identification of Hotspots, performed by identifying of all possible fragment binding sites on the full 3D structure of the protein followed by spatial clustering.

RESULTS

The SILCS-Hotspots approach identifies a large number of sites on the target protein, including many sites not accessible in experimental structures due to low binding affinities and binding sites on the protein interior. The identified sites are shown to recapitulate the location of known drug-like molecules in both allosteric and orthosteric binding sites on seven proteins including the androgen receptor, the CDK2 and Erk5 kinases, PTP1B phosphatase and three GPCRs; the β2-adrenergic, GPR40 fatty-acid binding and M2-muscarinic receptors. Analysis indicates the importance of considering all possible fragment binding sites, and not just those accessible to experimental methods, when identifying novel binding sites and performing ligand design versus just considering the most favorable sites. The approach is shown to identify a larger number of known binding sites of drug-like molecules versus the commonly used FTMap and Fpocket methods.

GENERAL SIGNIFICANCE

The present results indicate the potential utility of the SILCS-Hotspots approach for fragment-based rational design of ligands, including allosteric modulators.

Aromatic Guanylhydrazones for the Control of Heme-Induced Antibody Polyreactivity

In a healthy immune repertoire, there exists a fraction of polyreactive antibodies that can bind to a variety of unrelated self- and foreign antigens. Apart from naturally polyreactive antibodies, in every healthy individual, there is a fraction of antibody that can gain polyreactivity upon exposure to porphyrin cofactor heme. Molecular mechanisms and biological significance of the appearance of cryptic polyreactivity are not well understood. It is believed that heme acts as an interfacial cofactor between the antibody and the newly recognized antigens. To further test this claim and gain insight into the types of interactions involved in heme binding, we herein investigated the influence of a group of aromatic guanylhydrazone molecules on the heme-induced antibody polyreactivity. From the analysis of SAR and the results of UV-vis absorbance spectroscopy, it was concluded that the most probable mechanism by which the studied molecules inhibit heme-mediated polyreactivity of the antibody is the direct binding to heme, thus preventing heme from binding to antibody and/or antigen. The inhibitory capacity of the most potent compounds was substantially higher than that of chloroquine, a well-known heme binder. Some of the guanylhydrazone molecules were able to induce polyreactivity of the studied antibody themselves, possibly by a mechanism similar to heme. Results described here point to the conclusion that heme indeed must bind to an antibody to induce its polyreactivity, and that both π-stacking interactions and iron coordination contribute to the binding affinity, while certain structures, such as guanylhydrazones, can interfere with these processes.

Incorporating G-C Pair-Recognizing Guanidinium into PNAs for Sequence and Structure Specific Recognition of dsRNAs over dsDNAs and ssRNAs

Recognition of RNAs under physiological conditions is important for the development of chemical probes and therapeutic ligands. Nucleobase-modified dsRNA-binding PNAs (dbPNAs) are promising for the recognition of dsRNAs in a sequence and structure specific manner under near-physiological conditions. Guanidinium is often present in proteins and small molecules for the recognition of G bases in nucleic acids, in cell-penetrating carriers, and in bioactive drug molecules, which might be due to the fact that guanidinium is amphiphilic and has unique hydrogen bonding and stacking properties. We hypothesized that a simple guanidinium moiety can be directly incorporated into PNAs to facilitate enhanced molecular recognition of G-C pairs in dsRNAs and improved bioactivity. We grafted a guanidinium moiety directly into a PNA monomer (designated as R) using a two-carbon linker as guided by computational modeling studies. The synthetic scheme of the PNA R monomer is relatively simple compared to that of the previously reported L monomer. We incorporated the R residue into various dbPNAs for binding studies. dbPNAs incorporated with R residues are excellent in sequence specifically recognizing G-C pairs in dsRNAs over dsDNA and ssRNAs. We demonstrated that the R residue is compatible with unmodified T and C and previously developed modified L and Q residues in dbPNAs for targeting model dsRNAs, the influenza A viral panhandle duplex structure, and the HIV-1 frameshift site RNA hairpin. Furthermore, R residues enhance the cellular uptake of PNAs.

Heme uptake and utilization by hypervirulent Acinetobacter baumannii LAC-4 is dependent on a canonical heme oxygenase (abHemO)

Acinetobacter baumannii is an opportunistic pathogen that causes serious infections in critically ill and immune compromised patients. The ability to acquire iron from the hosts iron and heme containing proteins is critical to their survival and virulence. The majority of A. baumannii hypervirulent strains encode a heme uptake system that includes a putative heme oxygenase (hemO). Despite reports indicating A. baumannii can grow on heme direct evidence of extracellular heme uptake and metabolism has not been shown. Through isotopic labeling (¹³C-heme) we show the hypervirulent A. baumannii LAC-4 metabolizes heme to biliverdin IXα (BVIXα), whereas ATC 17978 that lacks the hemO gene cluster cannot efficiently utilize heme. Expression and purification of the protein encoded by the A. baumannii LAC-4 hemO gene confirmed catalytic conversion of heme to BVIX. We further show inhibition of abHemO with previously characterized P. aeruginosa HemO inhibitors in a fluorescence based assay that couples HemO catalytic activity to the BVIXα binding phytochrome IFP1.4. Furthermore, the hemO gene cluster encodes genes with homology to heme-dependent extra cytoplasmic function (ECF) σ factor systems. The hemophore-dependent ECF system in Pseudomonas aeruginosa has been shown to play a critical role in heme sensing and virulence within the host. The prevalence of a hemO gene cluster in A. baumannii LAC4 and other hypervirulent strains suggests it is required within the host to adapt and utilize heme and is a major contributor to virulence.

Optimization and Evaluation of Site-Identification by Ligand Competitive Saturation (SILCS) as a Tool for Target-Based Ligand Optimization

Chemical fragment cosolvent sampling techniques have become a versatile tool in ligand-protein binding prediction. Site-identification by ligand competitive saturation (SILCS) is one such method that maps the distribution of chemical fragments on a protein as free energy fields called FragMaps. Ligands are then simulated via Monte Carlo techniques in the field of the FragMaps (SILCS-MC) to predict their binding conformations and relative affinities for the target protein. Application of SILCS-MC using a number of different scoring schemes and MC sampling protocols against multiple protein targets was undertaken to evaluate and optimize the predictive capability of the method. Seven protein targets and 551 ligands with broad chemical variability were used to evaluate and optimize the model to maximize Pearson's correlation coefficient, Pearlman's predictive index, correct relative binding affinity, and root-mean-square error versus the absolute experimental binding affinities. Across the protein-ligand sets, the relative affinities of the ligands were predicted correctly an average of 69% of the time for the highest overall SILCS protocol. Training the FragMap weighting factors using a Bayesian machine learning (ML) algorithm led to an increase to an average 75% relative correct affinity predictions. Furthermore, once the optimal protocol is identified for a specific protein-ligand system average predictabilities of 76% are achieved. The ML algorithm is successful with small training sets of data (30 or more compounds) due to the use of physically correct FragMap weights as priors. Notably, the 76% correct relative prediction rate is similar to or better than free energy perturbation methods that are significantly computationally more expensive than SILCS. The results further support the utility of SILCS as a powerful and computationally accessible tool to support lead optimization and development in drug discovery.

Emerging Opportunities To Manipulate Metal Trafficking for Therapeutic Benefit

The indispensable requirement for metals in life processes has led to the evolution of sophisticated mechanisms that allow organisms to maintain dynamic equilibria of these ions. This dynamic control of the level, speciation, and availability of a variety of metal ions allows organisms to sustain biological processes while avoiding toxicity. When functioning properly, these mechanisms allow cells to return to their metal homeostatic set points following shifts in the metal availability or other stressors. These periods of transition, when cells are in a state of flux in which they work to regain homeostasis, present windows of opportunity to pharmacologically manipulate targets associated with metal-trafficking pathways in ways that could either facilitate a return to homeostasis and the recovery of cellular function or further push cells outside of homeostasis and into cellular distress. The purpose of this Viewpoint is to highlight emerging opportunities for chemists and chemical biologists to develop compounds to manipulate metal-trafficking processes for therapeutic benefit.

Small Molecule Inhibitors of the BfrB-Bfd Interaction Decrease Pseudomonas aeruginosa Fitness and Potentiate Fluoroquinolone Activity

The iron storage protein bacterioferritin (BfrB) is central to bacterial iron homeostasis. The mobilization of iron from BfrB, which requires binding by a cognate ferredoxin (Bfd), is essential to the regulation of cytosolic iron levels in P. aeruginosa. This paper describes the structure-guided development of small molecule inhibitors of the BfrB–Bfd protein–protein interaction. The process was initiated by screening a fragment library and followed by obtaining the structure of a fragment hit bound to BfrB. The structural insights were used to develop a series of 4-(benzylamino)- and 4-((3-phenylpropyl)amino)-isoindoline-1,3-dione analogs that selectively bind BfrB at the Bfd binding site. Challenging P. aeruginosa cells with the 4-substituted isoindoline analogs revealed a dose-dependent growth phenotype. Further investigation determined that the analogs elicit a pyoverdin hyperproduction phenotype that is consistent with blockade of the BfrB–Bfd interaction and ensuing irreversible accumulation of iron in BfrB, with concomitant depletion of iron in the cytosol. The irreversible accumulation of iron in BfrB prompted by the 4-substituted isoindoline analogs was confirmed by visualization of BfrB-iron in P. aeruginosa cell lysates separated on native PAGE gels and stained for iron with Ferene S. Challenging P. aeruginosa cultures with a combination of commercial fluoroquinolone and our isoindoline analogs results in significantly lower cell survival relative to treatment with either antibiotic or analog alone. Collectively, these findings furnish proof of concept for the usefulness of small molecule probes designed to dysregulate bacterial iron homeostasis by targeting a protein–protein interaction pivotal for iron storage in the bacterial cell.

Exploring protein-protein interactions using the site-identification by ligand competitive saturation methodology

Protein docking methods are powerful computational tools to study protein-protein interactions (PPI). While a significant number of docking algorithms have been developed, they are usually based on rigid protein models or with limited considerations of protein flexibility and the desolvation effect is rarely considered in docking energy functions, which may lower the accuracy of the predictions. To address these issues, we introduce a PPI energy function based on the site-identification by ligand competitive saturation (SILCS) framework and utilize the fast Fourier transform (FFT) correlation approach. The free energy content of the SILCS FragMaps represent an alternative to traditional energy grids and they can be efficiently utilized to guide FFT-based protein docking. Application of the approach to eight diverse test cases, including seven from Protein Docking Benchmark 5.0, showed the PPI prediction using SILCS approach (SILCS-PPI) to be competitive with several commonly used protein docking methods indicating that the method has the ability to both qualitatively and quantitatively inform the prediction of PPI. Results show the utility of the SILCS-PPI docking approach for determination of probability distributions of PPI interactions over the surface of both partner proteins, allowing for identification of alternate binding poses. Such binding poses are confirmed by experimental crystal contacts in our test cases. While more computationally demanding than available PPI docking technologies, we anticipate that the SILCS-PPI docking approach will offer an alternative methodology for improved evaluation of PPIs that could be used in a variety of fields from systems biology to excipient design for biologics-based drugs.

Azines: synthesis, structure, electronic structure and their applications

Azines (2,3-diaza-1,3-butadienes): structure, electronic structure, tautomerism, and their applications in organic synthesis, medicinal chemistry and materials chemistry.

Identification of Thiourea-Based Inhibitors of the B-Cell Lymphoma 6 BTB Domain via NMR-Based Fragment Screening and Computer-Aided Drug Design

Protein-protein interactions (PPI) between the transcriptional repressor B-cell lymphoma 6 (BCL6) BTB domain (BCL6^BTB) and its corepressors have emerged as a promising target for anticancer therapeutics. However, identification of potent, drug-like inhibitors of BCL6^BTB has remained challenging. Using NMR-based screening of a library of fragment-like small molecules, we have identified a thiourea compound (7CC5) that binds to BCL6^BTB. From this hit, the application of computer-aided drug design (CADD), medicinal chemistry, NMR spectroscopy, and X-ray crystallography has yielded an inhibitor, 15f, that demonstrated over 100-fold improved potency for BCL6^BTB. This gain in potency was achieved by a unique binding mode that mimics the binding mode of the corepressor SMRT in the aromatic and the HDCH sites. The structure-activity relationship based on these new inhibitors will have a significant impact on the rational design of novel BCL6 inhibitors, facilitating the identification of therapeutics for the treatment of BCL6-dependent tumors.

The Asp99-Arg188 salt bridge of the Pseudomonas aeruginosa HemO is critical in allowing conformational flexibility during catalysis

The P. aeruginosa iron-regulated heme oxygenase (HemO) is required within the host for the utilization of heme as an iron source. As iron is essential for survival and virulence, HemO represents a novel antimicrobial target. We recently characterized small molecule inhibitors that bind to an allosteric site distant from the heme pocket, and further proposed binding at this site disrupts a nearby salt bridge between D99 and R188. Herein, through a combination of site-directed mutagenesis and hydrogen-deuterium exchange mass spectrometry (HDX-MS), we determined that the disruption of the D99-R188 salt bridge leads to significant decrease in conformational flexibility within the distal and proximal helices that form the heme-binding site. The RR spectra of the resting state Fe(III) and reduced Fe(II)-deoxy heme-HemO D99A, R188A and D99/R188A complexes are virtually identical to those of wild-type HemO, indicating no significant change in the heme environment. Furthermore, mutation of D99 or R188 leads to a modest decrease in the stability of the Fe(II)-O₂ heme complex. Despite this slight difference in Fe(II)-O₂ stability, we observe complete loss of enzymatic activity. We conclude the loss of activity is a result of decreased conformational flexibility in helices previously shown to be critical in accommodating variation in the distal ligand and the resulting chemical intermediates generated during catalysis. Furthermore, this newly identified allosteric binding site on HemO represents a novel alternative drug-design strategy to that of competitive inhibition at the active site or via direct coordination of ligands to the heme iron.

Tight binding of heme to Staphylococcus aureus IsdG and IsdI precludes design of a competitive inhibitor

The micromolar equilibrium constants for heme dissociation from IsdG and IsdI reported in the literature call into question whether these enzymes are actually members of the iron-regulated surface determinant system of Staphylococcus aureus, which harvests heme iron from a host during infection. In order to address this question, the heme dissociation constants for IsdG and IsdI were reevaluated using three approaches. The heme dissociation equilibrium constants were measured using a UV/Vis absorption-detected assay analyzed with an assumption-free model, and using a newly developed fluorescence-detected assay. The heme dissociation rate constants were estimated using apomyoglobin competition assays. Analyses of the UV/Vis absorption data revealed a critical flaw in the previous measurements; heme is 99.9% protein-bound at the micromolar concentrations needed for UV/Vis absorption spectroscopy, which renders accurate equilibrium constant measurement nearly impossible. However, fluorescence can be measured for more dilute samples, and analyses of these data resulted in dissociation equilibrium constants of 1.4 ± 0.6 nM and 12.9 ± 1.3 nM for IsdG and IsdI, respectively. Analyses of the kinetic data obtained from apomyoglobin competition assays estimated heme dissociation rate constants of 0.022 ± 0.002 s^-1 for IsdG and 0.092 ± 0.008 s^-1 for IsdI. Based upon these data, and what is known regarding the post-translational regulation of IsdG and IsdI, it is proposed that only IsdG is a member of the heme iron acquisition pathway and IsdI regulates heme homeostasis. Furthermore, the nanomolar dissociation constants mean that heme is bound tightly by IsdG and indicates that competitive inhibition of this protein will be difficult. Instead, uncompetitive inhibition based upon a detailed understanding of enzyme mechanism is a more promising antibiotic development strategy.

Structural and mutational analyses of the Leptospira interrogans virulence-related heme oxygenase provide insights into its catalytic mechanism

Heme oxygenase from Leptospira interrogans is an important virulence factor. During catalysis, redox equivalents are provided to this enzyme by the plastidic-type ferredoxin-NADP⁺ reductase also found in L. interrogans. This process may have evolved to aid this bacterial pathogen to obtain heme-iron from their host and enable successful colonization. Herein we report the crystal structure of the heme oxygenase-heme complex at 1.73 Å resolution. The structure reveals several distinctive features related to its function. A hydrogen bonded network of structural water molecules that extends from the catalytic site to the protein surface was cleared observed. A depression on the surface appears to be the H⁺ network entrance from the aqueous environment to the catalytic site for O₂ activation, a key step in the heme oxygenase reaction. We have performed a mutational analysis of the F157, located at the above-mentioned depression. The mutant enzymes were unable to carry out the complete degradation of heme to biliverdin since the reaction was arrested at the verdoheme stage. We also observed that the stability of the oxyferrous complex, the efficiency of heme hydroxylation and the subsequent conversion to verdoheme was adversely affected. These findings underscore a long-range communication between the outer fringes of the hydrogen-bonded network of structural waters and the heme active site during catalysis. Finally, by analyzing the crystal structures of ferredoxin-NADP⁺ reductase and heme oxygenase, we propose a model for the productive association of these proteins.

Benzene Probes in Molecular Dynamics Simulations Reveal Novel Binding Sites for Ligand Design

Protein flexibility poses a major challenge in binding site identification. Several computational pocket detection methods that utilize small-molecule probes in molecular dynamics (MD) simulations have been developed to address this issue. Although they have proven hugely successful at reproducing experimental structural data, their ability to predict new binding sites that are yet to be identified and characterized has not been demonstrated. Here, we report the use of benzenes as probe molecules in ligand-mapping MD (LMMD) simulations to predict the existence of two novel binding sites on the surface of the oncoprotein MDM2. One of them was serendipitously confirmed by biophysical assays and X-ray crystallography to be important for the binding of a new family of hydrocarbon stapled peptides that were specifically designed to target the other putative site. These results highlight the predictive power of LMMD and suggest that predictions derived from LMMD simulations can serve as a reliable basis for the identification of novel ligand binding sites in structure-based drug design.