Relationship among ligand conformations in solution, in the solid state, and at the Hsp90 binding site: geldanamycin and radicicol

The unknown effects of a receptor's environment on a ligand's conformation presents a difficult challenge in predicting feasible bioactive conformations, particularly if the receptor is ill-defined. The primary hypothesis of this work is that a structure's conformational ensemble in solution presents viable candidates for protein binding. The experimental solution profile can be achieved with the NAMFIS (NMR analysis of molecular flexibility in solution) method, which deconvolutes the average NMR spectrum of small flexible molecules into individual contributing conformations with varying populations. Geldanamycin and radicicol are structurally different macrocycles determined by X-ray crystallography to bind to a common site on the cellular chaperone heat shock protein 90 (Hsp90). Without benefit of a receptor structure, NAMFIS has identified the bioactive conformers of geldanamycin and radicicol in CDCl3 solution with populations of 4% and 21%, respectively. Conversely, docking the set of NAMFIS conformers into the unliganded proteins with GLIDE followed by MM-GBSA scoring reproduces the experimental crystallographic binding poses.     

Chemistry and biology of resorcylic acid lactones

While resorcylic acid lactones (RALs) have been known for a long time, the more recent discoveries that radicicol is a potent and selective HSP90 inhibitor while other members such as hypothemycin, LL-Z1640-2 and LL-783,277 are potent kinase inhibitors have stimulated a renewed interest in this family of natural products. The recent developments regarding the chemistry and biology of RALs are reviewed.     

Quantum chemical calculations and mutational analysis suggest heat shock protein 90 catalyzes trans-cis isomerization of geldanamycin

The affinity of geldanamycin (GA) for binding to heat shock protein 90 (HSP90) is 50- to 100-fold weaker than is the affinity of the structurally distinct natural product radicicol. X-ray crystallography shows that although radicicol maintains its free conformation when bound to HSP90, the conformation of GA is dramatically altered from an extended conformation with a trans amide bond to a kinked shape in which the amide group in the ansa ring has the cis configuration. We have performed ab initio quantum chemical calculations to demonstrate that the trans-cis isomeriztion of GA in solution is both kinetically and thermodynamically unfavorable. Thus, we propose that HSP90 catalyzes the isomerization of GA. We identify Ser113, a conserved residue outside the ATP binding pocket, as essential for the isomerization of GA. In support of this model, we show that radicicol binds equally well to both wild-type HSP90 and the Ser113 mutant, whereas the binding of GA to the Ser113 mutant is decreased significantly from its binding to wild-type HSP90. Based on this finding, a mechanism of keto-enol tautomerization of GA catalyzed by HSP90 is proposed. The added requirement of isomerization prior to tight binding may explain the enhanced binding affinity of GA for HSP90 in a cell extract versus in a purified form.     

Synthesis of resorcinylic macrocycles related to radicicol via ring-closing metathesis

Novel resorcinylic macrolides, for example, 17, 24, were prepared via ring-closing metathesis as analogues of the HSP90 inhibitor radicicol.     

Concise synthesis of pochonin A, an HSP90 inhibitor

[chemical reaction: see text]. An expedient synthesis of (-)-pochonin A is reported (seven steps). This natural product is closely related to radicicol and was shown to be a 90 nM inhibitor of HSP90.     

Inhibition of Hsp90 with Resorcylic Acid Macrolactones: Synthesis and Binding Studies

A series of resorcylic acid macrolactones, analogues of the natural product radicicol has been prepared by chemical synthesis, and evaluated as inhibitors of heat shock protein 90 (Hsp90), an emerging attractive target for novel cancer therapeutic agents. The synthesis involves acylation of an ortho‐toluic acid dianion, esterification, followed by a ring‐closing metathesis to form the macrocycle. Subsequent manipulation of the protected hydroxymethyl side chain allows access to a range of new analogues following deprotection of the two phenolic groups. Co‐crystallization of one of the new macrolactones with the N‐terminal domain of yeast Hsp90 confirms that it binds in a similar way to the natural product radicicol and to our previous synthetic analogues, but that the introduction of the additional hydroxymethyl substituent appears to result in an unexpected change in conformation of the macrocyclic ring. As a result of this conformational change, the compounds bound less favorably to Hsp90.     

