A molecular mechanics (MM3) study of fluorinated hydrocarbons

A molecular mechanics study of small saturated hydrocarbons (up to C-6), substituted by up to six fluorines was carried out with the MM3 force field. Perfluorobutane and Teflon were also studied. A parameter set was 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 in good agreement with experiment when only one or two fluorines are present, but some rather large discrepancies were noted when the F/H ratio becomes high. These can be taken into account only by using a force field more complicated than MM3. Some of the requirements of such a force field are delineated. Some pertinent ab initio results are also reported in this article.     

Development of an internal searching algorithm for parameterization of the MM2/MM3 force fields

Application of Allinger's MM2/MM3 force fields to molecules of real interest is frequently hindered by the lack of parameters for various heterocyclic systems and for poly-functionalized molecules. A common approach to this problem is to manually choose missing parameters “by analogy” with those that are part of the force field's internal parameter set. Naturally, this is generally attempted only by those possessing extensive experience with force fields. In order to use the MM2/MM3 force fields to study herbicides, an algorithm has been developed to automate this process for the non MM2 specialist. Using a set of “relative cost” criteria for atom type replacement, the algorithm searches the force field parameter set and selects the most appropriate parameters for a given molecule whose MM2 output file contains “missing parameter” errors. The program selects parameter error messages from a standard MM2 output file, finds analogous parameters, asks the user to verify their appropriateness and creates a standard MM2 parameter deck for the molecule of interest.     

Parameter analysis and refinement toolkit system and its application in MM3 parameterization for phosphine and its derivatives

The multiparameter multistep relaxation (MPMSR) method, a routine within a new suite of parameterization programs entitled parameter analysis and refinement toolkit system (PARTS), was developed to assist in the development of molecular mechanics (MM3 and MM2) force field parameters and represents an ongoing effort in our laboratories to generate more accurate force fields in shorter times. In contrast to other computerized parameterization approaches, this method simulates intuition guided trial-and-error and has been used successfully within our laboratories to develop MM2 and MM3 force fields. The primary aim of this approach is to minimize human inspection time and effort, with simultaneous improvement in the efficiency and accuracy of the parameterization process. In an effort to validate the generality of the MPMSR method, a well parameterized data set of phosphine derivatives was reexamined. With the identical set of training molecules used in the original MM3 phosphine parameterization and with minimal human intervention, MPMSR shortened the process from several months to approximately five days. Although the previous phosphine force field is well parameterized, the newly generated MPMSR set of parameters has achieved an overall better fit to the experimentally observed data and ab initio calculations. © 1996 by John Wiley & Sons, Inc.     

从头算和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.     

Ab initio and molecular mechanics (MM2 and MM3) calculations of positively charged conjugated nitrogen-containing compounds

Ab initio calculations have been carried out on s-trans-N-vinylmethyleneammonium, pyridinium, and related compounds to obtain rotational barriers, structures, and vibrational frequencies. The restricted Hartree-Fock (RHF) level of theory with 6-31G** basis set was used for these calculations. In addition, the MM2(91) and MM3(94) force fields have been parameterized to calculate these positively charged nitrogen-containing compounds. A bond order term was incorporated in the force field to reproduce accurately the rotational barriers of s-trans-N-vinylmethyleneammonium and related compounds. Molecular mechanics geometries and vibrational frequencies compare well with those calculated by ab initio methods. © 1996 by John Wiley & Sons, Inc.     

Partial hessian fitting for determining force constant parameters in molecular mechanics

We present a new protocol for deriving force constant parameters that are used in molecular mechanics (MM) force fields to describe the bond‐stretching, angle‐bending, and dihedral terms. A 3 × 3 partial matrix is chosen from the MM Hessian matrix in Cartesian coordinates according to a simple rule and made as close as possible to the corresponding partial Hessian matrix computed using quantum mechanics (QM). This partial Hessian fitting (PHF) is done analytically and thus rapidly in a least‐squares sense, yielding force constant parameters as the output. We herein apply this approach to derive force constant parameters for the AMBER‐type energy expression. Test calculations on several different molecules show good performance of the PHF parameter sets in terms of how well they can reproduce QM‐calculated frequencies. When soft bonds are involved in the target molecule as in the case of secondary building units of metal‐organic frameworks, the MM‐optimized geometry sometimes deviates significantly from the QM‐optimized one. We show that this problem is rectified effectively by use of a simple procedure called Katachi that modifies the equilibrium bond distances and angles in bond‐stretching and angle‐bending terms. © 2016 Wiley Periodicals, Inc.     

