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.     

2D H and C NMR conformational studies of thienopyridines and carboline biarylic compounds

In this paper, we report on conformational studies of biarylic compounds, as prepared through the well-known aza-Wittig methodology. The conformational studies were mainly realized by bidimensional (2D) nuclear magnetic resonance spectroscopy (NMR) and NOESY experiments. The conformational behavior showed that these biarylic compounds display an orthogonal symmetry and adopt a characteristic arrangement around the pivotal bond. Molecular modeling calculations were performed to support structure conformations.     

Conformational studies on (17alpha,20Z)-21-(X-Phenyl)-19-norpregna-1, 3,5(10),20-tetraene-3,17beta-diols using 1D and 2D NMR spectroscopy and GIAO calculations of (13)C shieldings

Differences in agonist responses of the novel estrogen receptor ligands (17alpha,20Z)-(p-methoxyphenyl)vinyl estradiol (1), (17alpha, 20Z)-(o-alpha,alpha,alpha-trifluoromethylphenyl)vinyl estradiol (2), and (17alpha,20Z)-(o-hydroxymethylphenyl)vinyl estradiol (3) led us to investigate their solution conformation. In competitive binding assay studies, we observed that several phenyl-substituted (17alpha, 20E/Z)-(X-phenyl)vinyl estradiols exhibited significant estrogen receptor binding, but with variation (RBA (1) = 20; RBA (2) = 23; RBA (3) = 140 where estradiol RBA = 100) depending on the phenyl substitution pattern. Because the 17alpha-phenylvinyl substituent interacts with the key helix-12 of the ligand binding domain, we considered that differences in the preferred conformation of 1-3 could account for their varying binding affinity. 2D NMR experiments at 500 MHz allowed the complete assignment of the (13)C and (1)H spectra of 1-3. The conformations of these compounds in solution were established by 2D and 1D NOESY spectroscopy. A statistical approach of evaluating contributing conformers of 1-3 from predicted (13)C shifts correlated quite well with the NOE data. The 17alpha substituents of 1 and 2 exist in similar conformational equilibria with some differences in relative populations of conformers. In contrast, the 17alpha substituent of 3 exists in a different conformational equilibrium. The similarity in solution conformations of 1 and 2 suggests they occupy a similar receptor volume, consistent with similar RBA values of 20 and 23. Conversely, the different conformational equilibria of 3 may contribute to the significant binding affinity (RBA = 140) of this ligand.     

Conformational and configurational analysis of an N,N carbonyl dipyrrinone-derived oximate and nitrone by NMR and quantum chemical calculations

The geometries and relative energies of new N,N carbonyl dipyrrinone-derived oxime molecules (E/Z-s-cis 4a and E/Z-s-cis 4b) have been investigated. The calculated energies, molecular geometries, and (1) H/(13) C NMR chemical shifts agree with experimental data, and the results are presented herein. The E-s-cis conformations of 4a and 4b and the Z-s-cis conformation of 5b were found to be the thermodynamically most stable isomers with the oxime hydrogen atom or the methyl functional group adopting an anti-orientation with respect to the dipyrrinone group. This conformation was unambiguously supported by a number of 2D NMR experiments.     

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.     

Partially O-alkylated thiacalix[4]arenes: synthesis, molecular and crystal structures, conformational behavior

NMR spectroscopy, X-ray diffraction analysis, and quantum chemical calculations were used for conformational behavior study of partially alkylated thiacalix[4]arenes bearing methyl (1), ethyl (2), or propyl (3) groups at the lower rim. The conformational properties are governed by two basic effects: (i) stabilization by intramolecular hydrogen bonds, and (ii) sterical requirements of the alkoxy groups at the lower rim. While the monosubstituted derivatives 1a and 3a adopt the cone conformation in solution, distally disubstituted compounds 1b, 1'b, 2b, 2'b, 3b, and 3'b exhibit several interesting conformational features. They prefer pinched cone conformation in solution, and, except for 3'b, they form also 1,2-alternate conformation, which is flexible and undergoes rather fast transition between two identical structures. The crystal structures of the compounds 1b, 2b, 2'b, and 3b revealed yet quite rare 1,2-alternate conformation forming molecular channels held together by pi-pi interactions. Different channels-with hexagonal symmetry, 0.26 nm wide-are formed in the crystal structure of the pinched cone conformation of 3b. An uncommon hydrogen bonding pattern was found in dimethoxy and dipropoxy derivatives 1'b and 3'b that adopt distorted cone conformations in crystal. Trialkoxy-substituted compounds 1c and 3c adopt the partial cone conformation in solution. A higher mobility of methyl derivative 1c enables also existence of the cone conformer.     

