CPMD/GULP QM/MM interface for modeling periodic solids: Implementation and its application in the study of Y‐zeolite supported Rhn clusters

We report here the development of hybrid quantum mechanics/molecular mechanics (QM/MM) interface between the plane‐wave density functional theory based CPMD code and the empirical force‐field based GULP code for modeling periodic solids and surfaces. The hybrid QM/MM interface is based on the electrostatic coupling between QM and MM regions. The interface is designed for carrying out full relaxation of all the QM and MM atoms during geometry optimizations and molecular dynamics simulations, including the boundary atoms. Both Born–Oppenheimer and Car–Parrinello molecular dynamics schemes are enabled for the QM part during the QM/MM calculations. This interface has the advantage of parallelization of both the programs such that the QM and MM force evaluations can be carried out in parallel to model large systems. The interface program is first validated for total energy conservation and parallel scaling performance is benchmarked. Oxygen vacancy in α‐cristobalite is then studied in detail and the results are compared with a fully QM calculation and experimental data. Subsequently, we use our implementation to investigate the structure of rhodium cluster (Rh_n; n = 2 to 6) formed from Rh(C₂H₄)₂ complex adsorbed within a cavity of Y‐zeolite in a reducible atmosphere of H₂ gas. © 2016 Wiley Periodicals, Inc.     

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.     

Mixed ab initio QM/MM modeling using frozen orbitals and tests with alanine dipeptide and tetrapeptide

Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree–Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by “breaking” bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. Results from tests of the method to find relative conformational energies and geometries of alanine dipeptides and alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1468–1494, 1999     

NMR和分子力学法研究溶液中有机分子结构

在分子力学能量优化项中加入由ＮＭＲ实验得到的几何参数二面角（相应于~３Ｊ_（ＨＨ）耦合常数）及质子间距所构造的约束函数，使分子力学计算得到的能量极小点的分子构象，其耦合常数（~３Ｊ_（ＨＨ））及质子间距同ＮＭＲ实验结果相符合。文中以（±）－３－（４＇－甲苯基）－１－氮杂双环［２．２．２］－３－辛醇为例，采用自编的ＮＭＲ参数约束的分子力学计算程序ＭＭ２ＮＪ对其构象进行了计算。所得结果与Ｘ射线衍射法比较，有较好的一致性     

NMR spin-spin coupling constants in polymethine dyes as polarity indicators

Herein, we explore the use of spin-spin coupling constants (SSCCs) in merocyanine (MCYNE) dyes as indicators of polarity. For this purpose, we use Car-Parrinello hybrid quantum mechanics/molecular mechanics (QM/MM) to determine the structures of MCYNE in solvents of different polarity, followed by computations of the SSCCs by using QM/MM linear-response theory. The molecular geometry of MCYNE switches between neutral, cyanine-like, and zwitterionic depending on the polarity of the solvent. This structural variation is clearly reflected in the proton SSCCs in the polymethine backbone, which are highly sensitive to the dielectric nature of the environment; this mechanism can be used as a "polarity indicator" for different microenvironments. This result is highlighted by computing the SSCCs of the MCYNE probe in the cavity of the beta-lactoglobulin protein. The computed SSCCs clearly indicate a non-polar hydrophobic dielectric nature of this cavity.     

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.     

NMR Spin–Spin Coupling Constants in Polymethine Dyes as Polarity Indicators

Herein, we explore the use of spin–spin coupling constants (SSCCs) in merocyanine (MCYNE) dyes as indicators of polarity. For this purpose, we use Car–Parrinello hybrid quantum mechanics/molecular mechanics (QM/MM) to determine the structures of MCYNE in solvents of different polarity, followed by computations of the SSCCs by using QM/MM linear‐response theory. The molecular geometry of MCYNE switches between neutral, cyanine‐like, and zwitterionic depending on the polarity of the solvent. This structural variation is clearly reflected in the proton SSCCs in the polymethine backbone, which are highly sensitive to the dielectric nature of the environment; this mechanism can be used as a “polarity indicator” for different microenvironments. This result is highlighted by computing the SSCCs of the MCYNE probe in the cavity of the beta‐lactoglobulin protein. The computed SSCCs clearly indicate a non‐polar hydrophobic dielectric nature of this cavity.     

酸性磷化合物萃取金属反应中取代基空间效应的分子力学研究

本文利用分子力学计算从理论上估算了酸性有机磷化合物萃取金属离子反应中取代基的空间效应。二烷基磷(膦)酸叔丁酯(Ⅱ)作为萃合物分子的近似模型结构。计算结果表明, 羟基氧在自由萃取剂Ⅰ和Ⅱ中的局部空间张力能之差与烷基的结构有密切的关系, 可近似地用来作为烷基空间效应的量度(E_(s,ex))。并用于一系列有机磷(膦)酸萃取Ni、Co为稀土的多元回归分析。结果表明, E_(s,ex)基本上反映了该萃取过程中的空间效应。     

Electrostatic interactions in molecular mechanics (MM2) calculations via PEOE partial charges I. Haloalkanes

The PEOE (partial equalization of orbital electronegativity) procedure has been modified slightly and reparametrized for haloalkanes to calculate partial atomic charges suitable for evaluation of dipole moments and electrostatic energies in conjunction with molecular mechanics (MM2) calculations. Dipole moments of 66 haloalkanes are calculated with an average absolute deviation of 0.14 D from experimental values. The conformational energies of 40 compounds have been calculated and the agreement with experimental data is generally good and compares well with calculations by the IDME (induced dipole moment and energy) method. In addition, carbon and proton charges correlate well with C-1s core binding energies and ¹H-NMR (nuclear magnetic resonance) shifts for halomethanes. The most striking benefit of treating electrostatics through a set of partial charges compared to the standard MM2 bond dipole approach is demonstrated by calculations on 1,4-disubstituted cyclohexanes, for which standard MM2 fails to predict the most stable conformation.     

