A Test of the Single-Conformation Hypothesis in the Analysis of NMR Data for Small Polar Molecules: A Force Field Comparison

Ensemble determination by NMR data deconvolution

Nuclear magnetic resonance (NMR) is the spectroscopic technique of choice for determining molecular conformations in solution at atomic resolution. As solution NMR spectra are rich in structural and dynamic information, the way in which the data should be acquired and handled to deliver accurate ensembles is not trivial. This Review provides a guide to the NMR experiment selection and parametrization process, the generation of viable theoretical conformer pools and the deconvolution of time-averaged NMR data into a conformer ensemble that accurately represents a flexible molecule in solution. In addition to reviewing the key elements of solution ensemble determination of flexible mid-sized molecules, the feasibility and pitfalls of data deconvolution are discussed with a comparison of the performance of representative algorithms.

Going Viral: An Investigation into the Chameleonic Behaviour of Antiviral Compounds

The ability to adjust conformations in response to the polarity of the environment, i.e. molecular chameleonicity, is considered to be important for conferring both high aqueous solubility and high cell permeability to drugs in chemical space beyond Lipinski's rule of 5. We determined the conformational ensembles populated by the antiviral drugs asunaprevir, simeprevir, atazanavir and daclatasvir in polar (DMSO-d₆ ) and non-polar (chloroform) environments with NMR spectroscopy. Daclatasvir was fairly rigid, whereas the first three showed large flexibility in both environments, that translated into major differences in solvent accessible 3D polar surface area within each conformational ensemble. No significant differences in size and polar surface area were observed between the DMSO-d₆ and chloroform ensembles of these three drugs. We propose that such flexible compounds are characterized as "partial molecular chameleons" and hypothesize that their ability to adopt conformations with low polar surface area contributes to their membrane permeability and oral absorption.

Incorporating NOE-Derived Distances in Conformer Generation of Cyclic Peptides with Distance Geometry

Nuclear magnetic resonance (NMR) data from NOESY (nuclear Overhauser enhancement spectroscopy) and ROESY (rotating frame Overhauser enhancement spectroscopy) experiments can easily be combined with distance geometry (DG) based conformer generators by modifying the molecular distance bounds matrix. In this work, we extend the modern DG based conformer generator ETKDG, which has been shown to reproduce experimental crystal structures from small molecules to large macrocycles well, to include NOE-derived interproton distances. In noeETKDG, the experimentally derived interproton distances are incorporated into the distance bounds matrix as loose upper (or lower) bounds to generate large conformer sets. Various subselection techniques can subsequently be applied to yield a conformer bundle that best reproduces the NOE data. The approach is benchmarked using a set of 24 (mostly) cyclic peptides for which NOE-derived distances as well as reference solution structures obtained by other software are available. With respect to other packages currently available, the advantages of noeETKDG are its speed and that no prior force-field parametrization is required, which is especially useful for peptides with unnatural amino acids. The resulting conformer bundles can be further processed with the use of structural refinement techniques to improve the modeling of the intramolecular nonbonded interactions. The noeETKDG code is released as a fully open-source software package available at www.github.com/rinikerlab/customETKDG.

Nuclear magnetic resonance free ligand conformations and atomic resolution dynamics

Knowledge of free ligand conformational preferences (energy minima) and conformational dynamics (rotational energy barriers) of small molecules in solution can guide drug design hypotheses and help rank ideas to bias syntheses towards more active compounds. Visualization of conformational exchange dynamics around torsion angles, by replica exchange with solute tempering molecular dynamics (REST-MD), gives results in agreement with high-resolution 1 H nuclear magnetic resonance (NMR) spectra and complements free ligand conformational analyses. Rotational energy barriers around individual bonds are comparable between calculated and experimental values, making the in-silico method relevant to ranking prospective design ideas in drug discovery programs, particularly across a series of analogs. Prioritizing design ideas, based on calculations and analysis of measurements across a series, efficiently guides rational discovery towards the "right molecules" for effective medicines.

Cell Permeability of Isomeric Macrocycles: Predictions and NMR Studies

Conformation-dependent 3D descriptors have been shown to provide better predictions of the physicochemical properties of macrocycles than 2D descriptors. However, the computational identification of relevant conformations for macrocycles is nontrivial. Herein, we report that the Caco-2 cell permeability difference between a pair of diastereomeric macrocycles correlated with their solvent accessible 3D polar surface area and radius of gyration. The descriptors were calculated from the macrocycles’ solution-phase conformational ensembles and independently from ensembles obtained by conformational sampling. Calculation of the two descriptors for three other stereo- and regioisomeric macrocycles also allowed the correct ranking of their cell permeability. Methods for conformational sampling may thus allow ranking of passive permeability for moderately flexible macrocycles, thereby contributing to the prioritization of macrocycles for synthesis in lead optimization.

