Ab initio calculations on Pro-Ala and Pro-Gly dipeptides

The geometries of the dipeptides -Pro- -Ala, -Pro- -Ala and -Pro-Gly were investigated by a grid scan ab initio calculation. The 6-31G basis set was used to estimate the effect of the alanyl side-chain on the conformation of the peptide backbone and to provide a computational basis for the interpretation of known physical-chemical properties of larger peptides that contain these dipeptides. These calculations furnish a direct quantum mechanical assessment of the energetic consequences of a methyl side-chain in the i + 2 position of a turn. The results of the calculation support the current view that the presence of a -Ala residue in the i + 2 position favors a type II β-turn over a type I β-turn conformation, while -Ala has the opposite effect. Total and relative energies for all the optimized conformations identified by the grid search are given and geometric parameters (bond lengths, bond angles and dihedral angles) and net atomic charges have been calculated.     

Molecular mechanics (MM3) parameterization for oxocarbenium ions

The physical properties of a diverse group of 12 oxocarbenium ions have been studied with ab initio calculations at the MP2/6‐31+G* level of theory. Based on theoretically derived properties such as molecular equilibrium geometry, dipole moment, and vibrational frequencies, a molecular mechanics (MM3) force field has been developed with the assistance of the programs TORSMART and MPMSR, components of our artificial parameter development and refinement method. The MM3 force field is now able to reproduce bond lengths, bond angles, moments of inertia, dipole moments, torsional energy profiles, and vibrational frequencies of oxocarbenium ions, which will allow further studies of glycoside hydrolysis and their rates of reaction. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 329–339, 2000     

Molecular simulation of guanidinium-based ionic liquids

A new systematic all-atom force field was developed for cyclic guanidinium-based ionic liquids (ILs) based on the AMBER force field. Optimized molecular geometries and equilibrium bond lengths and angles were obtained by ab initio calculations, and charges were allocated to each atom center by fitting the ab initio electrostatic potential. Molecular dynamics simulations were performed for eleven kinds of ILs that are comprised of NO3(-) anions and cyclic guanidinium-based cations. Validation was carried out by comparing our simulated densities with experimental and calculated data from the literature. Transport properties such as self-diffusion coefficients, viscosities, and conductivities were calculated by molecular dynamic simulation, and their dependence on the length of the alkyl chains of cyclic guanidinium-based cations are discussed. Radial distribution functions and spatial distribution functions were investigated to depict the microscopic structures of the ILs, and the relationship between their properties and microstructures is also discussed.     

Study of peptide conformation in terms of the ABEEM/MM method

The ABEEM/MM model (atom-bond electronegativity equalization method fused into molecular mechanics) is applied to study of the polypeptide conformations. The Lennard-Jones and torsional parameters were optimized to be consistent with the ABEEM/MM fluctuating charge electrostatic potential. The hydrogen bond was specially treated with an electrostatic fitting function. Molecular dipole moments, dimerization energies, and hydrogen bond lengths of complexes are reasonably achieved by our model, compared to ab initio results. The ABEEM/MM fluctuating charge model reproduces both the peptide conformational energies and structures with satisfactory accuracy with low computer cost. The transferability is tested by applying the parameters of our model to the tetrapeptide of alanine and another four dipeptides. The overall RMS deviations in conformational energies and key dihedral angles for four di- or tetrapeptide, is 0.39 kcal/mol and 7.7 degrees . The current results agree well with those by the accurate ab initio method, and are comparable to those from the best existing force fields. The results make us believe that our fluctuating charge model can obtain more promising results in protein and macromolecular modeling with good accuracy but less computer cost.     

Theoretical studies on the structure and electronic properties of aryl sulfides and sulfones

Molecular orbital calculations are reported on the structure and electronic properties of diphenyl sulfide using both semiempirical and ab initio methods. Neither the MNDO nor AM1 methods give satisfactory structures, but better results are obtained with the PM3 method. At the ab initio level, the 4-31G basis set with polarization functions on sulfur alone (4-31G/S*) gives comparable results to those obtained with the 6-31G** basis set. The corresponding bond lengths and angles at the sulfur atom of 4-aminophenyl-4′-nitrophenyl sulfide and related derivatives of diphenyl sulfone, diphenyl disulfide, and phenylthiosulfonate calculated at the 4-31 G/S* level show a good correlation with crystallographic data where available. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 141–147, 1998     

The gas-phase molecular structure of 1,1-diethynylsilacyclobutane as determined by means of electron diffraction and ab initio calculations

