Calculating accurate proton chemical shifts of organic molecules with density functional methods and modest basis sets

The purpose of this paper is to convince practitioners of (1)H NMR spectroscopy to consider simple quantum chemical calculations as a viable option to aid them in the assignment of their spectra. To this end, it is demonstrated, on a test set of 80 conformationally stable molecules of various kinds carrying different functional groups, that, in contrast to what is claimed in the literature, large basis sets are not needed to obtain rather accurate predictions of (1)H NMR chemical shifts by quantum chemical calculations. On the other hand, modeling the solvent by an SCRF-type calculation may improve certain predictions significantly. The best accuracy/cost ratio is provided by GIAO calculations in chloroform as a solvent with the specially parametrized WP04 functional of Cramer et al. using the cc-pVDZ or 6-31G** basis set, closely followed by similar calculations with the ubiquitious B3LYP functional (both predict (1)H chemical shifts with an average deviation of ca. 0.12 ppm, if the results are scaled linearly). A slightly higher accuracy can be attained by adding diffuse functions to the basis set, but going to the triple-zeta basis sets which have invariably been used hitherto in calculations of chemical shifts does not lead to any improvement. The popular increment schemes such as those implemented in the ChemDraw or ACD programs do not do nearly as well and are often incapable of correctly distinguishing stereoisomers.

Optimized end-stacking provides specificity of N-methyl mesoporphyrin IX for human telomeric G-quadruplex DNA

N-methyl mesoporphyrin IX (NMM) is exceptionally selective for G-quadruplexes (GQ) relative to duplex DNA and, as such, has found a wide range of applications in biology and chemistry. In addition, NMM is selective for parallel versus antiparallel GQ folds, as was recently demonstrated in our laboratory. Here, we present the X-ray crystal structure of a complex between NMM and human telomeric DNA dAGGG(TTAGGG)(3), Tel22, determined in two space groups, P2(1)2(1)2 and P6, at 1.65 and 2.15 Å resolution, respectively. The former is the highest resolution structure of the human telomeric GQ DNA reported to date. The biological unit contains a Tel22 dimer of 5'-5' stacked parallel-stranded quadruplexes capped on both ends with NMM, supporting the spectroscopically determined 1:1 stoichiometry. NMM is capable of adjusting its macrocycle geometry to closely match that of the terminal G-tetrad required for efficient π-π stacking. The out-of-plane N-methyl group of NMM fits perfectly into the center of the parallel GQ core where it aligns with potassium ions. In contrast, the interaction of the N-methyl group with duplex DNA or antiparallel GQ would lead to steric clashes that prevent NMM from binding to these structures, thus explaining its unique selectivity. On the basis of the biochemical data, binding of NMM to Tel22 does not rely on relatively nonspecific electrostatic interactions, which characterize most canonical GQ ligands, but rather it is hydrophobic in nature. The structural features observed in the NMM-Tel22 complex described here will serve as guidelines for developing new quadruplex ligands that have excellent affinity and precisely defined selectivity.

Spontaneous transfer of chirality in an atropisomerically enriched two-axis system

One of the most well-recognized stereogenic elements in a chiral molecule is an sp(3)-hybridized carbon atom that is connected to four different substituents. Axes of chirality can also exist about bonds with hindered barriers of rotation; molecules containing such axes are known as atropisomers. Understanding the dynamics of these systems can be useful, for example, in the design of single-atropisomer drugs or molecular switches and motors. For molecules that exhibit a single axis of chirality, rotation about that axis leads to racemization as the system reaches equilibrium. Here we report a two-axis system for which an enantioselective reaction produces four stereoisomers (two enantiomeric pairs): following a catalytic asymmetric transformation, we observe a kinetically controlled product distribution that is perturbed from the system's equilibrium position. As the system undergoes isomerization, one of the diastereomeric pairs drifts spontaneously to a higher enantiomeric ratio. In a compensatory manner, the enantiomeric ratio of the other diastereomeric pair decreases. These observations are made for a class of unsymmetrical amides that exhibits two asymmetric axes--one axis is defined through a benzamide substructure, and the other axis is associated with differentially N,N-disubstituted amides. The stereodynamics of these substrates provides an opportunity to observe a curious interplay of kinetics and thermodynamics intrinsic to a system of stereoisomers that is constrained to a situation of partial equilibrium.

Atomic Charges

The problem of deriving atomic charges from the results of ab Initio MO calculations has been studied by the use of several reported methods: Mulliken population analysis, the minimal basis set (MBS) procedure, the natural population analysis (NPA), two electrostatic potential fitting methods, M-K and ChelpG, the Hirshfeld population analysis, and CM5 (charge model 5), which is related to the Hirshfeld method. The first set of studies were concerned with hydrogen charges. It was found that the MBS charges were linearly related to the Hirshfeld charges. The Hirshfeld, CM5, and MBS methods, but not the others, provided an excellent correlation for H atomic charge with the H-C-H bond angle, and with calculated gas-phase acidity. The two methods that were linearly related and gave hydrogen charges in agreement with an experimental study of partially deuterated methanes are MBS and Hirshfeld. In order to see which of the two methods is the more satisfactory, the methanol dimer was examined. The calculated H bond energy was 6.2 kcal/mol, which was in good agreement with studies of hydrogen bonds. The Coulombic interaction for the O···H bond was estimated using the MBS and Hirshfeld charges. The latter gave a calculated energy of 3-6 kcal/mol, whereas MBS gave an energy of ∼35 kcal/mol. Clearly, the Hirshfeld method is more satisfactory and should be the method of choice.

N-methylimidazole promotes the reaction of homophthalic anhydride with imines

The addition of N-methylimidazole (NMI) to the reaction of homophthalic anhydride with imines such as pyridine-3-carboxaldehyde-N-trifluoroethylimine (9) reduces the amount of elimination byproduct and improves the yield of the formal cycloadduct, tetrahydroisoquinolonic carboxylate 10. Carboxanilides of such compounds are of interest as potential antimalarial agents. A mechanism that rationalizes the role of NMI is proposed, and a gram-scale procedure for the synthesis and resolution of 10 is also described.

Are the enolates of amides and esters stabilized by electrostatics?

The fact that amides and esters form less stable enolates than ketones might be seen as evidence that electrostatic stabilization is unimportant in these anions. However, ab initio molecular orbital calculations show that electrostatic stabilization does in fact lie beneath the competing resonance effect that causes the decrease in acidity. The electrostatic contribution is revealed by examining torsionally twisted amide and ester structures in which the pi resonance interactions are largely inhibited. These twisted amides and esters have greater enolate acidity than the corresponding ketones. Qualitatively similar behavior is observed with respect to protonation, such that twisted amides and esters are generally less basic than the reference ketones, in striking contrast to their behavior in the normal geometries.