Solution-State Preorganization of Cyclic β-Hairpin Ligands Determines Binding Mechanism and Affinities for MDM2

Understanding mechanisms of protein folding and binding is crucial to designing their molecular function. Molecular dynamics (MD) simulations and Markov state model (MSM) approaches provide a powerful way to understand complex conformational change that occurs over long time scales. Such dynamics are important for the design of therapeutic peptidomimetic ligands, whose affinity and binding mechanism are dictated by a combination of folding and binding. To examine the role of preorganization in peptide binding to protein targets, we performed massively parallel explicit-solvent MD simulations of cyclic β-hairpin ligands designed to mimic the p53 transactivation domain and competitively bind mouse double minute 2 homologue (MDM2). Disrupting the MDM2-p53 interaction is a therapeutic strategy to prevent degradation of the p53 tumor suppressor in cancer cells. MSM analysis of over 3 ms of aggregate trajectory data enabled us to build a detailed mechanistic model of coupled folding and binding of four cyclic peptides which we compare to experimental binding affinities and rates. The results show a striking relationship between the relative preorganization of each ligand in solution and its affinity for MDM2. Specifically, changes in peptide conformational populations predicted by the MSMs suggest that entropy loss upon binding is the main factor influencing affinity. The MSMs also enable detailed examination of non-native interactions which lead to misfolded states and comparison of structural ensembles with experimental NMR measurements. In contrast to an MSM study of p53 transactivation domain (TAD) binding to MDM2, MSMs of cyclic β-hairpin binding show a conformational selection mechanism. Finally, we make progress toward predicting accurate off rates of cyclic peptides using multiensemble Markov models (MEMMs) constructed from unbiased and biased simulated trajectories.

Employing complementary spectroscopies to study the conformations of an epimeric pair of side-chain stapled peptides in aqueous solution

Understanding the conformational preferences of free ligands in solution is often necessary to rationalize structure-activity relationships in drug discovery. Herein, we examine the conformational behavior of an epimeric pair of side-chain stapled peptides that inhibit the FAD dependent amine oxidase lysine specific demethylase 1 (LSD1). The peptides differ only at a single stereocenter, but display a major difference in binding affinity. Their Raman optical activity (ROA) spectra are most likely dominated by the C-terminus, obscuring the analysis of the epimeric macrocycle. By employing NMR spectroscopy, we show a difference in conformational behavior between the two compounds and that the LSD1 bound conformation of the most potent compound is present to a measurable extent in aqueous solution. In addition, we illustrate that Molecular Dynamics (MD) simulations produce ensembles that include the most important solution conformations, but that it remains problematic to identify relevant conformations with no a priori knowledge from the large conformational pool. Furthermore, this work highlights the importance of understanding the scope and limitations of the available techniques for conducting conformational analyses. It also emphasizes the importance of conformational selection of a flexible ligand in molecular recognition.

Are Computational Methods Useful for Structure Elucidation of Large and Flexible Molecules? Belizentrin as a Case Study

Quantum mechanical NMR methods are progressively becoming decisive in structure elucidation. However, problems arise using low-level calculations for complex molecules, whereas methods using higher levels of theory are not practical for large molecules. This report outlines a synergistic effort employing computationally inexpensive quantum mechanical NMR calculations with conformer selection incorporating ³J_HH values as a way to solve the structure of large, complex, and highly flexible molecules using readily available computational resources with belizentrin as a case study.

Predicting the Permeability of Macrocycles from Conformational Sampling - Limitations of Molecular Flexibility

Macrocycles constitute superior ligands for targets that have flat binding sites but often require long synthetic routes, emphasizing the need for property prediction prior to synthesis. We have investigated the scope and limitations of machine learning classification models and of regression models for predicting the cell permeability of a set of denovo-designed, drug-like macrocycles. 2D-Based classification models, which are fast to calculate, discriminated between macrocycles that had low-medium and high permeability and may be used as virtual filters in early drug discovery projects. Importantly, stereo- and regioisomer were correctly classified. QSPR studies of two small sets of comparator drugs suggested that use of 3D descriptors, calculated from biologically relevant conformations, would allow development of more precise regression models for late phase drug projects. However, a 3D permeability model could only be developed for a rigid series of macrocycles. Comparison of NMR based conformational analysis with in silico conformational sampling indicated that this shortcoming originates from the inability of the molecular mechanics force field to identify the relevant conformations for flexible macrocycles. We speculate that a Kier flexibility index of ≤10 constitutes a current upper limit for reasonably accurate 3D prediction of macrocycle cell permeability.

