Derivation of dihedral angles from CH-CH dipolar-dipolar cross-correlated relaxation rates: a C-C torsion involving a quaternary carbon atom in epothilone A bound to tubulin

Conformational Properties of the Chemotherapeutic Drug Analogue Epothilone A: How to Model a Flexible Protein Ligand Using Scarcely Available Experimental Data

Epothilones are among the most potent chemotherapeutic drugs used for the treatment of cancer. Epothilone A (EpoA), a natural product, is a macrocyclic molecule containing 34 non-hydrogen atoms and a thiazole side chain. NMR studies of EpoA in aqueous solution, unbound as well as bound to αβ-tubulin, and unbound in dimethyl sulfoxide (DMSO) solution have delivered sets of nuclear Overhauser effect (NOE) atom-atom distance bounds, but no structures based on NMR data are present in structural data banks. X-ray diffraction of crystals has provided structures of EpoA unbound and bound to αβ-tubulin. Since both crystal structures derived from X-ray diffraction intensities do not completely satisfy the three available sets of NOE distance bounds for EpoA, molecular dynamics (MD) simulations have been employed to obtain conformational ensembles in aqueous and in DMSO solution that are compatible with the respective NOE data. It was found that EpoA displays a larger conformational variability in DMSO than in water and the two conformational ensembles show little overlap. Yet, they both provide conformational scaffolds that are energetically accessible at physiological temperature and pressure.

Synthesis and Biological Evaluation of 7-Deoxy-Epothilone Analogues

The synthesis of two deoxygenated analogues of potent epothilones is reported in an effort to analyze the relative importance of molecular conformation and ligand-target interactions to biological activity. 7-deoxy-epothilone D and 7-deoxy-(S)-14-methoxy-epothilone D were prepared through total synthesis and shown to maintain the conformational preferences of their biologically active parent congeners through computer modeling and nuclear magnetic resonance (NMR) studies. The significant decrease in observed potency for each compound suggests that a hydrogen bond between the C7-hydroxyl group and the tubulin binding site plays a critical role in the energetics of binding in the epothilone class of polyketides.

Molecular recognition of epothilones by microtubules and tubulin dimers revealed by biochemical and NMR approaches

The binding of epothilones to dimeric tubulin and to microtubules has been studied by means of biochemical and NMR techniques. We have determined the binding constants of epothilone A (EpoA) and B (EpoB) to dimeric tubulin, which are 4 orders of magnitude lower than those for microtubules, and we have elucidated the conformation and binding epitopes of EpoA and EpoB when bound to tubulin dimers and microtubules in solution. The determined conformation of epothilones when bound to dimeric tubulin is similar to that found by X-ray crystallographic techniques for the binding of EpoA to the Tubulin/RB3/TTL complex; it is markedly different from that reported for EpoA bound to zinc-induced sheets obtained by electron crystallography. Likewise, only the X-ray structure of EpoA bound to the Tubulin/RB3/TTL complex at the luminal site, but not the electron crystallography structure, is compatible with the results obtained by STD on the binding epitope of EpoA bound to dimeric tubulin, thus confirming that the allosteric change (structuring of the M-loop) is the biochemical mechanism of induction of tubulin assembly by epothilones. TR-NOESY signals of EpoA bound to microtubules have been obtained, supporting the interaction with a transient binding site with a fast exchange rate (pore site), consistent with the notion that epothilones access the luminal site through the pore site, as has also been observed for taxanes. Finally, the differences in the tubulin binding affinities of a series of epothilone analogues has been quantitatively explained using the newly determined binding pose and the COMBINE methodology.

Epothilone D and its 9-Methyl analogues: combinatorial syntheses, conformation, and biological activities

Epothilone D (Epo D) and its 9-Methyl conformational analogues were synthesized through a highly efficient combinatorial approach. The fragment E was synthesized in 11 total steps with 6 longest linear steps, and each aldehyde B was prepared via a 3-step sequence. Starting from the common precursor E and a suitable aldehydes B, each target molecule were obtained in only 4 steps. The 9-(S)-epo D and 9-(R)-epo D demonstrated significant difference in inhibition activities against cancer cell lines and in conformational analysis.

The interaction of microtubules with stabilizers characterized at biochemical and structural levels

Since the discovery of paclitaxel and its peculiar mechanism of cytotoxicity, which has made it and its analogues widely used antitumour drugs, great effort has been made to understand the way they produce their effect in microtubules and to find other products that share this effect without the undesired side effects of low solubility and development of multidrug resistance by tumour cells. This chapter reviews the actual knowledge about the biochemical and structural mechanisms of microtubule stabilization by microtubule stabilizing agents, and illustrates the way paclitaxel and its biomimetics induce microtubule assembly, the thermodynamics of their binding, the way they reach their binding site and the conformation they have when bound.