Dimerization of thymol blue in solution: Theoretical evidence

The possibility of dimerization of thymol blue was addressed by ab initio and force field calculations. In agreement with experimental information, a dimer forming symmetrical chemical environments for hydrogen bond formation was determined. This dimer is stable in vacuum and aqueous media and corresponds to the same protonated state proposed by the experiment. A comparison of the CVFF and MM3 force fields and ab initio results shows the suitability of CVFF to qualitatively describe this system.     

Ab initio study of the molecular structures, force fields, and vibrational spectra of alkaline metal fluoride dimers MM′F2 (M, M′=Li, Na, K)

Ab initio calculations of the equilibrium geometrical parameters, force constants, and IR vibration frequencies and intensities of Li₂F₂, Na₂F₂, K₂F₂, LiNaF₂, LiKF₂, and NaKF₂ are reported. The calculations use the Hartree-Fock-Roothaan method and second-order Möller-Plesset perturbation theory along with configuration interaction theory including singly and doubly excited configurations and corrections for quartic excitations with basis sets of grouped Gaussian functions: Li — (9s3p1d/4s3p1d), Na — (12s8p1d/6s4p1d), K — (14s11p3d/9s8p3d), F — (9s5p1d/4s2p1d). According to the results of calculations, the equilibrium structures of the molecules are planar cyclic structures of D_2h (for M₂F₂) and C_2v (for MM′F₂) symmetries. The linear configurations M-F-M′-F (of C_∞v symmetry) are 70–190 kJ/mole less stable than the cyclic ones; for all molecules except M-F-K-F, these configurations correspond to local minima on the potential energy surface. The role of correlation effects in ab initio calculations of the geometry, force fields, and IR characteristics of molecules with highly polar chemical bonds is discussed. The theoretical force fields of the molecules are represented in canonical form in a system of redundant natural vibrational coordinates. The force fields of MM′F₂ molecules are studied. The results of the ab initio calculations are compared with the experimental structural data and vibrational spectra available in the literature.     

Molecular mechanics (MM3) calculations on lithium amide compounds

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.     

Shapes of sulfur, oxygen, and nitrogen mustards

Thorough conformational analyses have been performed on representative sulfur, oxygen, and nitrogen mustards. A total of 23, 18, and 38 unique conformers have been located for SM, OM, and NM, respectively, at the MP2/aug-cc-pVDZ level of theory. Despite the fact that these molecules differ only in the identity of the central heteroatom, comparison of their low energy conformations reveals that the shapes they adopt are distinctive to each molecule. Potential energy surfaces for CH(2)-X (X = S, O, and N-CH(3)) and CH(2)-CH(2) bond rotations are presented and, where possible, compared with dihedral angle distributions observed in crystal structure data. These results were used to benchmark and improve the performance of the MM3 and MMFF94 force fields.     

A comparison of conformational energies calculated by molecular mechanics (MM2(85), Sybyl 5.1, Sybyl 5.21, and ChemX) and semiempirical (AM1 and PM3) methods

Conformational energies of different conformers have been calculated for a series of molecules using various molecular mechanics and semiempirical methods. The quality of the force fields has also been tested by calculating barriers to rotation about carbon-carbon bonds. The molecular mechanics force fields used are MM2(85), Sybyl 5.1, Sybyl 5.21, and ChemX, ver. Jan 89. The semiempirical methods used are AM1 and PM3. Molecules with different functional groups, for which good experimental data exist, have been selected. The semiempirical methods generally calculate barriers to rotation which are lower than the experimentally determined. The conformational energies for hydrocarbons are reasonably well reproduced by all tested methods although MM2(85) gives the quantitatively best agreement with experiments. For compounds containing oxygen, nitrogen and halogens MM2(85) gives results which are in best agreement with the experimentally determined values.     

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.     