As a continuation of our systematic investigation of the effect of substituents on the ring geometry and dynamics in silacyclobutanes and in order to explore the role of the silicon atom as a mediator for electronic interactions between the attached fragments, we studied the molecular structure of 1,1-diethynylsilacyclobutane (DESCB) by means of gas-phase electron diffraction and ab initio calculations. The structural refinement of the electron diffraction data yielded the following bond lengths (r_a) and bond angles (uncertainties are 3σ): r(Si–C)=1.874(2) Å, r(Si–C)=1.817(1) Å, (C–Si–C)=79.2(6)°, (C–Si–C)=106.5(6)°. The geminal Si–CC moieties were found to be bent outwards by 3.1(15)° and the puckering angle was determined to be 30.0(15)°. The evidently short Si–C bond length, which was also reproduced by the ab initio calculations, could be rationalized as being the consequence of the electronic interaction between the outer π charges of the triple bond and the 3pπ orbitals at the silicon atom. It is also likely that the conjugation of the geminal ethynyl groups leads to an enhancement of this bond contraction. Electrostatic interactions and the subsequent reduction of the covalent radius of the silicon atom may also contribute to this bond shortening. It has been found that the endocyclic Si–C bond length fits nicely within a scheme describing a monotonous decrease of the Si–C bond length with the increase of the electronegativity of the substituent in various geminally substituted silacyclobutanes.A series of related silacyclobutanes and acyclic diethynylsilanes have been studied by applying various ab initio methods and their optimized structures were compared to the structure of DESCB. Among these compounds are 1,1-dicyanosilacyclobutane (DCYSCB), which is isoelectronic to DESCB, 1,1-diethynylcyclobutane (DECB) which is isovalent to DESCB, monoethynylsilacyclobutane (MESCB) and monocyanosilacyclobutane (MCYSCB). Searching for reasonable support for the explanation of the structural results of DESCB we performed detailed natural population analysis as well as Mulliken population analysis (MPA) on DESCB and other related molecules. In contrast to the Mulliken charges, the natural atomic charges provided helpful information concerning the bonding properties in DESCB and the corresponding compounds. By varying the size of some basis sets, we could demonstrate the validity of the repeatedly discussed dependency of the Mulliken MPA on the basis set.For the performance of the quantum mechanical calculations we employed the following methods and basis sets: HF/6-31G(d,p), DFT/B3PW91/6-31G(d), DFT/B3PW91/6-311++G(d,p), MP2/6-31G(d,p) and MP2/6-311++G(d,p).     

Structure of a C26H28 alkene formed via titanium-promoted reductive dimerization of 4-(spirocyclopropyl)-pentacyclo[5.4.0.02,6.03,10.05,9]undecan-8-one

A spiro-fused C₂₆H₂₈ cage≓ dimer hydrocarbon contains two three-membered rings, two four-membered rings, eight five-membered rings, and two six-membered rings. The molecule is calculated by molecular mechanics to have 902 kJ/mol of strain energy distributed primarily between angle strain (457 kJ/mol) and torsional interactions (368 kJ/mol). Molecular mechanics calculations and a geometry-optimized ab initio calculation are used to analyze the bond lengths and bond angles in the molecule. There is one major discrepancy between observed and calculated distances.     

Investigation of C–HX (X=N, O, S) intramolecular hydrogen bond in 1-vinyl-2-(2′-heteroaryl)pyrroles by ab initio calculations

The C–HX (X=N, O, S) intramolecular hydrogen bond between the α-hydrogen of the vinyl group and the corresponding heteroatom in the series of 1-vinyl-2-(2′-heteroaryl)pyrroles was examined by ab initio calculations at the B3LYP/6-311(d,p) level. It was shown that the C–HN hydrogen bond is stronger than the C–HO hydrogen bond and the latter is, in turn, stronger than the C–HS hydrogen bond. This conclusion is supported by calculations of ¹H NMR chemical shieldings.     

Molecular mechanics (MM3) parameterization of hydroxylamine and methyl derivatives

Based on results of electron diffraction, gas phase infrared spectroscopy (IR), and MP2/6-31 + G* ab initio calculations, a set of molecular mechanics (MM3) parameters was developed for molecules containing the N(sp³)—O(sp³) moiety. Using this set of parameters, MM3 is able to reproduce structures (bond lengths and bond angles) and vibrational spectra satisfactorily. © 1994 by John Wiley & Sons, Inc.     

A molecular mechanics treatment of the anomeric effect

