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.     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

Force field parameters for sulfates and sulfamates based on ab initio calculations: Extensions of AMBER and CHARMm fields

Ab initio self-consistent field (SCF) Hartree-Fock calculations of sulfates R? O? SO₃(?1) (R = Me, Et, i-Pr) and sulfamates R? NHSO₃(?1) (R = H, Me, Et, i-Pr) were performed at the 4-31G(*S*N) //3-21G(*S*N) basis set levels, where asterisks indicate d functions on sulfur and nitrogen atoms. These standard levels were determined by comparing calculation results with several basis sets up to MP2/6-31G*//6-31G*. Several conformations per compound were studied to obtain molecular geometries, rotational barriers, and potential derived point charges. In methyl sulfate, the rotational barrier around the C? O bond is 1.6 kcal/mol at the MP2 level and 1.4 kcal/mol at the standard level. Its ground state has one of three HCOS torsion angles trans and one of three COSO torsion angles trans. Rotation over 60° around the single O? S bond in the sulfate group costs 2.5 kcal/mol at the MP2 and 2.1 kcal/mol at the standard level. For ethyl sulfate, the calculated rotational barrier in going from the ground state, which has its CCOS torsion angle trans, to the syn-periplanar conformation (CCOS torsion angle cis) is 4.8 kcal/mol. However, a much lower barrier of 0.7 kcal/mol leads to a secondary gauchelike conformation about 0.4 kcal/mol above the ground state, with the CCOS torsion angle at 87.6°. Again, one of the COSO torsion angles is trans in the ground state, and the rotational barrier for a 60° rotation of the sulfate group amounts to 1.8 kcal/mol. For methyl sulfamate, the rotational barriers are 2.5 kcal/mol around the C? N bond and 3.3 kcal/mol around the N? S bond. This is noteworthy because sulfamate itself has a calculated rotational barrier around the N? S bond of only 1.7 kcal/mol. These and other data were used to parameterize the well-known empirical force fields AMBER and CHARMm. When the new fields were tested by means of vibrational frequency calculations at the 6-31G*//6-31G* level for methyl sulfate, sulfamate, and methyl sulfamate ground states, the frequencies compared favorably with the AMBER and CHARMm calculated frequencies. The transferability of the force parameters to β-D -glucose-6-sulfate and isopropyl sulfate appears to be better than to isopropyl sulfamate. © 1995 by John Wiley & Sons, Inc.     

Synthesis of hydroxy and methoxy perylene quinones, their spectroscopic and computational characterization, and their antiviral activity

Hydroxy and methoxy perylene quinones are synthesized in an attempt to isolate the essential spectroscopic and biological features of light-induced antiviral agents such as hypericin and hypocrellin. Unlike their naturally occurring counterparts, these synthetic quinones bear the carbonyl, hydroxyl, and methoxy groups in the "bay region." The hydroxy and methoxy compounds have rich absorption spectra with broad features in the visible (approximately 450-800 nm) and relatively more intense and narrow features at wavelengths < or = 350 nm. High-level ab initio quantum mechanical calculations assign the features in the absorption spectra to electronic transitions from S0 to S2 and to higher-lying electronic states. The calculations indicate that in the ground state the trans dihydroxy isomer is 12.5 kcal/mol lower in energy than the cis dihydroxy isomer and is thus the only species present. The lowest-energy trans methoxy ground state isomer and the lowest-energy cis methoxy ground state isomer are found to be degenerate. An additional cis methoxy isomer 6.3 kcal/mol higher in energy than the global minimum is assumed to contribute to the spectrum and is also considered. Finally, the synthetic compounds exhibit similar light-induced antiviral activity to each other, but significantly less than that of hypericin.     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

Conformational analysis and rotational barriers of alkyl- and phenyl-substituted urea 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 methylurea, ethylurea, isopropylurea, tert-butylurea, and phenylurea, using the B3LYP/DZVP2 and MP2/aug-cc-pVDZ methods. The analysis of alkylureas reveals cis and (less stable) trans isomers that adopt anti geometries, whereas syn geometries do not correspond to stationary points. In contrast, the analysis of phenylurea reveals that the lowest energy form at the MP2 level is a trans isomer in a syn geometry. The fully optimized geometries are in good agreement with crystal structure data, and PESs are consistent with the experimental dihedral angle distribution. Rotation about the C(sp(2))-N bond in alkylureas and phenylurea is slightly more hindered (8.6-9.4 kcal/mol) than the analogous motion in the unsubstituted molecule (8.2 kcal/mol). At the MP2 level of theory, the maximum barriers to rotation for the methyl, ethyl, isopropyl, tert-butyl, and phenyl groups are predicted to be 0.9, 6.2, 6.0, 4.6, and 2.4 kcal/mol, respectively. The results are used to benchmark the performance of the MMFF94 force field. Systematic discrepancies between MMFF94 and MP2 results were improved by modification of several torsional parameters.     

