Nicolas Foloppe

Novel adenosine-derived inhibitors of 70 kDa heat shock protein, discovered through structure-based design

The design and synthesis of novel adenosine-derived inhibitors of HSP70, guided by modeling and X-ray crystallographic structures of these compounds in complex with HSC70/BAG-1, is described. Examples exhibited submicromolar affinity for HSP70, were highly selective over HSP90, and some displayed potency against HCT116 cells. Exposure of compound 12 to HCT116 cells caused significant reduction in cellular levels of Raf-1 and Her2 at concentrations similar to that which caused cell growth arrest.

Inhibition of the Heat Shock Protein 90 Molecular Chaperone in Vitro and in Vivo by Novel, Synthetic, Potent Resorcinylic Pyrazole/Isoxazole Amide Analogues.

Although the heat shock protein 90 (HSP90) inhibitor 17-allylamino-17-demethoxygeldanamycin (17-AAG) shows clinical promise, potential limitations encourage development of alternative chemotypes. We discovered the 3,4-diarylpyrazole resorcinol CCT018159 by high-throughput screening and used structure-based design to generate more potent pyrazole amide analogues, exemplified by VER-49009. Here, we describe the detailed biological properties of VER-49009 and the corresponding isoxazole VER-50589. X-ray crystallography showed a virtually identical HSP90 binding mode. However, the dissociation constant (Kd) of VER-50589 was 4.5 ± 2.2 nmol/L compared with 78.0 ± 10.4 nmol/L for VER-49009, attributable to higher enthalpy for VER-50589 binding. A competitive binding assay gave a lower IC50 of 21 ± 4 nmol/L for VER-50589 compared with 47 ± 9 nmol/L for VER-49009. Cellular uptake of VER-50589 was 4-fold greater than for VER-49009. Mean cellular antiproliferative GI50 values for VER-50589 and VER-49009 for a human cancer cell line panel were 78 ± 15 and 685 ± 119 nmol/L, respectively, showing a 9-fold potency gain for the isoxazole. Unlike 17-AAG, but as with CCT018159, cellular potency of these analogues was independent of NAD(P)H:quinone oxidoreductase 1/DT-diaphorase and P-glycoprotein expression. Consistent with HSP90 inhibition, VER-50589 and VER-49009 caused induction of HSP72 and HSP27 alongside depletion of client proteins, including C-RAF, B-RAF, and survivin, and the protein arginine methyltransferase PRMT5. Both caused cell cycle arrest and apoptosis. Extent and duration of pharmacodynamic changes in an orthotopic human ovarian carcinoma model confirmed the superiority of VER-50589 over VER-49009. VER-50589 accumulated in HCT116 human colon cancer xenografts at levels above the cellular GI50 for 24 h, resulting in 30% growth inhibition. The results indicate the therapeutic potential of the resorcinylic pyrazole/isoxazole amide analogues as HSP90 inhibitors. [Mol Cancer Ther 2007;6(4):1198–211]

Understanding the pKa of Redox Cysteines: The Key Role of Hydrogen Bonding

Many cellular functions involve cysteine chemistry via thiol-disulfide exchange pathways. The nucleophilic cysteines of the enzymes involved are activated as thiolate. A thiolate is much more reactive than a neutral thiol. Therefore, determining and understanding the pK(a)s of functional cysteines are important aspects of biochemistry and molecular biology with direct implications for redox signaling. Here, we describe the experimental and theoretical methods to determine cysteine pK(a) values, and we examine the factors that control these pK(a)s. Drawing largely on experience gained with the thioredoxin superfamily, we examine the roles of solvation, charge-charge, helix macrodipole, and hydrogen bonding interactions as pK(a)-modulating factors. The contributions of these factors in influencing cysteine pK(a)s and the associated chemistry, including the relevance for the reaction kinetics and thermodynamics, are discussed. This analysis highlights the critical role of direct hydrogen bonding to the cysteine sulfur as a key factor modulating the equilibrium between thiol S-H and thiolate S(-). This role is easily understood intuitively and provides a framework for biochemical functional insights.

Drug-like annotation and duplicate analysis of a 23-supplier chemical database totalling 2.7 million compounds

We have implemented five drug-like filters, based on 1D and 2D molecular descriptors, and applied them to characterize the drug-like properties of commercially available chemical compounds. In addition to previously published filters (Lipinski and Veber), we implemented a filter for medicinal chemistry tractability based on lists of chemical features drawn up by a panel of medicinal chemists. A filter based on the modeling of aqueous solubility (>1 microM) was derived in-house, as well as another based on the modeling of Caco-2 passive membrane permeability (>10 nm/s). A library of 2.7 million compounds was collated from the 23 compound suppliers and analyzed with these filters, highlighting a tendency toward highly lipophilic compounds. The library contains 1.6 M unique structures, of which 37% (607,223) passed all five drug-like filters. None of the 23 suppliers provides all the members of the drug-like subset, emphasizing the benefit of considering compounds from various compound suppliers as a source of diversity for drug discovery.

rDock: A Fast, Versatile and Open Source Program for Docking Ligands to Proteins and Nucleic Acids