Conformation of the Macrocyclic Drug Lorlatinib in Polar and Nonpolar Environments: A MD Simulation and NMR Study

The replica exchange molecular dynamics (REMD) simulation is demonstrated to readily predict the conformations of the macrocyclic drug lorlatinib, as validated by solution NMR studies. In aqueous solution, lorlatinib adopts a conformer identical to its target bound structure. This conformer is stabilized by an extensive hydrogen bond network to the solvents. In chloroform, lorlatinib populates two conformers with the second one being less polar, which may contribute to lorlatinib's ability to cross cell membranes.

Pyrrolidine ring puckering and prolyl amide bond configurations of 2-methyl-allo-hydroxyproline-based dipeptides

An expeditious method for the synthesis of homo and heterochiral dipeptides containing l-alanine and d/l 2-methyl allo-hydroxyl prolines was developed using direct aminolysis of bicyclic lactones derived from d/l alanine. The impact of C-2 methylation and its spatial orientation on the pyrrolidine ring puckering and prolyl amide bond configuration was ascertained by solution NMR studies. The present studies reveal that C-2 methylation causes the prolyl amide bond to exist exclusively in the trans geometry in both homo- and heterochiral dipeptides. However, the spatial orientation of the C-2 methyl group and its i + 2 position in appropriately capped model dipeptides may nucleate into a turn like structure.

Conformational Sampling of Macrocyclic Drugs in Different Environments: Can We Find the Relevant Conformations?

Conformational flexibility is a major determinant of the properties of macrocycles and other drugs in beyond rule of 5 (bRo5) space. Prediction of conformations is essential for design of drugs in this space, and we have evaluated three tools for conformational sampling of a set of 10 bRo5 drugs and clinical candidates in polar and apolar environments. The distance-geometry based OMEGA was found to yield ensembles spanning larger structure and property spaces than the ensembles obtained by MOE-LowModeMD (MOE) and MacroModel (MC). Both MC and OMEGA but not MOE generated different ensembles for polar and apolar environments. All three conformational search methods generated conformers similar to the crystal structure conformers for 9 of the 10 compounds, with OMEGA performing somewhat better than MOE and MC. MOE and OMEGA found all six conformers of roxithromycin that were identified by NMR in aqueous solutions, whereas only OMEGA sampled the three conformers observed in chloroform. We suggest that characterization of conformers using molecular descriptors, e.g., the radius of gyration and polar surface area, is preferred to energy- or root-mean-square deviation-based methods for selection of biologically relevant conformers in drug discovery in bRo5 space.

Flexibility is important for inhibition of the MDM2/p53 protein-protein interaction by cyclic β-hairpins

Protein-protein interactions that have large, flat and featureless binding sites are difficult drug targets. In the development of their modulators conventional drug discovery strategies are often unsuccessful. Gaining a detailed understanding of the binding mode of protein-protein interaction inhibitors is therefore of vast importance for their future pharmaceutical use. The MDM2/p53 protein pair is a highly promising target for cancer treatment. Disruption of the protein complex using p53 α-helix mimetics has been shown to be a successful strategy to control p53 activity. To gain further insight into the binding of inhibitors to MDM2, the flexibility of four cyclic β-hairpins that act as α-helical mimetics and potential MDM2/p53 interaction inhibitors was investigated in relation to their inhibitory activity. MDM2-binding of the mimetics was determined using fluorescence polarization and surface plasmon resonance assays, whereas their conformation and dynamics in solution was described by the combined experimental and computational NAMFIS analysis. Molecular flexibility was shown to be important for the activity of the cyclic β-hairpin based MDM2 inhibitors.

Halogen Bonding: A Powerful Tool for Modulation of Peptide Conformation

Halogen bonding is a weak chemical force that has so far mostly found applications in crystal engineering. Despite its potential for use in drug discovery, as a new molecular tool in the direction of molecular recognition events, it has rarely been assessed in biopolymers. Motivated by this fact, we have developed a peptide model system that permits the quantitative evaluation of weak forces in a biologically relevant proteinlike environment and have applied it for the assessment of a halogen bond formed between two amino acid side chains. The influence of a single weak force is measured by detection of the extent to which it modulates the conformation of a cooperatively folding system. We have optimized the amino acid sequence of the model peptide on analogues with a hydrogen bond-forming site as a model for the intramolecular halogen bond to be studied, demonstrating the ability of the technique to provide information about any type of weak secondary interaction. A combined solution nuclear magnetic resonance spectroscopic and computational investigation demonstrates that an interstrand halogen bond is capable of conformational stabilization of a β-hairpin foldamer comparable to an analogous hydrogen bond. This is the first report of incorporation of a conformation-stabilizing halogen bond into a peptide/protein system, and the first quantification of a chlorine-centered halogen bond in a biologically relevant system in solution.