Design, synthesis, and structure--activity relationships for chimeric inhibitors of Hsp90

Inhibition of the 90 kDa heat shock protein (Hsp90) family of molecular chaperones represents a promising new chemotherapeutic approach toward the treatment of several cancers. Previous studies have demonstrated that the natural products, radicicol and geldanamycin, are potent inhibitors of the Hsp90 N-terminal ATP binding site. The cocrystal structures of these molecules bound to Hsp90 have been determined, and through molecular modeling and superimposition of these ligands, hybrids of radicicol and geldanamycin have been designed. A series of macrocylic chimeras of radicicol and geldanamycin and the corresponding seco-agents have been prepared and evaluated for both antiproliferative activity and their ability to induce Hsp90-dependent client protein degradation.     

ALADDIN: An integrated tool for computer-assisted molecular design and pharmacophore recognition from geometric,steric, and substructure searching of three-dimensional molecular structures

Summary ALADDIN is a computer program for the design or recognition of compounds that meet geometric, steric, and substructural criteria. ALADDIN searches a database of three-dimensional structures, marks atoms that meet substructural criteria, evaluates geometric criteria, and prepares a number of files that are input for molecular modification and coordinate generation as well as for molecular graphics. Properties calculated from the three-dimensional structure are described by either properties calculated from the molecule itself or from the molecule as compared to a reference molecule and associated surfaces. ALADDIN was used to design analogues to probe a bioactive conformation of a small molecule and a peptide, to test alternative superposition rules for receptor mapping of the D2 dopamine receptor, to recognize unexpected D2 dopamine agonist activity of existing compounds, and to design compounds to fit a binding site on a protein of known structure. We have found that series designed by ALADDIN show much more subtle variation in shape than do those designed by traditional methods and that compounds can be designed to be very close matches to the objective.     

Flexible matching of test ligands to a 3D pharmacophore using a molecular superposition force field: Comparison of predicted and experimental conformations of inhibitors of three enzymes

Summary A computer procedure TFIT, which uses a molecular superposition force field to flexibly match test compounds to a 3D pharmacophore, was evaluated to find out whether it could reliably predict the bioactive conformations of flexible ligands. The program superposition force field optimizes the overlap of those atoms of the test ligand and template that are of similar chemical type, by applying an attractive force between atoms of the test ligand and template which are close together and of similar type (hydrogen bonding, charge, hydrophobicity). A procedure involving Monte Carlo torsion perturbations, followed by torsional energy minimization, is used to find conformations of the test ligand which cominimize the internal energy of the ligand and the superposition energy of ligand and template. The procedure was tested by applying it to a series of flexible ligands for which the bioactive conformation was known experimentally. The 15 molecules tested were inhibitors of thermolysin, HIV-1 protease or endothiapepsin for which X-ray structures of the bioactive conformation were available. For each enzyme, one of the molecules served as a template and the others, after being conformationally randomized, were fitted. The fitted conformation was then compared to the known binding geometry. The matching procedure was successful in predicting the bioactive conformations of many of the structures tested. Significant deviation from experimental results was found only for parts of molecules where it was readily apparent that the template did not contain sufficient information to accurately determine the bioactive conformation.     

Concise asymmetric syntheses of radicicol and monocillin I.

Radicicol (1) exhibits potent anticancer properties in vitro, which are likely to be mediated through its high affinity (20 nM) for the molecular chaperone Hsp90. Recently, we reported the results of a synthetic program targeting radicicol (1) and monocillin I (2), highlighted by the application of ring-closing metathesis to macrolide formation. These efforts resulted in a highly convergent synthesis of radicicol dimethyl ether but failed in the removal of the two aryl methyl ethers. Simple exchange of these methyl ethers with more labile functionalities disabled a key esterification in the initial route. Through extended experimentation, a successful route to both natural products was secured, along with some intriguing results that emphasize the implications of this design on a broad range of fused benzoaliphatic targets, including analogues of these natural products.     