Mechanism of proteolysis in matrix metalloproteinase‐2 revealed by QM/MM modeling

The mechanism of enzymatic peptide hydrolysis in matrix metalloproteinase‐2 (MMP‐2) was studied at atomic resolution through quantum mechanics/molecular mechanics (QM/MM) simulations. An all‐atom three‐dimensional molecular model was constructed on the basis of a crystal structure from the Protein Data Bank (ID: 1QIB), and the oligopeptide Ace‐Gln‐Gly～Ile‐Ala‐Gly‐Nme was considered as the substrate. Two QM/MM software packages and several computational protocols were employed to calculate QM/MM energy profiles for a four‐step mechanism involving an initial nucleophilic attack followed by hydrogen bond rearrangement, proton transfer, and C? N bond cleavage. These QM/MM calculations consistently yield rather low overall barriers for the chemical steps, in the range of 5–10 kcal/mol, for diverse QM treatments (PBE0, B3LYP, and BB1K density functionals as well as local coupled cluster treatments) and two MM force fields (CHARMM and AMBER). It, thus, seems likely that product release is the rate‐limiting step in MMP‐2 catalysis. This is supported by an exploration of various release channels through QM/MM reaction path calculations and steered molecular dynamics simulations. © 2015 Wiley Periodicals, Inc.     

Efficient treatment of out-of-plane bend and improper torsion interactions in MM2, MM3, and MM4 molecular mechanics calculations

Simple and very efficient formulas are presented for four-body out-of-plane bend (used in MM2 and MM3 force fields) and improper torsion (used in the MM4 force field) internal coordinates and their first and second derivatives. The use of a small set of bend and stretch intermediates allows for order of magnitude decreases in calculation time for potential energies and their first and second derivatives, which are required in molecular mechanics calculations. The formulas are eminently suitable for use in molecular simulations of systems with complicated bond networks. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1804–1811, 1997     

Solvation free energies calculated using the GB/SA model: Sensitivity of results on charge sets,protocols, and force fields

The sensitivity of aqueous solvation free energies (SFEs), estimated using the GB/SA continuum solvent model, on charge sets, protocols, and force fields, was studied. Simple energy calculations using the GB/SA solvent model were performed on 11 monofunctional organic compounds. Results indicate that calculated SFEs are strongly dependent on the charge sets. Charges derived from electrostatic potential fitting to high level ab initio wave functions using the CHELPG procedure and “class IV” charges from AM1/CM1a or PM3/CM1p calculations yielded better results than the corresponding Mulliken charges. Calculated SFEs were similar to MC/FEP energies obtained in the presence of explicit TIP4P water. Further improvements were obtained by using GVB/6-31G** and MP2/6-31+G** (CHELPG) charge sets that included correlation effects. SFEs calculated using charge sets assigned by the OPLSA* force field gave the best results of all standard force fields (MM2*, MM3*, MMFF, AMBER*, and OPLSA*) implemented in MacroModel. Comparison of relative and absolute SFEs computed using either the GB/SA continuum model or MC/FEP calculations in the presence of explicit TIP4P water showed that, in general, relative SFEs can be estimated with greater accuracy. A second set of 20 mono- and difunctional molecules was also studied and relative SFEs estimated using energy minimization and thermodynamic cycle perturbation (TCP) protocols. SFEs calculated from TCP calculations using the GB/SA model were sensitive to bond lengths of dummy bonds (i.e., bonds involving dummy atoms). In such cases, keeping the bond lengths of dummy bonds close to the corresponding bond lengths of the starting structures improved the agreement of TCP-calculated SFEs with energy minimization results. Overall, these results indicate that GB/SA solvation free energy estimates from simple energy minimization calculations are of similar accuracy and value to those obtained using more elaborate TCP protocols. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 769–780, 1998     