Solution Conformations of Curcumin in DMSO

NAMFIS (NMR Analysis of Molecular Flexibility In Solution) has been applied to curcumin dissolved in DMSO. Quantitative ¹H-¹H distance constraints reduce a pool of candidate conformations to a solution collection of four enol conformations-two of these match curcumin crystallized with human transthyretin, and one is closely related to a single-crystal structure of curcumin.

Dynamic Docking of Conformationally Constrained Macrocycles: Methods and Applications

Many natural products consist of large and flexible macrocycles that engage their targets via multiple contact points. This combination of contained flexibility and large contact area often allows natural products to bind at target surfaces rather than deep pockets, making them attractive scaffolds for inhibiting protein-protein interactions and other challenging therapeutic targets. The increasing ability to manipulate such compounds either biosynthetically or via semisynthetic modification means that these compounds can now be considered as starting points for medchem campaigns rather than solely as ends. Modern medchem benefits substantially from rational improvements made on the basis of molecular docking. As such, docking methods have been enhanced in recent years to deal with the complicated binding modalities and flexible scaffolds of macrocyclic natural products and natural product-like structures. Here, we comprehensively review methods for treating and docking these large macrocyclic scaffolds and discuss some of the resulting advances in medicinal chemistry.

Conformational analysis: ³JHCOC and ³JHCCC Karplus relationships for methylene ¹H nuclei

NAMFIS (NMR Analysis of Molecular Flexibility In Solution) was applied to 1-[2-(benzyloxy)phenyl]ethanone using quantitative (1)H-(1)H NOE distances and (3)J proton-carbon coupling constant (CC) restraints for averaged methylene proton (3)J(HCOC) and (3)J(HCCC) pathways H2-(3)J-X imposed by density functional theory-generated Karplus relationships. Comparison of the NOE-only versus the NOE + CC conformational selections illustrates that the experimentally measured average (3)J coupling constants of methylene protons can be used for solution conformational analysis, potentially valuable in the study of small-molecule drugs and natural products which lack the typically studied H1-(3)J-X Karplus relationships.

Antiviral atropisomers: conformational energy surfaces by NMR for host-directed myxovirus blockers

Biologically active organic molecules characterized by a high single bond torsional barrier generate isolable isomers (atropisomers) and offer a unique stereochemical component to the design of selective therapeutic agents. The present work presents a nanomolar active inhibitor of myxoviruses, which most likely acts by blocking one or more cellular host proteins but also, serendipitously, exhibits axial chirality with an energy barrier of ΔG((++)) ≥30 kcal/mol. The latter has been probed by variable temperature NMR and microwave irradiation and by high level DFT transition state analysis and force field calculations. Full conformational profiles of the corresponding (aR,S) and (aS,S) atropisomers at ambient temperature were derived by conformer deconvolution with NAMFIS (NMR Analysis by Molecular Flexibility In Solution) methodology to generate seven and eight individual conformations, each assigned a % population. An accurate evaluation of a key torsion angle at the center of the molecules associated with a (3)JC-S-C-H coupling constant was obtained by mapping the S-C bond rotation with the MPW1PW91/6-31G-d,p DFT method followed by fitting the resulting dihedral angles and J-values to a Karplus expression. Accordingly, we have developed a complete conformational profile of diastereomeric atropisomers consistent with both high and low rotational barriers. We expect this assessment to assist the rationalization of the selectivity of the two (aR,S) and (aS,S) forms against host proteins, while offering insights into their divergent toxicity behavior.

Quantification of free ligand conformational preferences by NMR and their relationship to the bioactive conformation

Accurate unbound solution 3D-structures of ligands provide unique opportunities for medicinal chemistry and, in particular, a context to understand binding thermodynamics and kinetics. Previous methods of deriving these 3D-structures have had neither the accuracy nor resolution needed for drug design and have not yet realized their potential. Here, we describe and apply a NMR methodology to the aminoglycoside streptomycin that can accurately quantify accessible 3D-space and rank the occupancy of observed conformers to a resolution that enables medicinal chemistry understanding and design. Importantly, it is based upon conventional small molecule NMR techniques and can be performed in physiologically-relevant solvents. The methodology uses multiple datasets, an order of magnitude more experimental data than previous NMR approaches and a dynamic model during refinement, is independent of computational chemistry and avoids the problem of virtual conformations. The refined set of solution 3D-shapes for streptomycin can be grouped into two major families, of which the most populated is almost identical to the 30S ribosomal subunit bioactive shape. We therefore propose that accurate unbound ligand solution conformations may, in some cases, provide a subsidiary route to bioactive shape without crystallography. This experimental technique opens up new opportunities for drug design and more so when complemented with protein co-crystal structures, SAR data and pharmacophore modeling.