Zampanolide, a potent new microtubule-stabilizing agent, covalently reacts with the taxane luminal site in tubulin α,β-heterodimers and microtubules

Zampanolide and its less active analog dactylolide compete with paclitaxel for binding to microtubules and represent a new class of microtubule-stabilizing agent (MSA). Mass spectrometry demonstrated that the mechanism of action of both compounds involved covalent binding to β-tubulin at residues N228 and H229 in the taxane site of the microtubule. Alkylation of N228 and H229 was also detected in α,β-tubulin dimers. However, unlike cyclostreptin, the other known MSA that alkylates β-tubulin, zampanolide was a strong MSA. Modeling the structure of the adducts, using the NMR-derived dactylolide conformation, indicated that the stabilizing activity of zampanolide is likely due to interactions with the M-loop. Our results strongly support the existence of the luminal taxane site of microtubules in tubulin dimers and suggest that microtubule nucleation induction by MSAs may proceed through an allosteric mechanism.

NMR in natural products: understanding conformation, configuration and receptor interactions

Covering: up to 2011. Natural products are of tremendous importance in both traditional and modern medicine. For medicinal chemistry natural products represent a challenge, as their chemical synthesis and modification are complex processes, which require many, often stereo-selective, synthetic steps. A prerequisite for the design of analogs of natural products, with more accessible synthetic routes, is the availability of their bioactive conformation. Nuclear Magnetic Resonance (NMR) spectroscopy and X-ray crystallography are the two techniques of choice to investigate the structure of natural products. In this review, I describe the most recent advances in NMR to study the conformation of natural products either free in solution or bound to their cellular receptors. In chapter 2, I focus on the use of residual dipolar couplings (RDC). On the basis of a few examples, I discuss the benefit of complementing classical NMR parameters, such as NOEs and scalar couplings, with dipolar couplings to simultaneously determine both the conformation and the relative configuration of natural products in solution. Chapter 3 is dedicated to the study of the structure of natural products in complex with their cellular receptors and is further divided in two sections. In the first section, I describe two solution-state NMR methodologies to investigate the binding mode of low-affinity ligands to macromolecular receptors. The first approach, INPHARMA (Interligand Noes for PHArmacophore Mapping), is based on the observation of interligand NOEs between two small molecules binding competitively to a common receptor. INPHARMA reveals the relative binding mode of the two ligands, thus allowing ligand superimposition. The second approach is based on paramagnetic relaxation enhancement (PRE) of ligand resonances in the presence of a receptor containing a paramagnetic center. In the second section, I focus on solid-state NMR spectroscopy as a tool to access the bioactive conformation of natural products in complex with macromolecular receptors.

An NMR-based scoring function improves the accuracy of binding pose predictions by docking by two orders of magnitude

Low-affinity ligands can be efficiently optimized into high-affinity drug leads by structure based drug design when atomic-resolution structural information on the protein/ligand complexes is available. In this work we show that the use of a few, easily obtainable, experimental restraints improves the accuracy of the docking experiments by two orders of magnitude. The experimental data are measured in nuclear magnetic resonance spectra and consist of protein-mediated NOEs between two competitively binding ligands. The methodology can be widely applied as the data are readily obtained for low-affinity ligands in the presence of non-labelled receptor at low concentration. The experimental inter-ligand NOEs are efficiently used to filter and rank complex model structures that have been pre-selected by docking protocols. This approach dramatically reduces the degeneracy and inaccuracy of the chosen model in docking experiments, is robust with respect to inaccuracy of the structural model used to represent the free receptor and is suitable for high-throughput docking campaigns.

Peptide ligands selected with CD4-induced epitopes on native dualtropic HIV-1 envelope proteins mimic extracellular coreceptor domains and bind to HIV-1 gp120 independently of coreceptor usage

