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.     

Alcohols,ethers, carbohydrates,and related compounds. I. The MM4 force field for simple compounds

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).     

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     

An improved empirical force field for saturated hydrocarbons

A new molecular mechanics force field for alkanes is presented. The force field aims to eliminate some identified failures of the well-known MM2 force field. The new energy function gives an improved prediction of the rotational barriers of highly congested molecules, a better calculation of short nonbonded contacts, and the correct reproduction of bond elongation in small torsion angles. The calculation of sublimation enthalpies is also improved. The standard deviation of the formation enthalpies for a set of 54 compounds is 0.63 kcal/mol; this compares with the reported value of 0.42 calculated with MM2 and MM3 for different sets. The force field parameters were obtained using a least squares method.     

Chlorosilanes: Development and application of MM2 force field parameters

The geometries, relative conformational energies, and dipole moments of mono and polychlorosilanes have been calculated using ab initio molecular orbital (MO) theory. Calculations at the HF/3–21G(*) level, with the exception of dipole moments, give reasonable agreement with experimental data. A new MM2 force field for chlorosilanes, which includes terms for bond length shortening and bond angle compression due to the attachment of electronegative Cl atoms, has been developed on the basis of experimental and ab initio results. The new force field is generally successful in predicting structural parameters, but is unable to reproduce the dipole moments of several model systems. While dipole moment predictions are not the authors' main interest, this failure defines a shortcoming in the MM2 method. The new parameters have been applied to problems in the prediction of stereochemistries of cyclic systems, and compared with experimental results where data are available.     

Molecular mechanics (MM4) calculations on carbonyl compounds part I: aldehydes

Aliphatic aldehydes have been studied with the aid of the MM4 force field. The structures, moments of inertia, vibrational spectra, conformational energies, barriers to internal rotation, and dipole moments have been examined for six compounds (nine conformations). MM4 parameters have been developed to fit the indicated quantities to the wide variety of experimental data. Ab initio (MP2) and density functional theory (B3LYP) calculations have been used to augment and/or replace experimental data, as appropriate. Because more, and to some extent, better, data have become available since MM3 was developed, it was anticipated that the overall accuracy of the information calculated with MM4 would be better than with MM3. The best single measure of the overall accuracy of a force field is the accuracy to which the moments of inertia of a set of compounds (from microwave spectroscopy) can be reproduced. For all of the 20 moments (seven conformations) experimentally known for the aldehyde compounds, the MM4 rms error is 0.30%, while with MM3, the most accurate force field presently available, the rms error over the same set is 1.01%. The calculation of the vibrational spectra was also improved overall. For the four aldehydes that were fully analyzed (over a total of 78 frequencies), the rms errors with MM4 and MM3 are 18 and 38 cm^?1, respectively. These improvements came from several sources, but the major ones were separate parameters involving the carbonyl carbon for formaldehyde, the alkyl aldehydes and the ketones, and new crossterms featured in the MM4 force field that are not present in the MM3 version. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1396–1425, 2001     

Development of a molecular mechanics (MM2) force field for α-chlorosilanes

Molecular mechanics (MM2) parameters for silanes which have a Si-C-Cl fragment have been developed based on available experimental data and ab initio molecular orbital (MO) calculations. Molecular properties, mainly rotational barriers and geometries, of α-chlorosilanes have been studied using our new MM2 parameter set. Changes in the Si-C bond lengths and several bond angles of α-chlorosilanes due to the additional attachment of polar atom(s) have been investigated utilizing ab initio calculations. An electronegativity correction to both bond lengths and angles helps MM2 to reproduce results from ab initio calculations. The new force field has been applied to the conformational analysis of l-(chloromethyl)-1,2-dimethylsilacyclopentane, a model used in our studies of rearrangements of α-halosilanes.     

Force‐field parameters of the Ψ and Φ around glycosidic bonds to oxygen and sulfur atoms

The Ψ and Φ torsion angles around glycosidic bonds in a glycoside chain are the most important determinants of the conformation of a glycoside chain. We determined force‐field parameters for Ψ and Φ torsion angles around a glycosidic bond bridged by a sulfur atom, as well as a bond bridged by an oxygen atom as a preparation for the next study, i.e., molecular dynamics free energy calculations for protein‐sugar and protein‐inhibitor complexes. First, we extracted the Ψ or Φ torsion energy component from a quantum mechanics (QM) total energy by subtracting all the molecular mechanics (MM) force‐field components except for the Ψ or Φ torsion angle. The Ψ and Φ energy components extracted (hereafter called “the remaining energy components”) were calculated for simple sugar models and plotted as functions of the Ψ and Φ angles. The remaining energy component curves of Ψ and Φ were well represented by the torsion force‐field functions consisting of four and three cosine functions, respectively. To confirm the reliability of the force‐field parameters and to confirm its compatibility with other force‐fields, we calculated adiabatic potential curves as functions of Ψ and Φ for the model glycosides by adopting the Ψ and Φ force‐field parameters obtained and by energetically optimizing other degrees of freedom. The MM potential energy curves obtained for Ψ and Φ well represented the QM adiabatic curves and also these curves' differences with regard to the glycosidic oxygen and sulfur atoms. Our Ψ and Φ force‐fields of glycosidic oxygen gave MM potential energy curves that more closely represented the respective QM curves than did those of the recently developed GLYCAM force‐field. © 2009 Wiley Periodicals, Inc., J Comput Chem, 2009     