Dynamics of the glycosidic bond: conformational space of lactose

The dynamics of the glycosidic bond of lactose was studied by a paramagnetic tagging-based NMR technique, which allowed the collection of an unusually large series of NMR data for a single compound. By the use of distance- and orientation-dependent residual dipolar couplings and pseudocontact shifts, the simultaneous fitting of the probabilities of computed conformations and the orientation of the magnetic susceptibility tensor of a series of lanthanide complexes of lactose show that its glycosidic bond samples syn/syn, anti/syn and syn/anti ϕ/ψ regions of the conformational space in water. The analysis indicates a higher reliability of pseudocontact shift data as compared to residual dipolar couplings with the presently available weakly orienting paramagnetic tagging technique. The method presented herein allows for an improved understanding of the dynamic behaviour of oligosaccharides.

Conformational ensembles of flexible beta-turn mimetics in DMSO-d6

Beta-turns play an important role in peptide and protein chemistry, biophysics, and bioinformatics. The aim of this research was to study short linear peptides that have a high propensity to form beta-turn structures in solution. In particular, we examined conformational ensembles of beta-turn forming peptides with a general sequence CBz-L-Ala-L-Xaa-Gly-L-Ala-OtBu. These tetrapeptides, APGA, A(4R)MePGA, and A(4S)MePGA, incorporate proline, (4R)-methylproline, and (4S)-methylproline, respectively, at the Xaa position. To determine the influence of the 4-methyl substituted prolines on the beta-turn populations, the NAMFIS (NMR analysis of molecular flexibility in solution) deconvolution analysis for these three peptides were performed in DMSO-d(6) solution. The NBO (natural bond orbital) method was employed to gain further insight into the results obtained from the NAMFIS analysis. The emphasis in the NBO analysis was to characterize remote intramolecular interactions that could influence the backbone-backbone interactions contributing to beta-turn stability. NAMFIS results indicate that the enantiospecific incorporation of the methyl substituent at the C(gamma) (C4) position of the proline residue can be used to selectively control the pyrrolidine ring puckering propensities and, consequently, the preferred varphi,psi angles associated with the proline residue in beta-turn forming peptides. The NAMFIS analyses show that the presence of (4S)-methylproline in A(4S)MePGA considerably increased the type II beta-turn population with respect to APGA and A(4R)MePGA. The NBO calculations suggest that this observation can be rationalized based on an n-->pi* interaction between the N-terminus alanine carbonyl oxygen and the proline carbonyl group. Several other interactions between remote orbitals in these peptides provide a more detailed explanation for the observed population distributions.

The discodermolide hairpin structure flows from conformationally stable modular motifs

(+)-Discodermolide (DDM), a polyketide macrolide from marine sponge, is a potent microtubule assembly promoter. Reported solid-state, solution, and protein-bound DDM conformations reveal the unusual result that a common hairpin conformational motif exists in all three microenvironments. No other flexible microtubule binding agent exhibits such constancy of conformation. In the present study, we combine force-field conformational searches with NMR deconvolution in different solvents to compare DDM conformers with those observed in other environments. While several conformational families are perceived, the hairpin form dominates. The stability of this motif is dictated primarily by steric factors arising from repeated modular segments in DDM composed of the C(Me)-CHX-C(Me) fragment. Furthermore, docking protocols were utilized to probe the DDM binding mode in beta-tubulin. A previously suggested pose is substantiated (Pose-1), while an alternative (Pose-2) has been identified. SAR analysis for DDM analogues differentiates the two poses and suggests that Pose-2 is better able to accommodate the biodata.

Conformational properties of the macrocyclic trichothecene mycotoxin verrucarin A in solution

Phase-sensitive nuclear Overhauser enhancement spectroscopy (NOESY) experiments, (3)J couplings and computational molecular modeling (MM2* and MMFF force fields) were employed to examine the conformational properties of verrucarin A in chloroform solutions. The MMFF force field calculations resulted in a family of 12 low-energy structures along with their populations, the latter being determined by the NMR analysis of molecular flexibility in solution(NAMFIS) deconvolution analysis. The concluded model was capable of reproducing successfully the experimental NOESY cross-peak volumes and the proton-coupling constants. Among the 12 conformers, the one which was similar to the structure of verrucarin A in the solid state was the predominant accounting for 75% of the total relative population, although other low-energy conformations contributed to a lesser degree in order to explain the experimental data.