The tendency of C? O bond lengths to change as a function of the torsional angles at an acetal carbon has been included in a new version of the molecular mechanics program MM 2(82), based on the observed behavior of molecules of this class as indicated by ab initio calculations and experimental structural data. The experimental geometries and energies are reasonably well reproduced.     

Electronic structure and conformational properties of the amide linkage Part 5. Internal rotation and inversion in 1-formylaziridine

Internal rotation and nitrogen inversion in 1-formylaziridine (1) have been investigated by quantum mechanical (ab initio and MNDO) calculations, especially with respect to the variation of the geometry of the aziridine ring. While conformational stability is mainly determined by the n(N)/π(CO) interaction, the bond lengths within the ring are affected by the amount of interaction between the π(CO) orbital and the Walsh orbital ω_A. To separate the two types of interaction, calculations were also performed on formylcyclopropane (9). The torsional potential of 1 has a minimum close to the perpendicular conformation 1b. The two bisected conformations, 1a and 1c, are transition states for internal rotation. For nitrogen inversion, a barrier of 1.44 kcal mol⁻¹ (ab initio) was calculated. Calculations on 1-cyanoaziridine (7) gave inversion barriers of 5.81 (ab initio) and 12.31 kcal mol⁻¹ (MNDO). Probably due to methodical reasons the ab initio values seem to be too low, as calculations with different basis sets for aziridine indicate.     

Interaction-induced dipoles of hydrogen molecules colliding with helium atoms: a new ab initio dipole surface for high-temperature applications

We report new ab initio results for the interaction-induced dipole moments Δμ of hydrogen molecules colliding with helium atoms. These results are needed in order to calculate collision-induced absorption spectra at high temperatures; applications include modeling the radiative profiles of very cool white dwarf stars, with temperatures from 3500 K to 9000 K. We have evaluated the dipoles based on finite-field calculations, with coupled cluster methods in MOLPRO 2006 and aug-cc-pV5Z (spdfg) basis sets for both the H and He centers. We have obtained values of Δμ for eight H(2) bond lengths ranging from 0.942 a.u. to 2.801 a.u., for 15 intermolecular separations R ranging from 2.0 a.u. to 10.0 a.u., and for 19 different relative orientations. In general, our values agree well with earlier ab initio results, for the geometrical configurations that are treated in common, but we have determined more points on the collision-induced dipole surface by an order of magnitude. These results make it possible to calculate transition probabilities for molecules in excited vibrational states, overtones, and rotational transitions with ΔJ > 4. We have cast our results in the symmetry-adapted form needed for absorption line shape calculations, by expressing Δμ as a series in the spherical harmonics of the orientation angles of the intermolecular vector and of a unit vector along the H(2) bond axis. The expansion coefficients depend on the H(2) bond length and the intermolecular distance R. For large separations R, we show that the ab initio values of the leading coefficients converge to the predictions from perturbation theory, including both classical multipole polarization and dispersion effects.     

Molecular mechanics studies of ketene derivatives and related structures

The MM2 and MM3 force fields have now been parameterized for ketene and its various derivatives. With the addition of the C_sp ? O bond stretching and C?C_sp?O bond bending parameters, calculations were performed on ketene and six substituted ketene compounds. The MM2 results are quite good with only minimal errors in the calculation of C? H bond lengths and H-C-H bond angles. Additionally, C? F bond parameters in MM2 have been re-adjusted to give better results in monofluorinated species, but, unfortunately, resulting in greater error in the polyfluorinated compounds. The results of geometry calculations by MM3 are similar to those obtained by MM2 with the exception of a significant improvement in the geometry of dimethylketene. The MM3 vibrational frequencies calculated in this study are also in good agreement with available experimental and ab initio results with the exception of a few low frequency in-and out-of-plane bending modes. © 1992 by John Wiley & Sons, Inc.     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

Unimolecular rearrangement of α-silylcarbenium ions. An ab initio study

