Rotation in biphenyls with a single ortho-substituent

Barriers to rotation in a range up to 15.4 kcal mol(-1) were determined by dynamic NMR spectroscopy for a series of biphenyl compounds 1a-1h and 2a-2d with a single ortho-substituent. Ab initio calculations reproduce these barriers satisfactorily and indicate the ground-state conformation of these molecules. Results are discussed in terms of the contribution of individual substituents to the barrier and of the buttressing of adjacent positions in a benzene ring by substituents.     

Methyl rotation barriers in proteins from 2H relaxation data. Implications for protein structure

Side-chain 2H and backbone 15N relaxation data have been collected at multiple temperatures in the samples of the SH3 domain from alpha-spectrin. Combined analyses of the data allowed for determination of the temperature-dependent correlation times tauf characterizing fast methyl motion. Molecular dynamics simulations confirmed that tauf are dominated by methyl rotation; the corresponding activation energies approximate methyl rotation barriers. For 33 methyl groups in the alpha-spectrin SH3 domain the average barrier height was thus determined to be 2.8 +/- 0.9 kcal/mol. This value is deemed representative of the "fluid" hydrophobic protein core where some barriers are increased and others are lowered because of the contacts with surrounding atoms, but there is no local order that could produce systematically higher (lower) barriers. For comparison, the MD simulation predicts the average barrier of 3.1 kcal/mol (calculated via the potential of mean force) or 3.4-3.5 kcal/mol (rigid barriers after appropriate averaging over multiple MD snapshots). The latter result prompted us to investigate rigid methyl rotation barriers in a series of NMR structures from the Protein Databank. In most cases the barriers proved to be higher than expected, 4-6 kcal/mol. To a certain degree, this is caused by tight packing of the side chains in the NMR structures and stems from the structure calculation procedure where the coordinates are first annealed toward the temperature of 0 K and then subjected to energy minimization. In several cases the barriers >10 kcal/mol are indicative of van der Waals violations. The notable exceptions are (i) the structures solved using the GROMOS force field where tight methyl packing is avoided (3.0-3.6 kcal/mol) and (ii) the structure solved by means of the dynamic ensemble refinement method (Lindorff-Larsen, K.; Best, R. B.; DePristo, M. A.; Dobson, C. M.; Vendruscolo, M. Nature 2005, 433, 128) (3.5 kcal/mol). These results demonstrate that methyl rotation barriers, derived from the experiments that are traditionally associated with studies of protein dynamics, can be also used in the context of structural work. This is particularly interesting in view of the recent efforts to incorporate dynamics data in the process of protein structure elucidation.     

Barriers to rotation in methyl formate by dynamic NMR spectroscopy and barriers to 1,3 oxygen-to-oxygen migration in methyl formate and trifluoromethyl formate by ab initio calculations

Free-energy barriers of 9.85 and 11.91 +/- 0.15 kcal/mol at -70.8 degrees C were found by dynamic NMR spectroscopy for the E-to-Z and Z-to-E conversions, respectively, of methyl formate (1) enriched in 13C to 99% for the carbonyl carbon [methyl formate 13C (2)]. These barriers are higher than the literature values reported for -53 degrees C. The free-energy barrier to 1,3 oxygen-to-oxygen migration of the methyl group in methyl formate was determined by ab initio calculations at several levels. The value of 58.7 kcal/mol obtained at the MP2/6-311+G (df,pd) level was compared to a literature barrier for this process (MINDO/3) and to barriers for related compounds. A free-energy barrier of 63.0 kcal/mol for the oxygen - to - oxygen migration of the CF3 group in trifluoromethyl formate (3) was calculated at the MP2/6-31+G level.     

Conformational studies by dynamic NMR. 99. Experimental and computed determination of rotation barriers in the crystalline state: the case of naphthylphenylsulfoxide

The (13)C NMR CP-MAS spectrum of 2-naphthylphenylsulfoxide in the solid state displays line broadening effects due to the restricted rotation about the Ph-S bond. Line shape simulation of the temperature-dependent traces allowed the corresponding barrier to be determined in the solids (14.7 kcal mol(-1)). By making use of the information obtained from single-crystal X-ray diffraction, this barrier could be satisfactorily reproduced by theoretical calculations (14.5 kcal mol(-1)) that take into account the correlated phenyl motion involving a large set of molecules in the crystalline state     