Evaluation of the tubulin-bound paclitaxel conformation: synthesis, biology, and SAR studies of C-4 to C-3' bridged paclitaxel analogues

The important anticancer drug paclitaxel binds to the beta-subunit of the alphabeta-tubulin dimer in the microtubule in a stoichiometric ratio, promoting microtubule polymerization and stability. The conformation of microtubule-bound drug has been the subject of intense study, and various suggestions have been proposed. In previous work we presented experimental and theoretical evidence that paclitaxel adopts a T-shaped conformation when it is bound to tubulin. In this study we report additional experimental data and calculations that delineate the allowable parameters for effective paclitaxel-tubulin interactions.

Structure activity relationships of 5-, 6-, and 7-methyl-substituted azepan-3-one cathepsin K inhibitors

The syntheses, in vitro characterizations, and rat and monkey in vivo pharmacokinetic profiles of a series of 5-, 6-, and 7-methyl-substituted azepanone-based cathepsin K inhibitors are described. Depending on the particular regiochemical substitution and stereochemical configuration, methyl-substituted azepanones were identified that had widely varied cathepsin K inhibitory potency as well as pharmacokinetic properties compared to the 4S-parent azepanone analogue, 1 (human cathepsin K, K(i,app) = 0.16 nM, rat oral bioavailability = 42%, rat in vivo clearance = 49.2 mL/min/kg). Of particular note, the 4S-7-cis-methylazepanone analogue, 10, had a K(i,app) = 0.041 nM vs human cathepsin K and 89% oral bioavailability and an in vivo clearance rate of 19.5 mL/min/kg in the rat. Hypotheses that rationalize some of the observed characteristics of these closely related analogues have been made using X-ray crystallography and conformational analysis. These examples demonstrate the potential for modulation of pharmacological properties of cathepsin inhibitors by substituting the azepanone core. The high potency for inhibition of cathepsin K coupled with the favorable rat and monkey pharmacokinetic characteristics of compound 10, also known as SB-462795 or relacatib, has made it the subject of considerable in vivo evaluation for safety and efficacy as an inhibitor of excessive bone resorption in rat, monkey, and human studies, which will be reported elsewhere.

Diastereomers of Dibromo-7-epi-10-deacetylcephalomannine: Crowded and Cytotoxic Taxanes Exhibit Halogen Bonds

The diastereomers of dibromo-7-epi-10-deacetylcephalomannine (6 and 7) have been isolated and characterized. Cytotoxicity and microtubule assembly assays demonstrate that cephalomannine analogue 6 possesses a potency profile very similar to that of Taxol, while isomer 7 is slightly less active. Solid state, solution, and tubulin-bound conformations of the two diastereomers were probed by using X-ray crystallography, 2-D NMR experiments in conjunction with the NAMFIS analysis, and the Glide docking protocol. In the crystal, isomer 7 exhibits an intermolecular halogen bond that may contribute to self-assembly. Neither crystal structure appears in the NAMFIS solution analysis, but both diastereomers are represented in solution by a T-shaped Taxol conformer. Glide docking demonstrates the latter to best fill the tubulin binding pocket, as has been shown for the parent Taxol drug. Each model of the bound complexes for 6 and 7 presents a single well-defined halogen bond from one of the ligand's bromines to Glu22 or Asp26 near the N-terminus of beta-tubulin, respectively. This first report of a halogen bond between taxanes and tubulin may prove useful in guiding the design and synthesis of other microtubule-stabilizing agents with a similar capacity.

Conformations of Laulimalide in DMSO-d6

Laulimalide is one of the newest naturally occurring macrolides known to act as a microtubule stabilizing agent with properties similar to Taxol. It also stands as being one of the most flexible with 18 rotatable bonds. This large number of rotatable bonds allows for approximately 3(18) potential conformers. To examine the conformational energy surface of laulimalide, we have performed an NAMFIS deconvolution analysis for laulimalide in DMSO-d6. The latter has been supplemented with a post-NAMFIS energy analysis at the Becke3LYP/6-31G level that examines the opposing effects of internal hydrogen bonding and syn-pentane interactions. In this way, we have identified 15 laulimalide conformations that can be classified into 5 different families: Supine, Convex, Cobra, Stretch, and Concave motifs.

Characterization of the preferred stereochemistry for the neuropharmacologic actions of antillatoxin