环多肽晶体的浮动电荷极化力场模拟

利用原子键电负性均衡结合分子力场方法(ABEEM/MM)对五种环多肽晶体进行了研究. 与传统力场相比, 该方法中的静电势包含了分子内和分子间的静电极化作用, 以及分子内电荷转移影响, 同时加入了化学键等非原子中心电荷位点, 合理地体现了分子中的电荷分布. 相对其他极化力场模型, 具有计算量较小的特点. 该模型下计算得到的环多肽分子单元相对实验测得的结构的原子位置、氢键长度和二面角的均方根偏差分别为0.009 nm、0.013 nm和5.16°, 能够很好地重复实验结果. 总体上, 其结果优于或相当于其他力场模型, 适用于对实际蛋白质体系的模拟和研究.     

Structural Dynamics of Some Large-Ring Cyclodextrins. A Molecular Dynamics Study: An Analysis of Force Field Performance

Molecular dynamics simulations were performed on four large ring cyclodextrins (LR-CDs), CD14, CD21, CD26 and CD28, both in gas phase and in water solution. Four different force fields (parm94, parm99, glycam2000a and MM3*) were used to examine their dependence with the results. The differences were not significant when parm99 and glycam2000a were used. Parm94 and MM3* results differ considerably due to inadequate bending parameters (parm94) or due to a strong stabilization by the intramolecular hydrogen bonds between hydroxyl groups. Simulations in the gas phase show that the CDs are elongated and twisted, and that cavities typical for the native CDs (CDn, n = 6, 7, 8) are not present. Simulations in water solution produced average structures that do not correspond to the conformations in the crystalline state. The LR-CDs are highly flexible and this could favor the formation of inclusion complexes through the intermediation of some other molecules that could fit to the smaller and more specific cavities they have.On sabbatical leave from the Institute of Organic Chemistry with Center of Phytochemistry, Bulgarian Academy of Science, ul. Acad. G. Bonchev, bloc 9, 1113 Sofia, Bulgaria     

Computational study on the ring distortions of 1,4-dicyanobenzene in the gas phase and in the crystal

A molecular mechanics software enhanced to perform empirical energy calculations on crystals (KESSHOU) was further developed to handle intermolecular electrostatic interactions as well. The packing of the molecules of 1,4-dicyanobenzene and 1,4-diisocyanobenzene in the crystal was studied. The role of the van der Waals and the electrostatic interactions in the balance of nonbonded atom-atom interactions is analyzed. The packing forces are dominated by van der Waals forces. The electrostatic interactions have higher stabilizing contribution for the dicyano isomer than for the diisocyano form. The dependence of the results on the size of the crystal, the molecular mechanics force field (MM2 vs MM3), and the dielectric constant are also assessed. Ab initio MP2/6–311G** geometries of the isolated molecules are in accordance with the observed benzene ring distortions determined by electron diffraction.     

QuanPol: A full spectrum and seamless QM/MM program

The quantum chemistry polarizable force field program (QuanPol) is implemented to perform combined quantum mechanical and molecular mechanical (QM/MM) calculations with induced dipole polarizable force fields and induced surface charge continuum solvation models. The QM methods include Hartree–Fock method, density functional theory method (DFT), generalized valence bond theory method, multiconfiguration self‐consistent field method, Møller–Plesset perturbation theory method, and time‐dependent DFT method. The induced dipoles of the MM atoms and the induced surface charges of the continuum solvation model are self‐consistently and variationally determined together with the QM wavefunction. The MM force field methods can be user specified, or a standard force field such as MMFF94, Chemistry at Harvard Molecular Mechanics (CHARMM), Assisted Model Building with Energy Refinement (AMBER), and Optimized Potentials for Liquid Simulations‐All Atom (OPLS‐AA). Analytic gradients for all of these methods are implemented so geometry optimization and molecular dynamics (MD) simulation can be performed. MD free energy perturbation and umbrella sampling methods are also implemented. © 2013 Wiley Periodicals, Inc.     

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