The isomerization equilibrium between cis and trans chloride ruthenium olefin metathesis catalysts from quantum mechanics calculations

The cis-trans chloride isomerization of a ruthenium olefin metathesis catalyst is studied using quantum mechanics (B3LYP DFT), including the Poisson-Boltzmann (PBF) continuum approximation. The predicted geometries agree with experiment. The energies in methylene chloride, lead to DeltaG = -0.70 kcal/mol and a cis:trans ratio of 76:24, quite close to the experimental value of DeltaG = -0.78 kcal/mol or c:t 78:22. In contrast, we predict that in benzene c:t = 4:96 in agreement with the experimental observation of only the trans isomer. Our calculated relative activation energies explain the observed difference in initiation rates and suggest that each isomer should be isolable in high ratio by simply changing solvent.     

Energetic preferences for alpha,beta versus beta,gamma unsaturation

Density functional theory has been applied at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) level to examine the energetics of alpha,beta- versus beta,gamma-unsaturation for some common organic functional groups. Specifically, the relative stabilities of allyl-X (H2C=CHCH2X) and 1-propenyl-X (H3CCH=CHX) isomers have been computed for X = methyl, vinyl, phenyl, formyl, acetyl, methoxy, methylthio, methylsulfinyl, methylsulfonyl, sulfamoyl, and methoxysulfonyl, and the results are compared to available experimental data. The intrinsic preference of 3 kcal/mol for the 1-propenyl isomer when X = CH3 is exceeded by 2-4 kcal/mol for first-row conjugating groups, but it is not met for the sulfur-containing groups. In particular, alpha,beta-unsaturation is favored by less than 1 kcal/mol for the sulfone and sulfonamide analogues, while it is preferred by 8 kcal/mol for the vinyl-substituted case. Detailed structural results and torsional energy profiles are also reported.     

Intramolecular and intermolecular contributions to the barriers for rotation of methyl groups in crystalline solids: electronic structure calculations and solid-state NMR relaxation measurements

The rotation barriers for 10 different methyl groups in five methyl-substituted phenanthrenes and three methyl-substituted naphthalenes were determined by ab initio electronic structure calculations, both for the isolated molecules and for the central molecules in clusters containing 8-13 molecules. These clusters were constructed computationally using the carbon positions obtained from the crystal structures of the eight compounds and the hydrogen positions obtained from electronic structure calculations. The calculated methyl rotation barriers in the clusters (E(clust)) range from 0.6 to 3.4 kcal/mol. Solid-state (1)H NMR spin-lattice relaxation rate measurements on the polycrystalline solids gave experimental activation energies (E(NMR)) for methyl rotation in the range from 0.4 to 3.2 kcal/mol. The energy differences E(clust) - E(NMR) for each of the ten methyl groups range from -0.2 kcal/mol to +0.7 kcal/mol, with a mean value of +0.2 kcal/mol and a standard deviation of 0.3 kcal/mol. The differences between each of the computed barriers in the clusters (E(clust)) and the corresponding computed barriers in the isolated molecules (E(isol)) provide an estimate of the intermolecular contributions to the rotation barriers in the clusters. The values of E(clust) - E(isol) range from 0.0 to 1.0 kcal/mol.     

Thiaporphyrin with an Inverted Thiophene Ring – DFT Studies

The structure and electronic energy have been investigated applyingthe density functional theory (DFT) for two idealized 2-thia-21-carbaporphyrin, (SCP)H and (3H-SCP), 2-thia-21-carbaporphyrin anion (SCP)^- and 21-thiaporphyrin (SP)H. The analysis of calculated total electronic energies, using the B3LYP/6-31-G^* approach, demonstrates that the energy difference between 21-thiaporphyrin and its inverted isomer equals 10.34 kcal/mol [(SCP)H] and 52.08 kcal/mol [(3H-SCP)]. In contrast to 21-thiaporphyrin thethiophene fragments in (SCP)H and SCP^- present a geometry resembling that one of an isolated thiophene molecule.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Photophysical studies of the trans to cis isomerization of the push-pull molecule: 1-(pyridin-4-yl)-2-(N-methylpyrrol-2-yl)ethene (mepepy)