Antillatoxin is a potent ichthyotoxin and cytotoxin previously discovered from the marine cyanobacterium Lyngbya majuscula. Ensuing studies of its mechanism of action showed it to activate the mammalian voltage-gated sodium channel at a pharmacological site that is distinct from any previously described. The structure of antillatoxin, initially formulated from spectroscopic information, was subsequently corrected at one stereocenter (C-4) as a result of synthesis of four different antillatoxin stereoisomers (all possible C-4 and C-5 diastereomers). In the current study these four stereoisomers, (4R,5R)-, (4S,5R)-, (4S,5S)-, and (4R,5S)-antillatoxin, were characterized in five different biological assay systems: ichthyotoxicity to goldfish, microphysiometry using cerebellar granule cells (CGCs), lactose dehydrogenase efflux from CGCs, monitoring of intracellular Ca(2+) concentrations in CGCs, and cytotoxicity to Neuro 2a cells. Across these various biological measures there was great consistency in that the natural antillatoxin (the 4R,5R-isomer) was greater than 25-fold more potent than any of the other stereoisomers. Detailed NMR studies provided a number of torsion and distance constraints that were modeled using the MM2 force field to yield predicted solution structures of the four antillatoxin stereoisomers. The macrocycle and side chain of natural (4R,5R)-antillatoxin present an overall "L-shaped" topology with an accumulation of polar substituents on the external surface of the macrocycle and a hydrogen bond between N(H)-7' and the C(O)-1 carbonyl. The decreased potency of the three non-naturally occurring antillatoxin stereoisomers is certainly a result of their dramatically altered overall molecular topologies.

α-Turn Mimetics: Short Peptide α-Helices Composed of Cyclic Metallopentapeptide Modules

Alpha-Helices are key structural components of proteins and important recognition motifs in biology. Short peptides (<or=15 residues) corresponding to these helical sequences are rarely helical away from their stabilizing protein environments. New techniques for stabilizing short peptide helices could be valuable for studying protein folding, modeling proteins, creating artificial proteins, and may aid the design of inhibitors or mimics of protein function. This study reports the facile incorporation of 3- and 4-alpha turns in 10-15 residue peptides through formation in situ of multiple cyclic metallopeptide modules [Pd(en)(H*XXXH*)](2+). The nonhelical peptides Ac-H*ELTH*H*VTDH*-NH(2) (1), Ac-H*ELTH*AVTDYH*ELTH*-NH(2) (2), and Ac-H*AAAH*HELTH*H*VTDH*-NH(2) (3) (H is histidine-methylated at imidazole-N3) react in N,N-dimethylformamide (DMF) or water with 2, 2, and 3 molar equivalents, respectively, of [Pd(en)(NO(3))(2)] to form exclusively [Pd(2)(en)(2)(Ac-H*ELTH*H*VTDH*-NH(2))](4+) (4), [Pd(2)(en)(2)(Ac-H*ELTH*AVTDYH*ELTH*-NH(2))](4+) (5), and [Pd(3)(en)(3)(Ac-H*AAAH*HELTH*H*VTDH*-NH(2))](6+) (6), characterized by mass spectrometry, 1D and 2D (1)H- and 1D (15)N-NMR spectroscopy. Despite the presence of multiple histidines and other possible metal-binding residues in these peptides, 2D (1)H NMR spectra reveal that Pd(en)(2+) is remarkably specific in coordinating to imidazole-N1 of only (i, i + 4) pairs of histidines (i.e., only those separated by three amino acids), resulting in 4-6 made up of cyclic metallopentapeptide modules ([Pd(en)(H*XXXH*)](2+))(n), n = 2, 2, 3, respectively, each cycle being a 22-membered ring. We have previously shown that a single metallopentapeptide can nucleate alpha-helicity (Kelso et al., Angew. Chem., Int. Ed. 2003, 42, 421-424.). We now demonstrate its use as an alpha-turn-mimicking module for the facile conversion of unstructured short peptides into helices of macrocycles and provide 1D and 2D NMR spectroscopic data, structure calculations via XPLOR and NMR analysis of molecular flexibility in solution (NAMFIS), and CD spectra in support of the alpha-helical nature of these monomeric metallopeptides in solution.

Pyridinedicarboxamide Strands Form Double Helices via an Activated Slippage Mechanism

The intertwining process of two strands of oligo-pyridinecarboxamides to form a double helix (Nature 2000, 407, 720) is found to consist of a series of discrete steps, where the tail of one of the strands proceeds inside the other single helix in an eddy-like process. While a plethora of minima can be located along the pathway, they exist only for a few, well-defined supramolecular arrangements of the two molecules. The initial transition state for the introduction of one molecule in the pitch of the other has the largest barrier and is therefore the rate-determining step of an activated slippage mechanism, which is characterized by a series of roller-coasting hills. Along the entire pathway, the intramolecular energy that stabilizes the single helices is slowly transformed into intermolecular energy that finally provides the necessary stabilization only near the end of the entwining process. Solvent or other chemical factors, such as the presence of ions, able to destabilize the full formation of the double helix may therefore drastically affect its formation.