Identification of chemical compounds with specific biological activities is an important step in both chemical biology and drug discovery. When the structure of the intended target is available, one approach is to use molecular docking programs to assess the chemical complementarity of small molecules with the target; such calculations provide a qualitative measure of affinity that can be used in virtual screening (VS) to rank order a list of compounds according to their potential to be active. rDock is a molecular docking program developed at Vernalis for high-throughput VS (HTVS) applications. Evolved from RiboDock, the program can be used against proteins and nucleic acids, is designed to be computationally very efficient and allows the user to incorporate additional constraints and information as a bias to guide docking. This article provides an overview of the program structure and features and compares rDock to two reference programs, AutoDock Vina (open source) and Schrödingers Glide (commercial). In terms of computational speed for VS, rDock is faster than Vina and comparable to Glide. For binding mode prediction, rDock and Vina are superior to Glide. The VS performance of rDock is significantly better than Vina, but inferior to Glide for most systems unless pharmacophore constraints are used; in that case rDock and Glide are of equal performance. The program is released under the Lesser General Public License and is freely available for download, together with the manuals, example files and the complete test sets, at http://rdock.sourceforge.net/

Mycoredoxin-1 is one of the missing links in the oxidative stress defence mechanism of Mycobacteria.

To survive hostile conditions, the bacterial pathogen Mycobacterium tuberculosis produces millimolar concentrations of mycothiol as a redox buffer against oxidative stress. The reductases that couple the reducing power of mycothiol to redox active proteins in the cell are not known. We report a novel mycothiol‐dependent reductase (mycoredoxin‐1) with a CGYC catalytic motif. With mycoredoxin‐1 and mycothiol deletion strains of Mycobacterium smegmatis, we show that mycoredoxin‐1 and mycothiol are involved in the protection against oxidative stress. Mycoredoxin‐1 acts as an oxidoreductase exclusively linked to the mycothiol electron transfer pathway and it can reduce S‐mycothiolated mixed disulphides. Moreover, we solved the solution structures of oxidized and reduced mycoredoxin‐1, revealing a thioredoxin fold with a putative mycothiol‐binding site. With HSQC snapshots during electron transport, we visualize the reduction of oxidized mycoredoxin‐1 as a function of time and find that mycoredoxin‐1 gets S‐mycothiolated on its N‐terminal nucleophilic cysteine. Mycoredoxin‐1 has a redox potential of −218 mV and hydrogen bonding with neighbouring residues lowers the pKa of its N‐terminal nucleophilic cysteine. Determination of the oxidized and reduced structures of mycoredoxin‐1, better understanding of mycothiol‐dependent reactions in general, will likely give new insights in how M. tuberculosis survives oxidative stress in human macrophages.

The bacterial ribosome, a promising focus for structure-based drug design.

Recent crystal structures of the bacterial ribosome have identified the complex molecular interactions involved in antibiotic-ribosome recognition. Insights into the binding of aminoglycosides, macrolides, tetracyclines and other antibiotics provide opportunities for computational, structure-based approaches to be used in the design of appropriate modifications to existing antibiotics as well as in the discovery of completely new drug classes.

Intrinsic flexibility of B-DNA: the experimental TRX scale

B-DNA flexibility, crucial for DNA–protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR 31P chemical shifts in solution. These measurements reflect the BI ↔ BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

Importance of accurate DNA structures in solution: the Jun-Fos model.

Understanding the recognition of DNA sequences by proteins requires an accurate description of the structural dynamics of free DNA, especially regarding indirect readout. This involves subtle sequence-dependent effects that are difficult to characterize in solution. To progress in this area, we applied NMR and extensive simulations to a DNA sequence relevant to the Jun-Fos system. The backbone and base behaviors demonstrate that unrestrained simulations with major force fields (Parm98, Parmbsc0, and CHARMM27) are not reliable enough for in silico predictions of detailed DNA structures. More realistic structures required molecular dynamics simulations supplemented by NMR restraints. A new methodological element involved restraints inferred from the phosphate chemical shifts and from the phosphate dynamics. This provided a detailed and dynamic view of the intrinsic properties of the free DNA sequence that can be related to its recognition, by comparison with a relevant DNA-protein complex. We show how to exploit the relationship between phosphate motions and helicoidal descriptors for structure determination toward an accurate description of DNA structures and dynamics in solution.

Targeting conserved water molecules: design of 4-aryl-5-cyanopyrrolo[2,3-d]pyrimidine Hsp90 inhibitors using fragment-based screening and structure-based optimization.

Inhibitors of the Hsp90 molecular chaperone are showing promise as anti-cancer agents. Here we describe a series of 4-aryl-5-cyanopyrrolo[2,3-d]pyrimidine ATP competitive Hsp90 inhibitors that were identified following structure-driven optimization of purine hits revealed by NMR based screening of a proprietary fragment library. Ligand-Hsp90 X-ray structures combined with molecular modeling led to the rational displacement of a conserved water molecule leading to enhanced affinity for Hsp90 as measured by fluorescence polarization, isothermal titration calorimetry and surface plasmon resonance assays. This displacement was achieved with a nitrile group, presenting an example of efficient gain in binding affinity with minimal increase in molecular weight. Some compounds in this chemical series inhibit the proliferation of human cancer cell lines in vitro and cause depletion of oncogenic Hsp90 client proteins and concomitant elevation of the co-chaperone Hsp70. In addition, one compound was demonstrated to be orally bioavailable in the mouse. This work demonstrates the power of structure-based design for the rapid evolution of potent Hsp90 inhibitors and the importance of considering conserved water molecules in drug design.