The MM3 force field for amides,polypeptides and proteins

The potential functions for simple amides, several peptides and a small protein have been worked out for the MM3 force field. Structures and energies were fit as previously with MM2, but additionally, we fit the vibrational spectra of the simple amides (average rms error over four compounds, 34 cm^?1), and examined more carefully electrostatic interactions, including charge-charge and charge-dipole interactions. The parameters were obtained and tested by examining four simple amides, five electrostatic model complexes, two dipeptides, six crystalline cyclic peptides, and the protein Crambin. The average root-mean-square deviation from the X-ray structures for the six cyclic peptide crystals was only 0.10 Å for the nonhydrogen atomic positions, and 0.011 Å, 1.0°, and 4.9° for bond lengths, bond angles, and torsional angles, respectively. The parameter set was then further tested by minimizing the high resolution crystal structure of the hydrophobic protein Crambin. The resultant root-mean-square deviations for the non-hydrogen atomic data, in the presence of the crystal lattice, are 0.22 Å, 0.023 Å, 2.0°, and 6.4° for coordinates, bond lengths, bond angles, and torsional angles, respectively.     

Molecular mechanics (MM4) study of amines

The MM4 force field has been extended to include aliphatic amines. About 20 amines have been examined to obtain a set of useful molecular mechanics parameters for this class. The vibrational spectra of seven amines (172 frequencies) calculated by MM4 have an overall rms error of 27 cm(-1), compared with corresponding MM4 value of 24 cm(-1) for alkanes. The rms and signed average errors of the moments of inertia of nine simple amines compared with the experimental data were 0.18% and -0.004%, respectively. The heats of formation of 30 amines were also studied. The MM4 weighted standard deviation is 0.41 kcal/mol, compared with experiment. Electronegativity effects occur in the hydrocarbon portion of an amine from the nitrogen, and are accounted for by including electronegativity induced changes in bond lengths and angles, and induced dipole-dipole interactions in the molecule. Negative hyperconjugation results from the presence of the lone pair of electrons on nitrogen, and leads to the Bohlmann bands in the infrared, and also to strong and unusual geometric changes in the molecules (Bohlmann effect), all of which are fairly well accounted for. The conformational energies in amines appear to be less straightforward than those for most other classes of molecules, apparently because of the Bohlmann effect, and these are probably not yet completely understood. In general, the agreement between the MM4 calculated results and the available data is reasonably good.     

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.     

Alcohols,ethers, carbohydrates,and related compounds. II. The anomeric effect

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.     

利用浮动电荷力场对N-甲基乙酰胺水溶液的分子动力学模拟

利用原子键电负性均衡结合分子力场方法(ABEEM/MM)对N-甲基乙酰胺(NMA)分子的水溶液体系进行了分子动力学模拟. 与经典的力场模型相比, 该方法中的静电势包含了分子内和分子间的静电极化作用, 以及分子内电荷转移影响, 同时加入了化学键等非原子中心电荷位点, 合理体现了分子中的电荷分布. 相对其它极化力场模型, 该模型具有计算量较小的特点. 在该模型下对NMA纯溶液和其水溶液体系进行了分子动力学模拟, 得到的径向分布函数、汽化热和偶极矩等物理量与实验值和其它极化力场方法符合很好, 合理描述了溶质与溶剂之间的静电极化和分子内的电荷转移.     

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.     

Molecular mechanics calculations on carbonyl compounds. II. Open‐chain ketones

Open‐chain aliphatic ketones were studied with the molecular mechanics (MM4) force field. A total of seven compounds were examined. Structures were well fit, including moments of inertia. Rotational barriers, vibrational spectra, and dipole moments were also well fit. The overall root mean square errors for MM3 and MM4 were 0.27 and 0.18%, respectively, for the six moments of inertia (known experimentally for two compounds) and 31 and 20 cm^?1, respectively, for the vibrational frequencies (over 99 weighted modes). © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1426–1450, 2001     