During HIV-1 entry, binding of the viral envelope glycoprotein gp120 to the cellular CD4 receptor triggers conformational changes resulting in exposure of new epitopes, the highly conserved CD4-induced (CD4i) epitopes that are essential for subsequent binding to chemokine receptor CCR5 or CXCR4. Due to their functional conservation, CD4i epitopes represent attractive viral targets for HIV-1 entry inhibition. The aim of the present study was to select peptide ligands for CD4i epitopes on native dualtropic (R5X4) HIV-1 envelope (Env) glycoproteins by phage display. Using CD4-activated retroviral particles carrying Env from the R5X4 HIV-1 89.6 strain as the target, we performed screenings of random peptide phage libraries under stringent selection conditions. Selected peptides showed partial identity with amino acids in the extracellular domains of CCR5/CXCR4, including motifs rich in tyrosines and aspartates at the N terminus known to be important for gp120 binding. A synthetic peptide derivative (XD3) corresponding to the most frequently selected phages was optimized for Env binding on peptide arrays. Interestingly, the optimized peptide could bind specifically to gp120 derived from HIV-1 strains with different coreceptor usage, competed with binding of CD4i-specific monoclonal antibody (MAb) 17b, and interfered with entry of both a CCR5 (R5)-tropic and a CXCR4 (X4)-tropic Env pseudotyped virus. This peptide ligand therefore points at unique properties of CD4i epitopes shared by gp120 with different coreceptor usage and could thus serve to provide new insight into the conserved structural details essential for coreceptor binding for further drug development.

Precise structural determination of weakly binding peptides by utilizing dihedral angle constraints

Structural determination of target-bound conformations of peptides is of primary importance for the optimization of peptide ligands and peptide-mimetic design. In the structural determination of weakly binding ligands, transferred nuclear Overhauser effect (TrNOE) methods have been widely used. However, not many distance constraints can be obtained from small peptide ligands by TrNOE, especially for peptides bound to a target molecule in an extended conformation. Therefore, for precise structural determination of weakly binding peptides, additional structural constraints are required. Here, we present a strategy to systematically introduce dihedral angle constraints obtained from multiple transferred cross-correlated relaxation experiments and demonstrate precise structures of weakly binding peptides. As a result, we could determine the bioactive conformations of phage-derived peptide ligands and define their core binding motifs.

Novel epothilone lactones by an unusual diversion of the Grubbs' metathesis reaction

An unusual reaction with Grubbs' catalyst during the synthesis of bridged epothilones yielded five-membered internal lactones instead of the expected metathesis products. Three of the lactones have activities comparable to epothilone D.

The INPHARMA technique for pharmacophore mapping: A theoretical guide to the method

During the process of drug discovery, INPHARMA can be used to derive the structure of receptor/lead compound complexes binding to each other with a K(d) in the microM to mM range. To be successful, the methodology needs adjustment of various parameters that depend on the physical constants of the binding event and on the receptor size. Here we present a thorough theoretical analysis of the INPHARMA interligand NOE effect in dependence of experimental parameters and physical constants. This analysis helps the experimentalist to choose the correct experimental parameters and consequentially to achieve optimal performance of the methodology.

De(side chain) model of epothilone: bioconformer interconversions DFT study

Using ab initio methods, we have studied conformations of the de(sidechain)de(dioxy)difluoroepothilone model to quantify the effect of stability change between the exo and endo conformers of the epoxy ring. The DFT minimization of the macrolactone ring reveals four low energy conformers, although MP2 predicted five stable structures. The model tested with DFT hybride functional (B3LYP/6-31+G(d,p)) exhibits the global minimum for one of the exo forms (C), experimentally observed in the solid state, but unexpectedly with the MP2 electron correlation method for the virtual endo form (W). Using the QST3 technique, several pathways were found for the conversion of the low energy conformers to the other low energy exo representatives, as well as within the endo analog subset. The potential energy relationships obtained for several exo forms suggest a high conformational mobility between three, experimentally observed, conformers. The high rotational barrier, however, excludes direct equilibrium with experimental EC-derived endo form S. The highest calculated transition state for the conversion of the most stable exo M interligand to the endo S form is approximately a 28 kcal/mol above the energy of the former. The two-step interconversion of the exo H conformer to the endo S requires at least 28 kcal/mol. Surprisingly, we found that the transition state energy of the H form to the virtual endo W has the acceptable value of about 9 kcal/mol and the next energy barrier for free interconversion of endo W to endo S is 13 kcal/mol.

Design, synthesis and biological evaluation of bridged epothilone D analogues

Six epothilone D analogues with a bridge between the C4-methyl and the C12-methyl carbons were prepared in an attempt to constrain epothilone D to its proposed tubulin-binding conformation. Ring-closing metathesis (RCM) was employed as the key step to build the C4-C26 bridge. In antiproliferative assays in the human ovarian cancer (A2780) and prostate cancer (PC3) cell lines, and also in tubulin assembly assay, all these compounds proved to be less active than epothilone D.

Conformational preferences of natural and C3-modified epothilones in aqueous solution