On the stability of a single-turn alpha-helix: the single versus multiconformation problem

The pentapeptide Ac-HAAAH-NH2, cyclized through its imidazoles by PdII to give [Pd(en)(peptide)]2+, has recently been evaluated by 2-D NMR and simulated annealing as a single alpha-helix conformation in solution. In the present work, we have questioned this assumption by developing Pd2+ parameters for AMBER*, performing an extensive conformational search for the [Pd(en)(peptide)]2+, and deconvoluting the averaged NMR data into eight rapidly equilibrating conformations with populations ranging from 2 to 55%. None of the latter correspond to the alpha-helix, although a 3% form possesses a related structure. As a critical component of interpreting an averaged NMR spectrum in terms of a single conformation, we advise testing this assumption with a method that permits conformational deconvolution.

The NH---FC dipole orientation effect for pendant exocyclic CH(2)F

[structure: see text] DFT and MMFF force field calculations for 2 (R = H) predict that two conformers dominate in water (>/=95%) and both sustain a geometry in which the C-F and H-N dipoles align oppositely in a near-planar arrangement. The (1)H NMR spectra (D(2)O and DMSO-d(6)) and X-ray structure for 2 (R = SO(2)NHEt) confirm the predictions in all essentials. A novel single-conformer system is proposed.

Synthesis and biological evaluation of novel macrocyclic paclitaxel analogues

[reaction: see text] This work describes the synthesis of two novel macrocyclic taxoid constructs by ring-closing olefin metathesis (RCM) and their biological evaluation. Computational studies examine conformational profiles of 1 and 2 for their fit to the beta-tubulin binding site determined by electron crystallography. The results support the hypothesis that paclitaxel binds to microtubules in a "T" conformation.

High-resolution structure and conformational dynamics of rigid, cofacially aligned porphyrin-bridge-quinone systems as determined by NMR spectroscopy and ab initio simulated annealing calculations

The high-resolution solution structure and dynamics of a cofacially aligned porphyrin--phenylene--quinone compound have been determined using (1)H NMR spectroscopy and simulated annealing calculations. Members of this class of pi-stacked assemblies feature a 1,8-naphthyl pillaring motif that enforces sub van der Waals interplanar separations between juxtaposed porphyryl, aromatic bridge, and quinonyl components of the donor--spacer--acceptor compound; this structural motif gives rise to a comprehensive set of structurally significant NOE signatures that can be used as constraints in quantitative structural calculations. Examination of such data using ab initio simulated annealing analytical methods shows that 5-[8'-(4' '-[8' "-(2' " ',5' "-benzoquinonyl)-1' "-naphthyl]-1' '-phenyl)-1'-naphthyl]-10,20-diphenylporphyrin displays an unusual degree of conformational homogeneity in the condensed phase, and represents a rare example where such an analysis determines unequivocally a single unique structure in solution.

Synthesis and NMR-driven conformational analysis of taxol analogues conformationally constrained on the C13 side chain

Analogues of Taxol (paclitaxel) with the side chain conformationally restricted by insertion of a carbon linker between the 2'-carbon and the ortho-position of the 3'-phenyl ring were synthesized. Biological evaluation of these new taxoids showed that activity was dependent on the length of the linker and the configuration at C2' and C3'. Two analogues in the homo series, 9a and 24a, showed tubulin binding and cytotoxicity comparable to that of Taxol. NAMFIS (NMR analysis of molecular flexibility in solution) deconvolution of the averaged 2-D NMR spectra for 9a yields seven conformations. Within the latter set, the hydrophobically collapsed "nonpolar" and "polar" classes are represented by one conformation each with predicted populations of 12-15%. The five remaining conformers, however, are extended, two of which correspond to the T-conformation (47% of the total population). The latter superimpose well with the recently proposed T-Taxol binding conformer in beta-tubulin. The results provide evidence for the existence of two previously unrecognized structural features that support Taxol-like activity: (1) a reduced torsion angle between C2' and C3' and (2) an orthogonal arrangement of the mean plane through C1', C2' and the 2'-hydroxyl and the 3'-phenyl plane, the latter ring bisected by the former plane. By contrast, epimerization at 2',3' and homologation of the tether to CH2-CH2 were both detrimental for activity. The decreased activity of these analogues is apparently due to configurational and steric factors, respectively.