Conformation, Inversion Barrier, and Solvent-Induced Conformational Shift in Exo- and Endo/Exo-Calix[4]arenes

Calixarenes 4a and 4b having hydroxyl groups in endo and exo positions and the ethanediyl-bridged exo-calixarene 5a were synthesized by a stepwise strategy. Single-crystal X-ray structures were obtained for 4a and for the exo-calixarene 3d, showing the molecules to exist in the 1,2-alternate conformation which is also found for 4a,b in solution. The inversion barriers of 4a and 4b (10.3 and 10.8 kcal mol(-1)) are similar to that determined for the endo-dihydroxycalixarene 12, indicating that the additional intramolecular hydrogen bond between the exo OH groups does not decrease the flexibility of the molecule. In CDCl(3) solution exo-calixarene 5a adopts a 1,2-alternate conformation with the methyl group at the bridge located in an axial position, while in DMSO-d(6) the conformation adopted is the partial cone. Similar solvent-induced conformational shifts were found for the exo-calixarenes 3b and 3d. MM3 calculations predict that the cone form is the lowest energy conformation of 4 and the exo-calixarenes 3 and 5. The calculations suggest that the conformational preferences of the methyl group at the bridge for either the axial or equatorial positions are in large part determined by the repulsive steric interactions with the hydroxyl groups. The inversion barrier of 4b is satisfactorily reproduced by calculations, which indicate that the rotation of the exo rings is less energetically demanding than the rotation of the endo rings.     

The quenching of isopropyl group rotation in van der Waals molecular solids

X-ray diffraction experiments are employed to determine the molecular and crystal structure of 3-isopropylchrysene. Based on this structure, electronic structure calculations are employed to calculate methyl group and isopropyl group rotational barriers in a central molecule of a ten-molecule cluster. The two slightly inequivalent methyl group barriers are found to be 12 and 15 kJ mol(-1) and the isopropyl group barrier is found to be about 240 kJ mol(-1), meaning that isopropyl group rotation is completely quenched in the solid state. For comparison, electronic structure calculations are also performed in the isolated molecule, determining both the structure and the rotational barriers, which are determined to be 15 kJ mol(-1) for both the isopropyl group and the two equivalent methyl groups. These calculations are compared with, and are consistent with, previously published NMR (1)H spin-lattice relaxation experiments where it was found that the barrier for methyl group rotation was 11+/-1 kJ mol(-1) and that the barrier for isopropyl group rotation was infinite on the solid state NMR time scale.     

Substituent effects on structural stability of formyl ketene and analysis of vibrational spectra of formyl haloketenes and formyl methylketene.