Organic molecules possessing intramolecular charge-transfer properties (D-pi-A type molecules) are of key interest particularly in the development of new optoelectronic materials as well as photoinduced magnetism. One such class of D-pi-A molecules that is of particular interest contains photoswitchable intramolecular charge-transfer states via a photoisomerizable pi-system linking the donor and acceptor groups. Here we report the photophysical and electronic properties of the trans to cis isomerization of 1-(pyridin-4-yl)-2-(N-methylpyrrol-2-yl)ethene ligand (mepepy) in aqueous solution using photoacoustic calorimetry (PAC) and theoretical methods. Density functional theory (DFT) calculations demonstrate a global energy difference between cis and trans isomers of mepepy to be 8 kcal mol(-1), while a slightly lower energy is observed between the local minima for the trans and cis isomers (7 kcal mol(-1)). Interestingly, the trans isomer appears to exhibit two ground-state minima separated by an energy barrier of approximately 9 kcal mol(-1). Results from the PAC studies indicate that the trans to cis isomerization results in a negligible volume change (0.9 +/- 0.4 mL mol(-1)) and an enthalpy change of 18 +/- 3 kcal mol(-1). The fact that the acoustic waves associated with the trans to cis transition of mepepy overlap in frequency with those of a calorimetric reference implies that the conformational transition occurs faster than the approximately 50 ns response time of the acoustic detector. Comparison of the experimental results with theoretical studies provide evidence for a mechanism in which the trans to cis isomerization of mepepy results in the loss of a hydrogen bond between a water molecule and the pyridine ring of mepepy.     

Conformational stability, vibrational assignmenents, barriers to internal rotations and ab initio calculations of 2-aminophenol (d 0 and d3)

The Raman (3700-100 cm(-1)) and infrared (4000-400 cm(-1)) spectra of solid 2-aminophenol (2AP) have been recorded. The internal rotation of both OH and NH2 moieties produce ten conformers with either Cs or C1 symmetry. However, the calculated energies as well as the imaginary vibrational frequencies reduce rotational isomerism to five isomers. The molecular geometry has been optimized without any constraints using RHF, MP2 and B3LYP levels of theory at 6-31G(d), 6-311+G(d) and 6-31++G(d,p) basis sets. All calculations predict 1 (cis; OH is directed towards NH2) to be the most stable conformation except RHF/6-31++G(d,p) basis set. The 1 (cis) isomer is found to be more stable than 8 (trans; OH is away from the NH2 moiety and the NH bonds are out-of-plane) by 1.7 kcal/mol (598 cm(-1)) as obtained from MP2/6-31G(d) calculations. Aided by experimental and theoretical vibrational spectra, cis and trans 2AP are coexist in solution but cis isomer is more likely present in the crystalline state. Aided by MP2 and B3LYP frequency calculations, molecular force fields, simulated vibrational spectra utilizing 6-31G(d) basis set as well as normal coordinate analysis, complete vibrational assignments for HOC6H4NH2 and DOC6H4ND2 have been proposed. Furthermore, we carried out potential surface scan, to determine the barriers to internal rotations of NH2 and OH groups. All results are reported herein and compared with similar molecules when appropriate.     

Proton migration and tautomerism in protonated triglycine.

Proton migration in protonated glycylglycylglycine (GGG) has been investigated by using density functional theory at the B3LYP/6-31++G(d,p) level of theory. On the protonated GGG energy hypersurface 19 critical points have been characterized, 11 as minima and 8 as first-order saddle points. Transition state structures for interconversion between eight of these minima are reported, starting from a structure in which there is protonation at the amino nitrogen of the N-terminal glycyl residue following the migration of the proton until there is fragmentation into protonated 2-aminomethyl-5-oxazolone (the b(2) ion) and glycine. Individual free energy barriers are small, ranging from 4.3 to 18.1 kcal mol(-)(1). The most favorable site of protonation on GGG is the carbonyl oxygen of the N-terminal residue. This isomer is stabilized by a hydrogen bond of the type O-H.N with the N-terminal nitrogen atom, resulting in a compact five-membered ring. Another oxygen-protonated isomer with hydrogen bonding of the type O-H.O, resulting in a seven-membered ring, is only 0.1 kcal mol(-)(1) higher in free energy. Protonation on the N-terminal nitrogen atom produces an isomer that is about 1 kcal mol(-)(1) higher in free energy than isomers resulting from protonation on the carbonyl oxygen of the N-terminal residue. The calculated energy barrier to generate the b(2) ion from protonated GGG is 32.5 kcal mol(-)(1) via TS(6-->7). The calculated basicity and proton affinity of GGG from our results are 216.3 and 223.8 kcal mol(-)(1), respectively. These values are 3-4 kcal mol(-)(1) lower than those from previous calculations and are in excellent agreement with recently revised experimental values.     