Adsorption energies for a nanoporous carbon from gas-solid chromatography and molecular mechanics

Gas-solid chromatography was used to obtain second gas-solid virial coefficients, B2s, in the temperature range 342-613 K for methane, ethane, propane, butane, 2-methylpropane, chloromethane, chlorodifluoromethane, dichloromethane, and dichlorodifluoromethane. The adsorbent used was Carbosieve S-III (Supelco), a carbon powder with fairly uniform, predominately 0.55 nm slit width pores and a N2 BET surface area of 995 m2/g. The temperature dependence of B2s was used to determine experimental values of the gas-solid interaction energy, E*, for each of these molecular adsorbates. MM2 and MM3 molecular mechanics calculations were used to determine the gas-solid interaction energy, E*(cal), for each of the molecules on various flat and nanoporous model surfaces. The flat model consisted of three parallel graphene layers with each graphene layer containing 127 interconnected benzene rings. The nanoporous model consisted of two sets of three parallel graphene layers adjacent to one another but separated to represent the pore diameter. A variety of calculated adsorption energies, E*(cal), were compared and correlated to the experimental E* values. It was determined that simple molecular mechanics could be used to calculate an attraction energy parameter between an adsorbed molecule and the carbon surface. The best correlation between the E*(cal) and E* values was provided by a 0.50 nm nanoporous model using MM2 parameters.     

The MMP2 calculational method

The molecular mechanics (MMP2) program and procedures for the treatment of conjugated hydrocarbons, and some of the results which they can achieve are described. The program is an updated version of the similar MMP1 program, but contains some differences. It is based on an SCF π system calculation, rather than on the VESCF method used earlier. All parameters are compatible with those in the MM2 program. Hence it is possible to calculate heats of formation, resonance energies, and structures for conjugated hydrocarbons in a way that is consistent with the calculations on non-conjugated molecules. The overall results as far as structure and energy are somewhat better than they were with the MMP1 program.     

An ab initio and molecular mechanical investigation of ureas and amide derivatives

Ab initio molecular orbital theory with the 6-31G* basis set has been used to investigate the geometries and preferred conformations for urea, derivatives of urea, and a few complicated amide derivatives. The results from the ab initio calculations provide insight into the gas-phase rotational barrier about the C? N bond and have been used to generate parameters for the MM2(87) molecular mechanics program. When applicable, theoretical structures are compared with corresponding previously reported experimental geometries. Urea is predicted to be nonplanar with pyramidal amino groups.     

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.     

The generalized molecular fractionation with conjugate caps/molecular mechanics method for direct calculation of protein energy

A generalized molecular fractionation with conjugate caps/molecular mechanics (GMFCC/MM) scheme is developed for efficient linear-scaling quantum mechanical calculation of protein energy. In this GMFCC/MM scheme, the interaction energy between neighboring residues as well as between non-neighboring residues that are spatially in close contact are computed by quantum mechanics while the rest of the interaction energy is computed by molecular mechanics. Numerical studies are carried out to calculate torsional energies of six polypeptides using the GMFCC/MM approach and the energies are shown to be in general good agreement with the full system quantum calculation. Among those we tested is a polypeptide containing 396 atoms whose energies are computed at the MP26-31G* level. Our study shows that using GMFCC/MM, it is possible to perform high level ab initio calculation such as MP2 for applications such as structural optimization of protein complex and molecular dynamics simulation.     

Prediction of the mechanisms of enzyme-catalysed reactions using hybrid quantum mechanical/molecular mechanical methods

The use of hybrid methods, involving both quantum mechanics and molecular mechanics, to model the mechanism of enzyme-catalysed reactions, is discussed. Two alternative approaches to treating the electrostatic interactions between the quantum mechanical and molecular mechanical regions are studied, involving either the inclusion of this term in the electronic Hamiltonian (QM/MM), or evaluating it purely classically (MO + MM). In the latter scheme, possible problems of using force fields that are standard for macromolecular modelling are identified. The use of QM/MM schemes to investigate the mechanism of the enzymes thymidine phosphorylase (ThdPase) and protein tyrosine phosphatase (PTP) is described. For both systems, transition states have been identified using a PM3 Hamiltonian. For ThdPase, concerted motion of the enzyme during the course of the reaction is suggested and, for PTP, a two-step dephosphorylation reaction is indicated, both with quite low barriers.     

GARLEEK: Adding an extra flavor to ONIOM

The ONIOM method, developed in the group of Keiji Morokuma, is one of the most successful examples of quantum mechanics/molecular mechanics (QM/MM) treatments, and of multilayer methods in general. Its implementation in the Gaussian program package is in particular widely used. This implementation has access to the wide variety of QM methods available in Gaussian, but is limited to only three specific force fields. The current article presents the GARLEEK interface, which expands the availability of molecular mechanics methods to the wide variety of force fields available in MM packages. The focus is in the simple installation and use. Two examples of the performance of the interface with selected systems are provided. GARLEEK is MIT-licensed and freely available at https://github.com/insilichem/garleek . © 2018 Wiley Periodicals, Inc.