The conformational behavior and the structural stability of formyl fluoroketene, formyl chloroketene and formyl methylketene were investigated by utilizing quantum mechanical DFT calculations at B3LYP/6-31I + + G** and ab initio calculations at MP2/6-311 + + G** levels. The three molecules were predicted to have a planar s-cis<-->s-trans conformational equilibrium. From the calculations, the direction of the conformational equilibrium was found to be dependent on the nature of the substituting group. In formyl haloketenes, the cis conformation, where the C=O group eclipses the ketenic group, was expected to be of lower energy than the trans conformer. In the case of formyl methylketene the conformational stability was reversed and the trans form (the aldehydic hydrogen eclipsing the ketenic group) was calculated to be about 2 kcal mol(-1) lower in energy than the cis form. The calculated cis-trans energy barrier was found to be in the order: fluoride (15.3 kcal mol(-1)) > chloride (13.1 kcal mol(-1)) > methyl (11.7 kcal mol(-1). Full optimization was performed at the ground and the transition states of the molecules. The vibrational frequencies for the stable conformers of the three ketenic systems were computed at the DFT-B3LYP level, and the zero-point corrections were included into the calculated rotational barriers. Complete vibrational assignments were made on the basis of both normal coordinate calculations and comparison with experimental results of similar molecules.     

Stereodynamics of 9,11-Diphenyl-10-azatetracyclo[6.3.0.0.(4,11)0.(5,9)]undecanes. Highly Restricted Nitrogen Inversion and Isolated Phenyl Rotation. X-ray Crystallographic, Dynamic NMR, and Molecular Mechanics Studies

The (1)H NMR spectra of 10-benzyl-9,11-diphenyl-10-azatetracyclo[6.3.0.0.(4,11)0.(5,9)]undecane (BnPh(2)()) and 10-methyl-9,11-diphenyl-10-azatetracyclo[6.3.0.0.(4,11)0.(5,9)]undecane (MePh(2)()) decoalesce due to slowing inversion at nitrogen and to slowing isolated bridgehead phenyl rotation. The high nitrogen inversion barriers in MePh(2)() (DeltaG() = 12.2 +/- 0.1 kcal/mol at 250 K) and BnPh(2)() (DeltaG() = 10.6 +/- 0.1 kcal/mol at 215 K) are typical of tertiary amines in which at least one C-N-C bond angle is constrained to a small value. Compared to the minuscule rotation barriers about sp(2)-sp(3) carbon-carbon bonds in simple molecular systems, the bridgehead phenyl rotation barriers in MePh(2)() (DeltaG() = 9.8 +/- 0.1 kcal/mol at 210 K) and BnPh(2)() (DeltaG() = 9.8 +/- 0.1 kcal/mol at 210 K) are unusually high. Molecular mechanics calculations (MMX force field) suggest that the origin of the high phenyl rotation barriers lies in the close passage of an o-phenyl proton and a methyl (or benzylmethylene) proton in the transition state. BnPh(2)() crystallized from hexane as white needles in the monoclinic system Pn. Unit cell dimensions are as follows: a = 12.198(1) ?, b = 6.1399(6) ?, c = 14.938(2) ?, beta = 107.470(4) degrees, V = 1067.1(2) ?(3), Z = 2. In the crystal molecular structure, the imine bridge CNC bond angle in BnPh(2)() is constrained to a small value (96 degrees ). The benzylic phenyl group is oriented gauche to the nitrogen lone pair.     

Structure, Conformation, and Stereodynamics of N-Nitroso-2,4-diaryl-3-azabicyclo[3.3.1]nonanes and N-Nitroso-2,4-diaryl-3-azabicyclo[3.3.1]nonan-9-ones(1)

The variable temperature (1)H, (13)C, and (19)F NMR spectra were measured for the title N-nitrosamines. The observed unusually low N-N rotation barriers (12-15 kcal/mol) result from a significant deviation of the nitrosamino system from planarity. A pyramidal character of the amino nitrogen was confirmed by the X-ray crystal structures of two compounds and by bathochromic shifts of the n-pi absorption bands in the UV spectra. The nonplanarity of the nitrosamino moiety is due to the strong pseudoallylic A((1,3)) strain caused by the steric interaction of the NNO group with the neighboring aryl substituents fixed in the equatorial positions of the bicyclic skeleton. In addition, the barriers to the C-C rotation of aryl groups were examined at temperatures lower than required to "freeze" the N-N rotation and different DeltaG() values were observed for the aryls oriented syn and anti to the nitroso oxygen.     

Conformational studies by dynamic NMR. 97. Structure, conformation, stereodynamics and enantioseparation of aryl substituted norbornanes

The structure of a 1,7,7-triaryl norbornane (compound 3) has been determined by X-ray diffraction and was found essentially equal to that predicted by molecular mechanics calculations. Restricted rotation of the aryl groups also has been observed by dynamic NMR spectroscopy in this compound and in a number of analogously substituted norbornanes. The aryl-norbornane bond rotation barriers were measured by line shape analysis of the (13)C NMR spectra obtained at temperatures lower than -100 degrees C and were found to cover the range 6.0 to 7.9 kcal mol(-1). An exception was the rotation involving the o-anisyl group in compound 5, which occurs near ambient temperature since the corresponding barrier is much higher (14.4 kcal mol(-1)). In one case (compound 4) configurational enantiomers could be separated by chiral HPLC and the corresponding CD spectra recorded.     

Conformation and stereodynamics of symmetrically ortho-disubstituted aryl carbinols and aryl ethers

By making use of low-temperature dynamic NMR spectroscopy, the rotation barriers about the sp3-sp2 bond have been determined in a number of hindered benzyl alcohols symmetrically substituted in the ortho positions, the substituents being F, Cl, Br, and Me. The free energies of activation covered the range 4.6-10.1 kcal mol-1. Ab initio computations matched satisfactorily the trend of these values and predicted the conformation adopted by these compounds. In one case, this result could be also confirmed by the X-ray diffraction structure. In the case of the corresponding methyl ethers two barriers could be measured, corresponding to the passage across two distinguishable transition states: the higher barriers covered the range 5.0-8.1 kcal mol-1 and the lower ones the range 4.7-6.2 kcal mol-1.     

Importance of correlated motions on the low barrier rotational potentials of crystalline molecular gyroscopes

The energetic and structural changes taking place upon rotation of the central phenylene of 1,4-bis(3,3,3-triphenylpropynyl)benzene in the solid state were computed using molecular mechanics calculations. Pseudopolymorphic crystals of a benzene clathrate (1A) and a desolvated form (1B) were analyzed with models that account for varying degrees of freedom within the corresponding lattices. The calculated rotational barriers in a rigid lattice approximation, 78 kcal/mol for 1A and 72 kcal/mol for 1B, are about 5 times greater than those previously measured by variable-temperature 13C CPMAS NMR and quadrupolar echo 2H NMR line-shape analysis: 12.8 kcal/mol for 1A and 14.6 kcal/mol for 1B. The potential energy barriers calculated with a model that restricts whole body rotation and translational motions but allows for internal rotations give results that are near the experimental free-energy barriers. The calculated barriers for 1A and 1B are 15.5 and 16.2 kcal/mol, respectively. The differences between the rigid and partially relaxed models are attributed to the effect of correlated motions of the lattice and the rotating group, which are evident from the structural analysis of the atomic position data as a function of the dihedral angle of the rotator. The displacements of neighboring molecules near the rotary transition states for 1A and 1B can be as large as 2.7 and 1.1 A, respectively. The displacement and oscillation (C2) of interpenetrating phenyl rings from neighboring rotors proximal to the event are significant for both 1A and 1B. In addition, 6-fold (C6) benzene rotations in clathrate 1A were found to be directly correlated to the rotation of the phenylene rotator.     

The Curtius rearrangement of cyclopropyl and cyclopropenoyl azides. A combined theoretical and experimental mechanistic study

A combined experimental and theoretical study addresses the concertedness of the thermal Curtius rearrangement. The kinetics of the Curtius rearrangements of methyl 1-azidocarbonyl cycloprop-2-ene-1-carboxylate and methyl 1-azidocarbonyl cyclopropane-1-carboxylate were studied by (1)H NMR spectroscopy, and there is close agreement between calculated and experimental enthalpies and entropies of activation. Density functional theory (DFT) calculations (B3LYP/6-311+G(d,p)) on these same acyl azides suggest gas phase barriers of 27.8 and 25.1 kcal/mol. By comparison, gas phase activation barriers for the rearrangement of acetyl, pivaloyl, and phenyl azides are 27.6, 27.4, and 30.0 kcal/mol, respectively. The barrier for the concerted Curtius reaction of acetyl azide at the CCSD(T)/6-311+G(d,p) level exhibited a comparable activation energy of 26.3 kcal/mol. Intrinsic reaction coordinate (IRC) analyses suggest that all of the rearrangements occur by a concerted pathway with the concomitant loss of N2. The lower activation energy for the rearrangement of methyl 1-azidocarbonyl cycloprop-2-ene-1-carboxylate relative to methyl 1-azidocarbonyl cyclopropane-1-carboxylate was attributed to a weaker bond between the carbonyl carbon and the three-membered ring in the former compound. Calculations on the rearrangement of cycloprop-2-ene-1-oyl azides do not support pi-stabilization of the transition state by the cyclopropene double bond. A comparison of reaction pathways at the CBS-QB3 level for the Curtius rearrangement versus the loss of N2 to form a nitrene intermediate provides strong evidence that the concerted Curtius rearrangement is the dominant process.     

Synthesis,structure, and dynamic behavior of cyclopentadienyl-lithium, -sodium,and -potassium annelated with bicyclo[2.2.2]octene units: a systematic study on site exchange of alkali metals on a cyclopentadienyl ring in tetrahydrofuran

Novel cyclopentadienyl (Cp)-alkali metal complexes 1-M and 2-M (M = Li, Na, K), in which the Cp ring is annelated with two bicyclo[2.2.2]octene units and substituted with a phenyl group for 1 and a tert-butyl group for 2, were synthesized, and their structures and dynamic behaviors were investigated by means of X-ray crystallography, dynamic (13)C NMR, and DFT calculations. The X-ray crystallography results indicated that 1-Li, 1-Na, and 2-Na form monomeric contact ion pairs (CIP) with three THF molecules coordinated to the metal atom. Also, in THF-d(8), all of the 1-M and 2-M form monomeric CIP in the ground state. However, variable-temperature (13)C NMR measurements of 1-M and 2-M in THF-d(8) demonstrated dynamic behavior in which the metal ion exchanges positions between the upper and lower faces of the Cp ring. From a study of the concentration dependence of the dynamic behavior, the exchange was found to proceed principally as an intramolecular process at concentration ranges lower than 0.2 M. The experimentally observed deltaG values for the intramolecular exchange process for all the 1-M and 2-M (except for 2-Li, whose intramolecular process was too slow to observe) were found to be quite similar in THF-d(8) solution and to fall within the range of 12-14 kcal mol(-)(1). Within this range, a tendency was observed for the deltaG values to increase as the size of the metal decreased. Theoretical calculations (B3LYP/6-31G(d)) afforded considerably large values as the gas-phase dissociation energy for 1-M (162.7 kcal mol(-)(1) for M = Li; 131.6 kcal mol(-)(1) for M = Na; 110.9 kcal mol(-)(1) for M = K) and for 2-M (170.0 kcal mol(-)(1) for M = Li; 137.5 kcal mol(-)(1) for M = Na; 115.4 kcal mol(-)(1) for M = K). These values should be compensated for by a decrease in the solvation energies for the metal ions with increasing size, as exemplified by the calculated solvation energy for M(+)(Me(2)O)(4), which serves as a model for metal ions solvated with four molecules of THF (-122.9 kcal mol(-)(1) for M = Li; -94.7 kcal mol(-)(1) for M = Na; -67.7 kcal mol(-)(1) for M = K). This compensation results in a small difference in the overall energy for dissociation of 1-M or 2-M in ethereal solutions, thus supporting the similar deltaG values observed for the intramolecular metal exchange.     

Conformational studies by dynamic NMR. 74.(1) stereomutations of the conformational enantiomers in peri-substituted 1-acylnaphthalenes

Naphthalenes bearing an acyl and a phenyl group in a peri relationship give rise to a pair of enantiomers in the temperature range where the rotations of the acyl group are slow. Such enantiomers were observed by means of low temperature NMR spectra in chiral environments. The barrier to rotation for the acyl substituents, that causes the interconversion of the enantiomers, was demonstrated to be lower than that for the phenyl group. In an appropriately synthesized derivative it was possible to measure the two barriers that were found equal to 10.4 and 15.9 kcal mol(-)(1), respectively. The barriers for the acyl group rotation increase regularly (from 9.5 to 13.2 kcal mol(-)(1)) with the increasing dimension of the RCO groups (R = Me, Et, Pr(i), Bu(t)). When a bromine atom replaces the phenyl group, the enantiomerization barrier for the corresponding acyl derivatives increases significantly.     

NH tautomerism in the natural chlorin derivatives

The NH tautomerism of five Mg-free chlorophyll a and b derivatives 2-6 was studied utilizing NMR spectroscopy and molecular modeling. The results from the dynamic NMR measurements of the chlorins revealed that substituent effects contribute crucially to the free energy of activation (DeltaG(double dagger)) in the NH tautomeric processes. An intermediate tautomer for the total tautomeric NH exchange in a chlorin was observed for the first time, when the (1)H NMR spectra of chlorin e(6) TME (3) and rhodin g(7) TME (4) (TME = trimethyl ester) were measured at lower temperatures. The lower energy barriers (DeltaG(1)(double dagger)) obtained for the formation of the intermediate tautomers of 3 and 4, assigned to the N(22)-H, N(24)-H trans-tautomer, were 10.8 and 10.6 kcal/mol, respectively. The energy barrier (DeltaG(2)(double dagger) value) for the total tautomeric NH exchange in the five chlorins was found to vary from 13.6 kcal/mol to values higher than 18 kcal/mol. The lowest DeltaG(2)(double dagger) value (13.6 kcal/mol) was obtained for rhodochlorin XV dimethyl ester (2), which was the only chlorophyll derivative lacking the C(15) substituent. In the case of chlorins 4 and 5, the steric crowding around the methoxycarbonylmethyl group at C(15) raised the DeltaG(2)(double dagger) activation free-energy to 17.1 kcal/mol. However, the highest energy barrier with DeltaG(2)(double dagger) > 18 kcal/mol was observed for the NH exchange of pyropheophorbide a methyl ester (6), possessing the macrocycle rigidifying isocyclic ring E. Our results demonstrate that the steric strain, arising either from the steric crowding around the bulky substituent at C(15) or the macrocycle rigidifying isocyclic ring E, slows down the NH tautomeric process. We suggest that deformations in the chlorin skeleton are closely connected to the NH tautomeric exchange and that the exchange occurs by a stepwise proton-transfer mechanism via a hydrogen bridge.     

13C NMR, infrared, solvation and theoretical investigation of the conformational isomerism in 1-haloacetones (X = Cl, Br and I)

The solvent dependence of the 13C NMR spectra of chloroacetone (CA), bromoacetone (BA) and iodoacetone (IA) are reported and the 3J(CH) couplings analysed using ab initio calculations and solvation theory. In CA the energy difference (E(cis) - E(gauche)) between the cis (Cl-C-C=O 0 degrees) and gauche (Cl-C-C=O 155 degrees) conformers is 1.7 kcal mol(-1) in the vapour, decreasing to 0.8 kcal mol(-1) in CCl4 solution and to -1.0 kcal mol(-1) in the pure liquid. The conformational equilibrium, in BA, is between the more polar cis (Br-C-C=O 0 degrees) and gauche (Br-C-C=O 132 degrees) conformations. The energy difference (E(cis) - E(gauche)) is 1.8 kcal mol(-1) in the vapour, decreasing to 0.9 kcal mol(-1) in CCl4 solution and to -0.4 kcal mol(-1) in the pure liquid. The energy difference (E(cis) - E(gauche)), in IA, between the cis (I-C-C=O 0 degrees) and gauche (I-C-C=O 104 degrees) conformers is 1.1 kcal mol(-1) in the vapour phase, decreasing to 0.5 kcal mol(-1) in CCl4 solution and to -0.5 kcal mol(-1) in the pure liquid. The vapour state energy difference for BA [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)] and for IA [1.6 kcal mol(-1) at B3LYP/6-311++G(d,p)/LANL2DZ)] are in very good agreement with the above values. For CA the agreement is also satisfactory [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)].     

