Biomedical subjects
Publications and source records attributed to Kuo-Hsiang Chen.
Cycloundecane (1) was shown to exist at -183.1 degrees C as a mixture of the [12323] (approximately 59%) and [335] (approximately 41%) conformations. Populations were determined from the (13)C NMR spectrum, and assignments were based on the (13)C spectra, calculated free energies and chemical shifts, and information from the literature, including X-ray studies of solid derivatives and calculated barriers.
A molecular mechanics study of small saturated hydrocarbons (up to C-6) substituted by up to six fluorines has been carried out with the MM4 force field. A parameter set has been developed for use in the calculation of bond lengths, bond angles, torsion angles, conformational energies, barriers to rotation, dipole moments, moments of inertia, and vibrational frequencies for these compounds. The results are mostly in fair to good agreement with experiment and ab initio calculations. The high electronegativity of fluorine leads to serious geometric consequences in these compounds, but these consequences can be dealt with adequately by suitable cross-terms in the force constant matrix, and by recognizing that some of the reference bond lengths and angles (l(0), theta(0)) and the corresponding stretching and bending constant parameters (k(s), k(theta)) that are usually thought of as constants must in fact be treated as functions of the electronegativity of the substituents. Additionally, the heavy mass of the fluorine (relative to the mass of hydrogen in alkanes) leads to large values for other cross-terms that were found to be unimportant in hydrocarbons. Conformational equilibria for polyfluorinated compounds are affected by the delta-two effect well-known in carbohydrates. A few larger fluorinated and polyfluorinated alkanes, including perfluoropropane, perfluorobutane, and Teflon, have also been studied.
The rotational barrier for a methyl group at the end of an anomeric system is sometimes lower than we might have anticipated. Thus, in the trans-trans conformation of dimethoxymethane, the barrier to methyl rotation is calculated (B3LYP/6-311++G(2d,2p)) to be 2.22 kcal/mol, just slightly smaller than the corresponding barrier to rotation of the methyl group in methyl propyl ether of 2.32 kcal/mol. However, if the methyl being rotated in dimethoxymethane is placed into a gauche conformation, that rotational barrier is reduced to 1.52 kcal/mol. This substantial (0.80 kcal/mol relative to methyl propyl ether) reduction in barrier height in the latter case is attributed mainly to the change in the bond order of the C-O bond to which the methyl is attached, as a function of conformation, which in turn is a result of the anomeric effect. We have called this barrier lowering the external-anomeric torsional effect. This effect is apparently widespread in carbohydrates, and it results in the changing of conformational energies by up to about 2 kcal/mol. If polysaccharide potential surfaces are to be accurately mapped by molecular mechanics, this effect clearly needs to be accounted for.
The MM3 force field has been extended to deal with the lithium amide molecules that are widely used as efficient catalysts for stereoselective asymmetric synthesis. The MM3 force field parameters have been determined on the basis of the ab initio MP2/6-31G* and/or DFT (B3LYP/6-31G*, B3-PW91/6-31G*) geometry optimization calculations. To evaluate the electronic interactions specific to the lithium amides derived from the diamine molecules properly, the Lewis bonding potential term for the interaction between the lithium atom and the nonbonded adjacent electronegative atom such as nitrogen was introduced into the MM3 force field. The bond dipoles were evaluated correctly from the electronic charges on the atoms calculated by fitting to the electrostatic potential at points selected. The MM3 results on the molecular structures, conformational energies, and vibrational spectra show good agreement with those from the quantum mechanical calculations.
Simple alcohols and ethers have been studied with the MM4 force field. The structures of 13 molecules have been well fit using the MM4 force field. Moments of inertia have been fit with rms percentage errors as indicated: 18 moments for ethers, 0.28%; 21 moments for alcohols, 0.22%. Rotational barriers and conformational equilibria have also been examined, and the experimental and ab initio results are reproduced substantially better with MM4 than they were with MM3. Much of the improvement comes from the use of additional interaction terms in the force constant matrix, of which the torsion-bend and torsion-torsion are particularly important. Induced dipoles are included in the calculation, and dipole moments are reasonably well fit. It has been possible for the first time to fit conformational energetic data for both open chain and cyclic alcohols (e.g., propanol and cyclohexanol) with the same parameter set. For vibrational spectra, over a total of 82 frequencies, the rms error is 27 cm(-1), as opposed to 38 cm(-1) with MM3. Both the alpha and beta bond shortening resulting from the presence of the electronegative oxygen atom in the molecule are well reproduced. The electronegativity of the oxygen is sufficient that one must also include not only the alpha and beta electronegativity effects on bond lengths, but also on angle distortions, if structures are to be well reproduced. The heats of formation of 32 alcohols and ethers were fit overall to within experimental error (weighted standard deviation error 0.26 kcal/mol).
The anomeric effect has been studied for a variety of compounds using the MM4 force field, and also using MP2/6-311++G(2d,2p) ab initio calculations and experimental data for reference purposes. Geometries and energies, including conformational, rotational barriers, and heats of formation were examined. Overall, the agreement of MM4 with the experimental and ab initio data is good, and significantly better than the agreement obtained with the MM3 force field. The anomeric effect is represented in MM4 by various explicit terms in the force constant matrix. The bond length changes are accounted for with torsion-stretch elements. The angle changes are accounted for with torsion-bend elements. The energies are taken into account with a number of torsional terms in the usual way. A torsion-torsion interaction is also of some importance. With all of these elements included in the calculation, the MM4 results now appear to be adequately accurate. The heats of formation were examined for a total of 12 anomeric compounds, and the experimental values were fit by MM4 with an RMS error of 0.42 kcal/mol.
Ethylene glycol, its dimethyl ether, and some related compounds have been studied using the MM4 molecular mechanics force field. The MM4 calculated structural and energetic results have been brought into satisfactory agreement with a considerable number of experimental data and MP2/6-311++G(2d,2p) ab initio calculations. The heats of formation of these compounds are also well calculated. The MM4 ethylene glycol conformations in particular are in good agreement, both geometrically and in terms of energy, with those from the ab initio calculations. The corresponding dimethyl ether is of special interest, because it has been suggested that the trans-gauche conformation is unusually stable due to the hydrogen bonding of a hydrogen on a methyl group with the more distant oxygen. It is shown in the present work that while this conformation is more stable than might have been expected, the energy is adequately calculated by MM4 without using any hydrogen bonding between the Cbond;H bond and the oxygen. If such hydrogen bonding occurs, it amounts to no more than about 0.5 kcal/mol in energy, and is too small to detect with certainty. Additionally, energetic relationships in trans-1,2-dimethoxycyclohexane, 1,3,5,7-tetraoxadecalin, and 3-methoxytetrahydropyran have been studied, and the calculated results are compared with experimental information, which is adequately reproduced.
Ab initio calculations [B3LYP/6-311++G(2d,2p)] have been carried out on 84 conformations of 12 different sugars (hexoses), in both pyranose and furanose forms, with the idea of generating a data base for carbohydrate structural energies that may be used for developing the predictive value of molecular mechanics calculations for carbohydrates. The average value for the apparent gas phase anomeric effect for a series of 31 pairs of pyranose conformations was found to be 1.83 kcal/mol (vs. 2.67 kcal/mol with a smaller basis set used in earlier calculations). In developing MM4 to reproduce these data, it was necessary first to have good energies for simple alcohols and ethers, together with an adequate treatment of hydrogen bonding, and then to include the anomeric effect, and the ethylene glycol type system, as was previously recognized. It was also found that the so-called delta-2 effect, long recognized in carbohydrates, must be explicitly included, in order to obtain acceptable results. When a force field that included all of these items as developed from the small molecules based on the MM4 hydrocarbon force field was applied without any parameter adjustment to the set of hexopyranose and furanose conformations mentioned earlier, the E(beta) - E(alpha) was found to have an average value of 1.88 kcal/mol, versus 1.74 for the quantum calculations. The signed average and RMS deviations of the MM4 from the QM results were +0.15 and 0.87 kcal/mol.
Neutral in-plane tris-homoaromaticity is evaluated in tris(bismethano)benzene (15) and modifications of this parent structure in which the pi-orbitals might interact in the plane established by the unsaturated carbon atoms (in-plane conjugation). On the basis of magnetic susceptibility exaltation, nucleus-independent shift (NICS), and aromatic stabilization energy (ASE, evaluated via homodesmotic and isodesmic equations using B3LYP/6-311+G + ZPVE energies, as well as by MM3 and MM4 force field computations), we identified triene 17, a triply bridged analogue of 15, as the system where homoaromaticity is most effective. The NICS(total) in the center of 17 is -30.1 ppm and the diatropic pi-contribution is -18.0 ppm. This structure possesses more than one-third of the aromatic stabilization of benzene and is the best candidate for neutral tris-homoaromaticity ever proposed. The previously described tris-(bismethano)-benzene (15) also shows homoaromaticity but to a smaller extent compared to 17. Structure 18, which is closely related to 17, also is significantly homoaromatic but, as evaluated by MM3, strain partially counteracts the stabilizing effects from homoconjugation. Such a counteracting increase in strain largely cancels or even overwhelms the stabilization from homoconjugation in all other species considered in this study.