Design, synthesis, and computational studies of inhibitors of Bcl-XL

One of the primary objectives in the design of protein inhibitors is to shape the three-dimensional structures of small molecules to be complementary to the binding site of a target protein. In the course of our efforts to discover potent inhibitors of Bcl-2 family proteins, we found a unique folded conformation adopted by tethered aromatic groups in the ligand that significantly enhanced binding affinity to Bcl-XL. This finding led us to design compounds that were biased by nonbonding interactions present in a urea tether to adopt this bioactive, folded motif. To characterize the key interactions that induce the desired conformational bias, a series of substituted N,N'-diarylureas were prepared and analyzed using X-ray crystallography and quantum mechanical calculations. Stabilizing pi-stacking interactions and destabilizing steric interactions were predicted to work in concert in two of the substitution patterns to promote the bioactive conformation as a global energy minimum and result in a high target binding affinity. Conversely, intramolecular hydrogen bonding present in the third substitution motif promotes a less active, extended conformer as the energetically favored geometry. These findings were corroborated when the inhibition constant of binding to Bcl-XL was determined for fully elaborated analogues bearing these structural motifs. Finally, we obtained the NMR solution structure of the disubstituted N,N'-diarylurea bound to Bcl-XL demonstrating the folded conformation of the urea motif engaged in extensive pi-interactions with the protein.     

Development of a conformational search strategy for flexible ligands: A study of the potent μ-selective opioid analgesic fentanyl

Summary An extensive conformational search of the potent opioid analgesic, fentanyl, was performed using the semiempirical quantum mechanical method AM1 and the CHARMm potential energy function. A combination of two procedures was used to search the conformational space for fentanyl, which included nested dihedral scans, geometry optimization and molecular dynamics simulation at different temperatures. In addition, the effect of a continuum solvent environment was taken into account by use of appropriate values for the dielectric constant in the CHARMm computations.The results of the conformational search allowed the determination of the probable conformation of fentanyl in polar and nonpolar solvents and of three candidate conformers for its bioactive form.     

Compare-Conformer: a program for the rapid comparison of molecular conformers based on interatomic distances and torsion angles.

A computer program for comparison of the conformations of a number of related molecular structures is described. The comparisons are performed on either interatomic distances or torsion angles. The comparisons are accomplished on ordered pairs of distances or torsion angles, and the distance comparisons can be performed in a manner that allows permutation of the distance pairs being compared. The algorithm utilizes bit-string Boolean operations that allow the comparisons to be performed rapidly. The program should be useful for computer-assisted molecular modeling studies in which the viable conformers of bioactive analogues are compared in order to locate those conformers that place key substituents in the same spatial orientation.     

Novel inhibitor discovery through virtual screening against multiple protein conformations generated via ligand-directed modeling: a maternal embryonic leucine zipper kinase example

Kinase targets have been demonstrated to undergo major conformational reorganization upon ligand binding. Such protein conformational plasticity remains a significant challenge in structure-based virtual screening methodology and may be approximated by screening against an ensemble of diverse protein conformations. Maternal embryonic leucine zipper kinase (MELK), a member of serine-threonine kinase family, has been recently found to be involved in the tumerogenic state of glioblastoma, breast, ovarian, and colon cancers. We therefore modeled several conformers of MELK utilizing the available chemogenomic and crystallographic data of homologous kinases. We carried out docking pose prediction and virtual screening enrichment studies with these conformers. The performances of the ensembles were evaluated by their ability to reproduce known inhibitor bioactive conformations and to efficiently recover known active compounds early in the virtual screen when seeded with decoy sets. A few of the individual MELK conformers performed satisfactorily in reproducing the native protein-ligand pharmacophoric interactions up to 50% of the cases. By selecting an ensemble of a few representative conformational states, most of the known inhibitor binding poses could be rationalized. For example, a four conformer ensemble is able to recover 95% of the studied actives, especially with imperfect scoring function(s). The virtual screening enrichment varied considerably among different MELK conformers. Enrichment appears to improve by selection of a proper protein conformation. For example, several holo and unliganded active conformations are better to accommodate diverse chemotypes than ATP-bound conformer. These results prove that using an ensemble of diverse conformations could give a better performance. Applying this approach, we were able to screen a commercially available library of half a million compounds against three conformers to discover three novel inhibitors of MELK, one from each template. Among the three compounds validated via experimental enzyme inhibition assays, one is relatively potent (15; K(d) = 0.37 μM), one moderately active (12; K(d) = 3.2 μM), and one weak but very selective (9; K(d) = 18 μM). These novel hits may be utilized to assist in the development of small molecule therapeutic agents useful in diseases caused by deregulated MELK, and perhaps more importantly, the approach demonstrates the advantages of choosing an appropriate ensemble of a few conformers in pursuing compound potency, selectivity, and novel chemotypes over using single target conformation for structure-based drug design in general.     