Probing the effects of heterogeneity on delocalized pi...pi interaction energies

Dimers composed of benzene (Bz), 1,3,5-triazine (Tz), cyanogen (Cy) and diacetylene (Di) are used to examine the effects of heterogeneity at the molecular level and at the cluster level on pi...pi stacking energies. The MP2 complete basis set (CBS) limits for the interaction energies (E(int)) of these model systems were determined with extrapolation techniques designed for correlation consistent basis sets. CCSD(T) calculations were used to correct for higher-order correlation effects (deltaE(CCSD)(T)(MP2)) which were as large as +2.81 kcal mol(-1). The introduction of nitrogen atoms into the parallel-slipped dimers of the aforementioned molecules causes significant changes to E(int). The CCSD(T)/CBS E(int) for Di-Cy is -2.47 kcal mol(-1) which is substantially larger than either Cy-Cy (-1.69 kcal mol(-1)) or Di-Di (-1.42 kcal mol(-1)). Similarly, the heteroaromatic Bz-Tz dimer has an E(int) of -3.75 kcal mol(-1) which is much larger than either Tz-Tz (-3.03 kcal mol(-1)) or Bz-Bz (-2.78 kcal mol(-1)). Symmetry-adapted perturbation theory calculations reveal a correlation between the electrostatic component of E(int) and the large increase in the interaction energy for the mixed dimers. However, all components (exchange, induction, dispersion) must be considered to rationalize the observed trend. Another significant conclusion of this work is that basis-set superposition error has a negligible impact on the popular deltaE(CCSD)(T)(MP2) correction, which indicates that counterpoise corrections are not necessary when computing higher-order correlation effects on E(int). Spin-component-scaled MP2 (SCS-MP2 and SCSN-MP2) calculations with a correlation-consistent triple-zeta basis set reproduce the trends in the interaction energies despite overestimating the CCSD(T)/CBS E(int) of Bz-Tz by 20-30%.     