Conformation and stereodynamics of 2,2'-disubstituted N,N'-diaryl ureas

Except in the most hindered of cases, N,N'-diaryl N,N'-dimethyl ureas adopt a conformation with the two aryl rings disposed cis to one another. Variable temperature NMR studies reveal the rate at which the Ar-N bonds rotate as well as the conformational preference of ortho disubstituted ureas in which more than one cis orientation is possible. In general, a conformation in which the aryl rings lie close in space but with their most bulky 2-substituents aligned anti is preferred, but with particularly bulky 2-substituents, conformations in which one of the aryl rings points away from the other may also be populated.     

Konformative Beweglichkeit Flexibler Ringsysteme—XIII: Untersuchungen mit Hilfe der Protonenresonanzspektroskopie die Konformation des Ungesättigten Siebenringes im Benzocyclohepten und in symmetrisch substituierten derivaten

The confromations of the unsaturated seven membered ring in 4,4,6,6-tetradeuterium-1,2-benzocycloheptene-(1) ( 1 ) and five benzocycloheptene derivatives were determined by NMR spectroscopy. For all investigated compounds at ?80°C only one conformer was present in detectable quantity. By analysis of the NMR data – molecular symmetry, coupling constants and chemical shift – it can be shown that the conformation is always the chair form. The free conformational enthalpy of both the other conformations with boat or twist form of the ring is for all six compounds more than 1.8 kcal/mole. The experimental results agree with those from model calculations: thus for benzocycloheptene, the 5,5-dimethyl derivative ( 2 ) and the 4,4,6,6-tetramethyl derivative ( 4 ) the lowest energy was found for the chair conformation; the second most stable conformations were found to be the boat for 1 and 4 , and the twist form for 2 .     

Synthesis and NMR spectral studies of N-chloroacetyl-2,6-diarylpiperidin-4-ones

A series of 2,6-diarylpiperidin-4-ones having electron withdrawing chloroacetyl group at the heterocyclic nitrogen were synthesized. Unambiguous characterizations of the synthesized compounds were achieved by one-dimensional ((1)H NMR and (13)C NMR) and two-dimensional (HOMOCOSY, NOESY and HSQC spectra for compounds 8 and 9 and HOMOCOSY spectrum only for 10) NMR spectroscopic data. The conformational preferences of N-chloroacetyl-2,6-diarylpiperidin-4-ones with and without alkyl substituent at C-3 and C-5 (8-14) have also been discussed using the spectral studies. The spectral data and extracted coupling constant values suggest that the compounds 8, 12 and 14 adopt flattened boat conformation whereas the remaining compounds exist in twist-boat conformations in solution with coplanar orientation of the chloroacetyl moiety present at the heterocyclic nitrogen. The substituent parameters for the chloroacetyl moiety on the heterocyclic ring carbons have also been derived and discussed elaborately on the basis of their steric, electronic and gamma-eclipsing interaction. This substituent at the nitrogen causes a substantial change on the chemical shifts of ring carbons and the associated protons.     

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.     

NMR‐based conformational analysis of perezone and analogues

Complete assignment of the ¹H NMR chemical shift and coupling constant values of perezone (1), O‐methylperezone (2) and 6‐hydroxyperezone (3) was carried out by total‐line‐shape‐fitting calculations using the PERCH iterative spectra analysis software (PERCH Solutions Ltd., Kuopio, Finland). The resulting simulated spectra for the three compounds showed strong similarity to their corresponding experimental spectra. Particularly, all vicinal, allylic and homoallylic coupling constant values for the side chain of the three compounds were very similar, thus revealing that the conformation of these three molecules in solution is indeed almost identical. This fact is in agreement with extended side chain conformations over folded chain conformations because 1, 2 and 3 undergo completely different intramolecular cycloaddition reactions. In addition, results of double pulsed field gradient spin echo NOESY 1D experiments performed on perezone (1) were unable to provide evidence for folded conformers. Copyright © 2013 John Wiley & Sons, Ltd.     

Molecular modelling prediction of ligand binding site flexibility

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.     

The origin of the conformational preference of N,N'-diaryl-N,N'-dimethyl ureas