Stabilization of zwitterions in solution: phosphinic and phosphonic acid GABA analogues

The influence of treatment of electron correlation, size of basis set and choice of solvation model on the predicted stabilization of zwitterionic phosphinic and phosphonic acid gamma-aminobutyric acid (GABA) analogues is investigated using ab initio molecular orbital methods and density functional theory. Density functional theory with the B3LYP functional and a composite basis set composed of the 6-31+G(2df) basis for phosphorus atoms and the 6-31+G(d,p) basis set for all other atoms is found to give an acceptable tradeoff between accuracy and computational expense. Either directly optimizing zwitterionic conformers within the conductor-like screening solvation model (COSMO) or pre-optimizing as dihydrates and then applying the COSMO solvation model give an acceptable treatment of solvation in terms of determining stable solvated structures, although directly optimizing within COSMO is simpler and less computationally expensive. With this protocol, cis-constrained phosphinic and phosphonic acid GABA analogues, which exhibit lower affinities for GABAC receptors, are found to possess only folded, intramolecularly hydrogen bonded low energy conformers. Trans-constrained analogues, on the other hand, exhibit higher affinities for GABAC receptors and are found to exist only as partially folded stable conformers. Conformationally flexible analogues can attain folded, partially folded or fully extended conformations and also have high biological activity. These results imply that the partially folded conformation is likely to be the most biologically active.     

Conformational analysis of TMC114, a novel HIV-1 protease inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.     

群多普利溶液构象的NMR研究

采用1D和2D NMR技术对血管紧张素转化酶抑制剂群多普利在CDCl₃溶液中存在的两种构象进行了结构解析, 并结合分子动力学和密度泛函理论方法对其进行了结构的几何优化和能量计算. 结果表明, 群多普利因分子中酰胺键的旋转而形成反式构象A和顺式构象B, 两种构象的能量差为6.35 kJ/mol, 且顺式构象为该药物的优势构象.     

Development of synthetic routes to macrocyclic compounds based on the HSP90 inhibitor radicicol

Short routes are reported to novel macrolides (e.g., 9, 12, 20) related to the HSP90 inhibitor radicicol.     

Solution-state conformational study of the hevamine inhibitor allosamidin and six potential inhibitor analogues by NMR spectroscopy and molecular modeling

The solution-state conformations of the hevamine inhibitor allosamidin and six potential inhibitor analogues were studied by various NMR spectroscopic techniques and molecular modeling using force field calculations. Determination solely of the global energy minimum conformation was found to be insufficient for consensus with the NMR results, and agreement between the NMR experimental data and the theoretical calculations was only reached by assessing the structures as population-weighted average conformers on the basis of Boltzmann distributions derived from the calculated relative energies. The conformations of the glycosidic linkages in the compounds were found to be similar when the sugar residues were the same, but differences were markedly evident otherwise and also for the various heterocyclic group linkages. The binding of the compounds to hevamine, which may also complex to chitinases in general, was assessed using HMQC, transfer-NOESY, and both 1-D and 2-D saturation transfer difference NMR experiments. Under the conditions employed, only allosamidin was implicated to be bound to hevamine, and then only by HMQC with the dipolar coupling-based experiments failing to substantiate the formation of the complex. However, the results are consistent with the biochemical activities of the compounds whereby only allosamidin has been shown to act as a competitive inhibitor.