Conformational preferences and internal rotation in alkyl- and phenyl-substituted thiourea derivatives

Potential energy surfaces (PES) for rotation about the N-C(sp(3)) or N-C(aryl) bond and energies of stationary points on PES for rotation about the C(sp(2))-N bond are reported for methylthiourea, ethylthiourea, isopropylthiourea, tert-butylthiourea, and phenylurea, using the MP2/aug-cc-pVDZ method. Analysis of alkylthioureas shows that conformations, with alkyl groups cis to the sulfur atom, are more stable (by 0.4-1.5 kcal/mol) than the trans forms. All minima adopt anti configurations with respect to nitrogen pyramidalization, whereas syn configurations are not stationary points on the MP2 potential surface. In contrast, analysis of phenylthiourea reveals that a trans isomer in a syn geometry is the global minimum, whereas a cis isomer in an anti geometry is a local minimum with a relative energy of 2.7 kcal/mol. Rotation about the C(sp(2))-N bond in alkyl and phenyl thioureas is slightly more hindered (9.1-10.2 kcal/mol) than the analogous motion in the unsubstituted molecule (8.6 kcal/mol). The maximum barriers to rotation for the methyl, ethyl, isopropyl, tert-butyl, and phenyl substituents are predicted to be 1.2, 8.9, 8.6, 5.3, and 0.9 kcal/mol, respectively. Corresponding PESs are consistent with the experimental dihedral angle distribution observed in crystal structures. The results of the electronic structure calculations are used to benchmark the performance of the MMFF94 force field. Systematic discrepancies between MMFF94 and MP2 results were improved by modification of selected torsion parameters and one of the van der Waals parameters for sulfur.     