从头算和ABEEMσπ/MM对(CH_3OH)_n(n=3～12)和[Na(CH_3OH)_n]~+(n=3～6)体系的研究

采用从头算方法(ab initio)和原子-键电负性均衡浮动电荷分子力场方法(ABEEMσπ/MM),对甲醇团簇(CH_3OH)n(n=3~12)和[Na(CH_3OH)_n]~+(n=3~6)体系的结构、电荷分布和结合能进行研究.依据从头算结果构建上述体系的ABEEMσπ/MM浮动电荷势能函数,并确定相关参数.结果表明,ABEEMσπ/MM所获得的结构和结合能等均优于OPLS/AA力场,并与从头算结果相符,其中键长的平均绝对偏差(AAD)小于0.004 nm,键长、键角和结合能的相对均方根偏差(RRMSD)分别小于3.8%,1.7%和6.8%;电荷分布与从头算结果的线性相关系数均大于0.99.     

DELPHI有机共轭体系MM力场的优化和应用

运用我们建议的SC-HMO方法优化了有机共轭烯烃的DELPHIMM力场。将优化所得力场和参数用于计算各类型101个共轭化合物的分子几何构型、生成热、偶极矩和IR频率等, 结果可与MM3法相比。     

Comparative performance of MM3(92) and two TINKER MM3 versions for the modeling of carbohydrates

The 1992 version of MM3 was largely used for modeling mono-, di-, and trisaccharides. In later versions of MM3 improvements were made in some parameters that may be important for carbohydrates. This corrected MM3 force field is part of the Tinker package, freely available (as its 4.1 version), and included in the Chem 3D Ultra 8.0 package (as the 3.7 version). The latter version lacks the corrections to the standard bond lengths produced by electronegativity and anomeric effects, whereas the Tinker 4.1 version only lacks the latter correction. The present work compares the performance of the three MM3 versions (and in some cases, DFT and/or HF/ab initio procedures) on several carbohydrate model problems as the chair and rotamer equilibria in 2-hydroxy- and 2-methoxytetrahydropyran, hydrogen bonding in cis-2,3-dihydroxytetrahydropyran, and the potential energy surfaces around the glycosidic bonds of two sulfated disaccharides and two trisaccharides. Tinker MM3 can be used accurately to estimate carbohydrate energies and geometries, and-with the help of some programming-to pursue studies on the potential energy surfaces of di- and trisaccharides. In most cases results obtained using the three MM3 versions are similar, although large energy differences are obtained when comparing a rotameric distribution around a O-C-O-H dihedral, which is almost forced to the exo-anomeric position by the Tinker versions. In other systems smaller energy differences are found, but they can nevertheless lead to a different global minimum when comparing conformers of similar energy. MM3(92) establishes better the differences between the bond lengths in both anomers, as an expected expression of the anomeric correction.     

A New Method for Molecular Mechanics

A promising new method for optimizing molecular structures is described. In place of the terms involving bond angles and torsion angles, used in the force fields of conventional molecular mechanics, two-body central forces between atoms are used exclusively, resulting in a considerable computational advantage. The program STRFIT, using this method has been tested by comparing geometries obtained with those found using the popular molecular mechanics program MM2 (Allinger) for a variety of cyclic and acyclic molecules. For unstrained molecules, the difference in steric energy between geometries refined by STRFIT and MM2 is only a few tenths of a kilocalorie and up to about a kilocalorie for strained molecules. Geometry optimization with STRFIT, to a structure that is slightly higher in energy than the structure arrived at by MM2 starting from the same initial starting geometry, is three to eight times faster. A complete new package of programs for conveniently and rapidly doing molecular mechanics calculations is described. It includes a convenient algorithm for the input of approximate molecular structures, a rapid structure-optimizing module using a pure Central force-field approach, and a drawing program designed for use with a dot-matrix printer or a laser printer.     

Molecular mechanics conformational analysis of tylosin

The conformations of the 16-membered macrolide antibiotic tylosin were studied with molecular mechanics (AMBER* force field) including modelling of the effect of the solvent on the conformational preferences (GB/SA). A Monte Carlo conformational search procedure was used for finding the most probable low-energy conformations. The present study provides complementary data to recently reported analysis of the conformations of tylosin based on NMR techniques. A search for the low-energy conformations of protynolide, a 16-membered lactone containing the same aglycone as tylosin, was also carried out, and the results were compared with the observed conformation in the crystal as well as with the most probable conformations of the macrocyclic ring of tylosin. The dependence of the results on force field was also studied by utilizing the MM3 force field. Some particular conformations were computed with the semiempirical molecular orbital methods AM1 and PM3.