The unimolecular rearrangements of hydrogen, methyl and phenyl groups at the Si atom in α-silylcarbenium ions have been investigated using an ab initio molecular orbital method. MP2/6–31 + G^*//HF/6–31G^* calculations predict that all three groups migrate from the Si to an adjacent C_α with no energy barrier. Thus, the silicenium ion is the only stable species in each potential energy surface. The conformation of the benzylsilicenium ion, (C₆H₅)CH₂−SiH₂⁺, indicates that the phenyl ring is significantly bent toward the silyl cationic center in order to interact with the vacant 3p(Si⁺) orbital. In contrast to MP2 results, Hartree-Fuck calculations (both HF/3–21G^* and HF/6–31G^* levels) predict small energy barriers for 1,2-migrations of H and Me (1.4 kcal mol⁻¹ for H migration, and 1.5 kcal mol⁻¹ for Me migration, respectively, at the HF/6–31G^* level). This difference provides convincing evidence that the incorporation of electron correlation is of particular importance in describing the potential energy surface for the rearrangement of α-silylcarbenium ions to silicenium ions. The results of the calculations have also been applied to the possible rearrangement mechanism of α-chlorosilanes to chlorosilanes, assuming that the experimental conditions are favorable toward the generation of ionic species. Various factors which may govern the migratory aptitudes of various R groups, i.e. (1) activation energies, (2) overall reaction energies and (3) the conformational preference of reactants have been investigated. The calculated activation energy obtained, namely the energy for the generation of the silicenium ion and the C⁻¹ ion from an α-chlorosilane, is consistent with the experimental migratory aptitude in the gas phase observed in mass spectrometers.     

Experimental and Ab Initio Equilibrium Structures of cis-Thionylimide,HNSO: Estimation of the Laurie Correction

The equilibrium structure of cis-thionylimide, HNSO has been determined using high-level ab initio calculations and various experimental procedures. The Laurie correction to the N-H bond length is discussed and it is shown that it can be estimated by ab initio methods. The results are found to be in good agreement and the best equilibrium structure is (in Å for the bond lengths): r(N–H) = 1.020(1), r(N–S) = 1.510(2), r(S–O) = 1.448(1), (HNS) = 115.94(39)°, and (NSO) = 120.44(10)°.     

DTMM and COSMIC molecular mechanics parameters for alkylsilanes

Harmonic DTMM and COSMIC molecular mechanical force fields were extended to cover aliphatic silanes. Structures of the thirty-five molecules were examined and they generally fit approximately within experimental error. For these, root mean square deviations between the calculated and the experimentally determined geometrical data were found to be less than 0.05 Å and 1.7 deg. for bond lengths and bond angles, respectively. The accuracy of these force fields was also verified in comparison with the other force fields (MM2, MM3), semi-empirical (PM3), ab initio (HF,MP2) and density functional (B3LYP, BP86) quantum chemical methods. The relative stability of the seventeen conformations of seven simple molecules was investigated and the limitations of the DTMM software is discussed.     

Conformational study of cis-1,4-di-tert-butylcyclohexane by dynamic NMR spectroscopy and computational methods. Observation of chair and twist-boat conformations

[reaction: see text] Low-temperature 13C NMR spectra of cis-1,4-di-tert-butylcyclohexane (1) showed signals for the twist-boat (1a) and chair (1b) conformations. 13C NMR signals were assigned to specific carbons based on the different populations, different symmetries (time-averaged C(2v) for 1a and time-averaged C(s) for 1b), and calculated chemical shifts (GIAO, HF/6-311+G*). In addition to slow ring inversion and interconversion of the chair and twist-boat conformations, slow rotation of the tert-butyl groups was found. Most of the expected 13C peaks were observed. Free-energy barriers of 6.83 and 6.35 kcal/mol were found for interconversion of 1a (major) and 1b (minor) at -148.1 degrees C. Conformational space was searched with Allinger's MM3 and MM4 programs, and free energies were obtained for several low-energy conformations 1a-c. Calculations were repeated with ab initio methods up to the HF/6-311+G* level. Molecular symmetries, relative free energies, relative enthalpies and entropies, frequencies, and NMR chemical shifts were obtained. A boat conformation (1d; C(2v) symmetry) was generated and optimized as a transition state by ab initio, MM3, and MM4 calculations.     

Molecular mechanics (MM4) study of fluorinated hydrocarbons

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.     

Ab initio and molecular mechanics (MM3) calculations of protonated-neutral diamine hydrogen bonds

Ab initio calculations of cation-neutral diamine complexes have been carried out at the MP2/6-311+G** level. The geometry and energetics of the charge-reinforced hydrogen bond are analyzed with respect to the alkyl substitution of both the protonated and neutral nitrogen atoms, and these results have been used to improve the quality of the MM3(2000) force field. In addition, specialized hydrogen bond parameters optimized for MM3(2000) are presented. These parameters allow very accurate gas-phase modeling of the charge-neutral diamine environment. Molecular mechanics calculations can model effectively protonated amine-neutral amine hydrogen bonds in the gas phase and solution (continuum dielectric) through a combination of charge-dipole interactions and explicit hydrogen-bonding terms.