Poly aromatic ureas and poly aromatic amides are important classes of foldamers-oligomers with well defined conformations. We have explored the origins of the conformational preference of some N,N'-diaryl-N,N'-dimethyl ureas by a combination of NMR spectroscopy and electronic structure calculations using both a recently developed density functional (M06-2X) and a DFT approach (DFT-D) having empirical corrections for dispersive interactions. We have validated the DFT-D approach for structures of this type using high level wavefunction calculations, (CCSD(T)), of the unsubstituted N,N'-diphenyl-N,N'-dimethyl urea. For the N,N'-diaryl-N,N'-dimethyl ureas we have identified a number of 'endo' conformers (i.e. having an E,E geometrical conformation about the two urea C-N bonds), both π- and tert-butyl-stacked, as well as 'exo' structures (having a Z geometrical conformation about at least one of the C-N bonds), and have computed the relative energies of these conformers as well as the barriers for their interconversion. We find that the relative energies of the 'endo' structures closely follow the relative values of the dispersive interactions. The calculations have allowed us to associate different conformers with the various peaks in the NMR spectra, which point to relatively small differences in energy between the conformers. Somewhat larger energy differences are predicted by the two computational approaches, with the M06-2X functional performing the better of the two. It is suggested that the continuum model employed may not be sufficiently accurate to reflect the solvation of the various conformers.     

NMR Conformational Analysis and Theoretical Calculations for 2‐Aryl‐ 1,3‐dihydroxy‐4,4,5,5‐tetramethylimidazolidines

Conformational studies of 1,3‐dihydroxy‐4,4,5,5‐tetramethyl‐2‐(pyridin‐1‐yl)imidazolidine ( 1a ) and 1,3‐dihydroxy‐4,4,5,5‐tetramethyl‐2‐(pyridin‐3‐yl)imidazolidine ( 1b ), carried out by using 1D ¹H‐ and ¹³C‐NMR and 2D HMQC, HMBC, and NOESY experiments and with the aid of theoretical calculations, indicate that the OH groups are trans to the pyridinyl substituent. Because the two ¹H‐NMR signals of the Me groups are distinguishable and do not change between 290 and 380 K, it is proposed that 1a and 1b have each only one conformation in this temperature range. This behavior was not found with 1,3‐dihydroxy‐4,4,5,5‐tetramethyl‐2‐(pyridin‐2‐yl)imidazolidine ( 1c ) because its Me ¹H‐NMR signals cross over at 300 K. Hence, more than one conformation must be present, beyond those produced by simple inversions. Theoretical calculations including temperature and solvent effects were performed to provide further information on the conformational analysis and to help to assign the NMR data. The combination of NMR measurements and quantum‐chemical calculations is shown to be a very promising strategy for conformational analysis studies in solution.     

Systematic approach to understanding macrolide-ribosome interactions: NMR and modeling studies of oleandomycin and its derivatives

The three-dimensional structures of oleandomycin (1) and its derivatives oleandomycin-9-oxime (2) and 10,11-anhydrooleandomycin (3) were determined in different solvents by the combined use of NMR and molecular modeling methods. The experimental NMR data were compared with the results of molecular modeling and known crystal structures of the related molecules. It was shown that the dominant conformation of the lactone ring is the folded-out conformation with some amounts of the folded-in one depending on the solvent and temperature, while desosamine and cladinose sugars adopt the usual chair conformations. Modeling calculations provided evidence for conformational changes in the upper lactone region as well. Saturation transfer difference (STD) NMR experiments have provided information on the binding epitopes of 1-3 in complexes with E. coli ribosomes. The obtained molecular surfaces in close contact with ribosomes were compared with recently available 3D structures of the related macrolide-ribosome complexes, and the observed differences were discussed. The knowledge gained from this study can serve as a platform for the design of novel macrolides with an improved biological profile.     

Conformational changes of functionalised indole receptors upon their interaction with anions

The conformational analysis of four C2-amido and C7-ureido functionalised indole anion receptors was performed by a combination of heteronuclear NMR spectroscopy and ab initio quantum mechanical calculations. NOE experiments showed that anti–anti conformation across C2–C2α and C7–N7α bonds is predominant in acetone solution in the absence of anions. Upon anion binding to receptors, syn–syn conformation becomes predominant. The conformational changes upon anion binding are in good agreement with energetic preferences established by ab initio calculations. Chemical shift changes induced by interaction of anions suggest that binding of chloride and bromide anions occurs primarily to H1 and H7α protons. Nitrate anions favour interaction with H7α and H7γ ureido protons, whereas acetate anions interact strongly with all four available hydrogen bond donor groups.     

Comprehensive structural studies of 2',3'-difluorinated nucleosides: comparison of theory, solution, and solid state