The conformational properties of the microtubule-stabilizing agent epothilone A ( 1a) and its 3-deoxy and 3-deoxy-2,3-didehydro derivatives 2 and 3 have been investigated in aqueous solution by a combination of NMR spectroscopic methods, Monte Carlo conformational searches, and NAMFIS calculations. The tubulin-bound conformation of epothilone A ( 1a), as previously proposed on the basis of solution NMR data, was found to represent a significant fraction of the ensemble of conformations present for the free ligands in aqueous solution.

NMR-Based Strategies to Elucidate Bioactive Conformations of Weakly Binding Ligands

Key processes in molecular biology are regulated by interactions between biomolecules. Protein-proteinand protein-ligand interactions, e.g., in signal transduction pathways, rely on the subtle interactionsbetween atoms at the binding interface of the involved molecules. Because biomolecules often havemany interacting partners, these interactions are not necessarily strong. The study of molecularrecognition gives insight into the complex network of signaling in life and is the basis of structure-baseddrug design.In the situation where the interaction is weak, one of the traditional methods that can be appliedto obtain structural information (internuclear distances) of the bound ligand is the so-called transferredNOE (trNOE) method. Recently, it became possible to use transferred cross-correlated relaxation (trCCR)to directly measure dihedral angles. The combined use of these two techniques significantly improvesthe precision of the structure determination of ligands weakly bound to macromolecules.The application of these techniques will be discussed in detail for a peptide derived fromIKKβ bound to the protein NEMO that plays an important rolein the NFκB pathway.

NMR determination of the bioactive conformation of peloruside A bound to microtubules

We report here on the determination of the conformation of Peloruside A bound to biochemically stabilized microtubules, by using TR-NOESY NMR experiments. As a previous step, the conformation of the free molecule in water solution has also been deduced. Despite the large size of the ring, Peloruside A mainly adopts two conformations in water solution. A conformational selection process takes place, and the microtubules-bound conformer is one of those present in the water solution, different than that existing in chloroform medium. A model of the binding mode to tubulin has also been proposed, by docking the bioactive conformation of peloruside, which involves the alpha-tubulin monomer, in contrast with taxol, which binds to the beta-monomer.

Pairwise NMR experiments for the determination of protein backbone dihedral angle Phi based on cross-correlated spin relaxation

Novel cross-correlated spin relaxation (CCR) experiments are described, which measure pairwise CCR rates for obtaining peptide dihedral angles Phi. The experiments utilize intra-HNCA type coherence transfer to refocus 2-bond JNCalpha coupling evolution and generate the N(i)-Calpha(i) or C'(i-1)-Calpha(i) multiple quantum coherences which are required for measuring the desired CCR rates. The contribution from other coherences is also discussed and an appropriate setting of the evolution delays is presented. These CCR experiments were applied to 15N- and 13C-labeled human ubiquitin. The relevant CCR rates showed a high degree of correlation with the Phi angles observed in the X-ray structure. By utilizing these CCR experiments in combination with those previously established for obtaining dihedral angle Psi, we can determine high resolution structures of peptides that bind weakly to large target molecules.

Ligand-Target Interactions: What Can We Learn from NMR?

The conformation of the ligand in complex with a macromolecular target can be studied by nuclear magnetic resonance (NMR) in solution for both tightly and weakly forming complexes. In the weak binding regime (k(off) > 10(4) Hz), the structure of the bound ligand is accessible also for very large complexes (>100 kDa), which are not amenable to NMR studies in the tight binding regime. Here I review the state-of-the-art NMR methodology used for screening ligands and for the structural investigation of bound ligand conformations, in both tight and weak binding regimes. The advantages and disadvantages of each approach are critically described. The NMR methodology used to investigate transiently forming complexes has expanded considerably in the past few years, opening new possibilities for a detailed description of ligand-target interactions. Novel methods for the determination of the bound ligand conformation, in particular transferred cross-correlated relaxation, are thoroughly reviewed, and their advantages with respect to established methodology are discussed, using the epothilone-tubulin complex as a primary example.

Computational comparison of microtubule-stabilising agents laulimalide and peloruside with taxol and colchicine

Microtubule-stabilising agents laulimalide and peloruside have been compared with tubulin-interacting drugs paclitaxel and colchicine by different computational approaches. Docking and QSAR-based programs point to a favourable interaction with the beta tubulin paclitaxel binding site, although an additional, preferred binding site has been found at the alpha subunit of tubulin. All together provides a plausible rationalisation of the singular binding features of these microtubule stabilisers and paves the way for future structural studies.