Stereomutations of atropisomers of sterically hindered salophen ligands

[structure: see text] The stereomutations in nonsymmetrical salophen ligands 1-4 were studied by means of dynamic NMR and HPLC methods. DNMR experiments showed that in DMSO-d(6) hindered ligands 2-4 exist in two chiral conformations, depending on whether the imine carbon atoms are in a cis or trans disposition with respect to the plane of the central o-phenylenediamine ring, the latter being more stable by 1.0 kcal mol(-1). Owing to its higher dipole moment, in the apolar solvent C(6)D(6) the cis conformer is destabilized with respect to the trans one, in agreement with the results of ab initio calculations. In DMSO-d(6) solution the two conformers are in equilibrium through the less hindered rotation about the C6-N7 bond aligned to the a(6,7) axis, and the interconversion barriers range from 18.4 to 19.3 kcal mol(-1). The enantiomerization process is a two step-process that implies sequential rotations around the C6-N7 and the C1-N8 bonds, so that the rate determining step is the slower rotation around the more hindered C1-N8 bond aligned to the a(1,8) axis, and the energy barriers range from 21.4 to 21.9 kcal mol(-1). These values compare well with those determined by chromatography on an enantioselective HPLC column at low temperature, thus confirming that DNMR and DHPLC can be conveniently employed as complementary techniques.