The conformations of three 2',3'-difluoro uridine nucleosides were studied by X-ray crystallography, NMR spectroscopy, and ab initio calculations in an attempt to define the roles that the two vicinal fluorine atoms play in the puckering preferences of the furanose ring. Two of the compounds examined contained fluorine atoms in either the arabino or xylo dispositions at C2' and C3' of a 2',3'-dideoxyuridine system. The third compound also incorporated fluorine atoms in the xylo configuration on the furanose ring but was substituted with a 6-azauracil base in place of uracil. A battery of NMR experiments in D 2O solution was used to identify conformational preferences primarily from coupling constant and NOE data. Both (1)H and (19)F NMR data were used to ascertain the preferred sugar pucker of the furanose ring through the use of the program PSEUROT. Compound-dependent parameters used in the PSEUROT calculations were newly derived from complete sets of conformations calculated from high-level ab initio methods. The solution and theoretical data were compared to the conformations of each molecule in the solid state. It was shown that both gauche and antiperiplanar effects may be operative to maintain a pseudodiaxial arrangement of the C2' and C3' vicinal fluorine atoms. These data, along with previously reported data by us and others concerning monofluorinated nucleoside conformations, were used to propose a model of how fluorine influences different aspects of nucleoside conformations.     

Synthesis, 1 H, 13 C, 15 N, and 113 Cd NMR, ESI-TOF MS, Semiempirical MO (PM3), ab initio/HF and Cation/Anion Binding Studies of N-deoxycholyl-l-tryptophan

The synthesis, structural characterization, and cation/anion binding properties ofa new bile acid-amino acid conjugate, N-deoxycholyl-l-tryptophan, aredescribed. The structures of the ligand and its cadmium adduct at different pHconditions and various cadmium perchlorate concentrations were determined bymodern multinuclear magnetic resonance spectroscopic as well as ESI-TOF MStechniques. Also semiempirical PM3 and ab initio/HF molecular modellingstudies were performed. Based on ¹H,¹H NOESY measurementsN-deoxycholyl-l-tryptophan in alkaline conditions was found to appearin a bent conformation which was clearly different from the conformations in neutraland acidic solutions. According to molecular modelling in its minimum energy structurethe tryptophan backbone of the ligand was folded close to the deoxycholic acid skeletonand the structure was stabilized by an intramolecular hydrogen bond. The multinuclearmagnetic resonance experiments indicated that Cd²⁺-cation was bound with theligand in neutral and alkaline conditions while in acidic conditions protons block thebinding site. ESI-TOF MS revealed clearly a competition between sodium and cadmiumions, the ligand having a stronger affinity for sodium. Cadmium binding occurred onlywhen excess of cadmium was used. Further, ESI-TOF MS spectra showed that variouschlorine oxyanions originated from perchlorate anion formed together with cationsdifferent adducts with the ligand.     

Configuration and conformational equilibrium of (±)-trans-1-oxo-3-thiophen-2-yl-isochroman-4-carboxylic acid methyl ester

The X-ray analysis of 1-oxo-3-thiophen-2-yl-isochroman-4-carboxylic acid methyl ester 1 confirmed its trans-configuration and a conformation with diaxial H-3 and H-4 atoms in solid state. NMR experiments indicated that trans-1 exists in solution in both expected conformers. In CDCl₃ and especially in CD₃OD or DMSO, the conformational equilibrium is shifted towards the conformer with diequatorial H-3 and H-4, which was also determined by 2D NOESY experiments. The shift is due to the greater polarity of that conformer deduced by ab initio calculations.     

Conformational energy penalties of protein-bound ligands

The conformational energies required for ligands to adopt their bioactive conformations were calculated for 33 ligand–protein complexes including 28 different ligands. In order to monitor the force field dependence of the results, two force fields, MM3 and AMBER, were employed for the calculations. Conformational analyses were performed in vacuo and in aqueous solution by using the generalized Born/solvent accessible surface (GB/SA) solvation model. The protein-bound conformations were relaxed by using flat-bottomed Cartesian constraints. For about 70% of the ligand–protein complexes studied, the conformational energies of the bioactive conformations were calculated to be 3 kcal/mol. It is demonstrated that the aqueous conformational ensemble for the unbound ligand must be used as a reference state in this type of calculations. The calculations for the ligand–protein complexes with conformational energy penalties of the ligand calculated to be larger than 3 kcal/mol suffer from uncertainties in the interpretation of the experimental data or limitations of the computational methods. For example, in the case of long-chain flexible ligands (e.g. fatty acids), it is demonstrated that several conformations may be found which are very similar to the conformation determined by X-ray crystallography and which display significantly lower conformational energy penalties for binding than obtained by using the experimental conformation. For strongly polar molecules, e.g. amino acids, the results indicate that further developments of the force fields and of the dielectric continuum solvation model are required for reliable calculations on the conformational properties of this type of compounds.