10]Annulene: bond shifting and conformational mechanisms for automerization

We report density-functional and coupled-cluster calculations on conformation change and degenerate bond shifting in [10]annulene isomers 1-5. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, conversion of the twist (1) to the heart (2) has a barrier of 10.1 kcal/mol, compared to Ea = 16.2 kcal/mol for degenerate "two-twist" bond shifting in 1. Pseudorotation in the all-cis boat isomer (3) proceeds with a negligible barrier. The naphthalene-like isomer 4 has a 3.9 kcal/mol barrier to degenerate bond shifting. The azulene-like isomer 5 is the only species for which the nature of the bond-equalized form (5-eq) depends on the method. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, 5-eq is 1.2 kcal/mol more stable than the bond-alternating form 5-alt. Conversion of 5-eq to 4 has a barrier of 12.6 kcal/mol. Despite being significantly nonplanar, both 5-eq and the transition state for bond shifting in 4 are highly aromatic based on magnetic susceptibility exaltations. On the basis of a detailed consideration of these mechanisms and barriers, we can now, with greater confidence, rule out 4 and 5 as candidates to explain the NMR spectra observed by Masamune. Our results support Masamune's original assignments for both isolated isomers.     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

Kinetic Study of 1-β-Cyanoethyl Aziridine for Ring-Opening Polymerization Initiated with Alkyl Sulfate

The kinetics of the cationic polymerization of 1-β-cyanoethyl aziridine initiated 3-hydroxy- 1-propane sulfonic acid sultone and methyl tosylate have been studied. The course of polymerization involved the propagation stage and termination reaction due to the reaction between the growing chain and imino groups in the polymer chain. The propagation constants and termination constants were obtained. The enthalpies of activation for the propagation and termination reactions are ΔH_p? = 12.9 kcal/mol and ΔH_t? = 12.4 kcal/mol, and the entropies of activation are ΔS_p? = -31 cal/deg·mol and ΔS_t? = -39 cal/deg·mol. Otherwise, the polymerization initiated with methyl tosylate involved an early stage which was initiated very quickly.     

Electronic Factors Influencing the Decarboxylation of beta-Keto Acids. A Model Enzyme Study

A theoretical study of the mechanism of decarboxylation of beta-keto acids is described. A cyclic transition structure was found with essentially complete proton transfer from the carboxylic acid to the beta-carbonyl group. The activation barrier for decarboxylation of formylacetic acid is predicted to be 28.6 kcal/mol (MP4SDTQ/6-31G//MP2/6-31G) while loss of CO(2) from its anion exhibits a barrier of only 20.6 kcal/mol (MP4SDTQ/6-31+G//MP2/6-31+G). Barrier heights of decarboxylation of malonic acid and alpha,alpha-dimethylacetoacetic acid are predicted to be 33.2 and 26.7 kcal/mol, respectively. Model enzyme studies using a thio methyl ester of malonate anion suggests that the role of malonyl-CoA is to afford a polarizable sulfur atom to stabilize the developing enolate anion in the transition structure for decarboxylation. Adjacent positively charged ammonium ions are also observed to stabilize the loss of CO(2) from a carboxylate anion by through-bond Coulombic stabilization of the transition structure.     

Conformational selectivity in the diels-alder cycloaddition: predictions for reactions of s-trans-1,3-butadiene

Diels-Alder cycloaddition of s-trans-1,3-butadiene (1) should yield trans-cyclohexene (7), just as reaction of the s-cis conformer gives cis-cyclohexene (9). Investigation of this long-overlooked process with Hartree-Fock, Moller-Plesset, CASSCF, and DFT methods yielded in every case a C(2)-symmetric concerted transition state. At the B3LYP/6-31G (+ZPVE) level, this structure is predicted to be 42.6 kcal/mol above reactants, while the overall reaction is endothermic by 16.7 kcal/mol. A stepwise diradical process has been studied by UBLYP/6-31G theory and found to have barriers of 35.5 and 17.7 kcal/mol for the two steps. Spin correction lowers these values to 30.1 and 13.0 kcal/mol. The barrier to pi-bond rotation in cis-cyclohexene (9) is predicted (B3LYP theory) to be 62.4 kcal/mol, with trans-cyclohexene (7) lying 53.3 kcal/mol above cis isomer 9. Results suggest that pi-bond isomerization and concerted reaction may provide competitive routes for Diels-Alder cycloreversion. It is concluded that full understanding of the Diels-Alder reaction requires consideration of both conformers of 1,3-butadiene.