RM1: a reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I

Twenty years ago, the landmark AM1 was introduced, and has since had an increasingly wide following among chemists due to its consistently good results and time-tested reliability--being presently available in countless computational quantum chemistry programs. However, semiempirical molecular orbital models still are of limited accuracy and need to be improved if the full potential of new linear scaling techniques, such as MOZYME and LocalSCF, is to be realized. Accordingly, in this article we present RM1 (Recife Model 1): a reparameterization of AM1. As before, the properties used in the parameterization procedure were: heats of formation, dipole moments, ionization potentials and geometric variables (bond lengths and angles). Considering that the vast majority of molecules of importance to life can be assembled by using only six elements: C, H, N, O, P, and S, and that by adding the halogens we can now build most molecules of importance to pharmaceutical research, our training set consisted of 1736 molecules, representative of organic and biochemistry, containing C, H, N, O, P, S, F, Cl, Br, and I atoms. Unlike AM1, and similar to PM3, all RM1 parameters have been optimized. For enthalpies of formation, dipole moments, ionization potentials, and interatomic distances, the average errors in RM1, for the 1736 molecules, are less than those for AM1, PM3, and PM5. Indeed, the average errors in kcal x mol(-1) of the enthalpies of formation for AM1, PM3, and PM5 are 11.15, 7.98, and 6.03, whereas for RM1 this value is 5.77. The errors, in Debye, of the dipole moments for AM1, PM3, PM5, and RM1 are, respectively, 0.37, 0.38, 0.50, and 0.34. Likewise, the respective errors for the ionization potentials, in eV, are 0.60, 0.55, 0.48, and 0.45, and the respective errors, in angstroms, for the interatomic distances are 0.036, 0.029, 0.037, and 0.027. The RM1 average error in bond angles of 6.82 degrees is only slightly higher than the AM1 figure of 5.88 degrees, and both are much smaller than the PM3 and PM5 figures of 6.98 degrees and 9.83 degrees, respectively. Moreover, a known error in PM3 nitrogen charges is corrected in RM1. Therefore, RM1 represents an improvement over AM1 and its similar successor PM3, and is probably very competitive with PM5, which is a somewhat different model, and not fully disclosed. RM1 possesses the same analytical construct and the same number of parameters for each atom as AM1, and, therefore, can be easily implemented in any software that already has AM1, not requiring any change in any line of code, with the sole exception of the values of the parameters themselves.     

AM1-SM2 and PM3-SM3 parameterized SCF solvation models for free energies in aqueous solution

Summary Two new continuum solvation models have been presented recently, and in this paper they are explained and reviewed in detail with further examples. Solvation Model 2 (AM1-SM2) is based on the Austin Model 1 and Solvation Model 3 (PM3-SM3) on the Parameterized Model 3 semiempirical Hamiltonian. In addition to the incorporation of phosphorus parameters, both of these new models address specific deficiencies in the original Solvation Model 1 (AM1-SM1), viz., (1) more accurate account is taken of the hydrophobic effect of hydrocarbons, (2) assignment of heavy-atom surface tensions is based on the presence or absence of bonded hydrogen atoms, and (3) the treatment of specific hydration-shell water molecules is more consistent. The new models offer considerably improved performance compared to AM1-SM1 for neutral molecules and essentially equivalent performance for ions. The solute charges within the Parameterized Model 3 Hamiltonian limit the utility of PM3-SM3 for compounds containing nitrogen and possibly phosphorus. For other systems both AM1-SM2 and PM3-SM3 give realistic results, but AM1-SM2 in general outperforms PM3-SM3. Key features of the models are discussed with respect to alternative approaches.     

Comparison of semiempirical calculations for silicon compounds

A comparison is made of the performance of the MINDO/3, MNDO, AM1, and PM3 methods in calculating the nature of the dimer reconstruction observed on the silicon (100) crystal surface. Based on this case study we conclude that MINDO/3 gives the most realistic results, with PM3 calculations being quite similar but both MNDO and AM1 missing some key features of this system and giving rather unrealistic charge distributions. Hence use of PM3 is recommended for Si containing molecules where a lack of parameters or other restrictions prevent the use of MINDO/3.     

Simulation of liquid water using semiempirical Hamiltonians and the divide and conquer approach

This work examines the ability of semiempirical methods to describe the structure of liquid water. Particularly, the standard AM1 and PM3 methods together with recently developed PM3-PIF and PM3-MAIS parametrizations have been considered. We perform molecular dynamics simulations for a system consisting of 64 or 216 water molecules in a periodic cubic box. The whole system is described quantum mechanically. Calculations with 64 molecules have been carried out using standard SCF techniques whereas calculations with 216 molecules have been done using the divide and conquer approach. This method has also been used in one simulation case with 64 molecules for test purposes. Within this scope, the molecular dynamics program ROAR have been coupled with a linear scaling semiempirical code (DivCon) implemented in a parallel MPI version. The predicted liquid water structure using either AM1 or PM3 is shown to be very poor due to well-known limitations of these methods describing hydrogen bonds. In contrast, PM3-PIF and PM3-MAIS calculations lead to results in reasonably good agreement with experimental data. The best results for the heat of vaporization are obtained with the PM3-PIF method. The average induced dipole moment of the water molecule in the liquid is underestimated by all semiempirical techniques, which seems to be related to the NDDO approximation and to the use of minimal basis sets. A brief discussion on charge-transfer effects in liquid water is also presented.     

Evaluation of PM3, AM1, and MNDO for calculation of higher energy ionization potentials

Higher ionization energies were calculated with PM3, AM1, and MNDO for three series of molecules, representative small molecules, molecules containing heteroatoms, and sterically congested alkenes. Values from PM3, AM1, and MNDO were compared to experimental values. In most instances, the semiempirical calculations correctly predict the ordering of higher ionization energies. In the absence of steric hindrance, MNDO is the method of choice. Within groups of molecules, AM1 performs better on hydrocarbons, especially twisted hydrocarbons, than PM3. PM3 commonly gives sigma orbitals which are too high in energy compared to related pi orbitals. PM3 performed better than AM1 with molecules containing oxygen, but failed to give the correct geometry for hydrogen peroxide.     

Comparison of SCC-DFTB and NDDO-based semiempirical molecular orbital methods for organic molecules

Extensive testing of the SCC-DFTB method has been performed, permitting direct comparison to data available for NDDO-based semiempirical methods. For 34 diverse isomerizations of neutral molecules containing the elements C, H, N, and O, the mean absolute errors (MAE) for the enthalpy changes are 2.7, 3.2, 5.0, 5.1, and 7.2 kcal/mol from PDDG/PM3, B3LYP/6-31G(d), PM3, SCC-DFTB, and AM1, respectively. A more comprehensive test was then performed by computing heats of formation for 622 neutral, closed-shell H, C, N, and O-containing molecules; the MAE of 5.8 kcal/mol for SCC-DFTB is intermediate between AM1 (6.8 kcal/mol) and PM3 (4.4 kcal/mol) and significantly higher than for PDDG/PM3 (3.2 kcal/mol). Similarly, SCC-DFTB is found to be less accurate for heats of formation of ions and radicals; however, it is more accurate for conformational energetics and intermolecular interaction energies, though none of the methods perform well for hydrogen bonds with strengths under ca. 7 kcal/mol. SCC-DFTB and the NDDO methods all reproduce MP2/cc-pVTZ molecular geometries with average errors for bond lengths, bond angles, and dihedral angles of only ca. 0.01 A, 1.5 degrees , and 3 degrees . Testing was also carried out for sulfur containing molecules; SCC-DFTB currently yields much less accurate heats of formation in this case than the NDDO-based methods due to the over-stabilization of molecules containing an SO bond.     

Calculations of molecular vibrational frequencies using semiempirical methods

MINDO/3, MNDO, AM1, and PM3 calculations of molecular vibrational frequencies are reported for 61 molecules. All techniques were applied to both well-behaved and badly behaved systems. Overall, MINDO/3 and MNDO were found to contain rather large errors whereas AM1 and PM3 were relatively accurate. Since no technique does well for all molecules, the technique used should be chosen based on the molecular vibration of interest. In general, AM1 and PM3 together provide fairly accurate results.     

NDDO semiempirical approximations coupled with Green's function technique—a reliable approach for calculating ionization potentials

The ionization potentials of different molecules have been calculated with the outer valence Green's function (OVGF) technique, coupled with semiempirical MNDO, AM1 and PM3 methods. It is found that the OVGF method gives significantly better agreement with the experimental data than do results obtained with semiempirical calculations using Koopman's theorem including a new SAM1 and MNDO/d methods. Of the three semiempirical methods tested (MNDO, AM1, PM3) the OVGF (AM1) method gives the best agreement with experiment.     

Accurate parametrized variational calculations of the molecular electronic polarizability by NDDO‐based methods

The variational method for the calculation of the electronic polarizability of molecules within the NDDO‐based semiempirical MO methods MNDO, AM1, and PM3 was parametrized to improve its accuracy. A training set of 156 compounds was used to fit 34 parameters simultaneously for 12 elements using a simplex optimization. The resulting parameters were tested for a test set of 83 molecules and the calculated polarizabilities compared with the experimental data. For AM1, the RMS deviation between experimental and calculated polarizabilities was reduced from 2.99 (using the original variational treatment) to 0.70 ų for the test set and from 2.81 to 0.40 ų for the training set. MNDO and PM3 gave similar improvements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 17–31, 1999     

Accuracy of free energies of hydration using CM1 and CM3 atomic charges

Absolute free energies of hydration (DeltaGhyd) have been computed for 25 diverse organic molecules using partial atomic charges derived from AM1 and PM3 wave functions via the CM1 and CM3 procedures of Cramer, Truhlar, and coworkers. Comparisons are made with results using charges fit to the electrostatic potential surface (EPS) from ab initio 6-31G* wave functions and from the OPLS-AA force field. OPLS Lennard-Jones parameters for the organic molecules were used together with the TIP4P water model in Monte Carlo simulations with free energy perturbation theory. Absolute free energies of hydration were computed for OPLS united-atom and all-atom methane by annihilating the solutes in water and in the gas phase, and absolute DeltaGhyd values for all other molecules were computed via transformation to one of these references. Optimal charge scaling factors were determined by minimizing the unsigned average error between experimental and calculated hydration free energies. The PM3-based charge models do not lead to lower average errors than obtained with the EPS charges for the subset of 13 molecules in the original study. However, improvement is obtained by scaling the CM1A partial charges by 1.14 and the CM3A charges by 1.15, which leads to average errors of 1.0 and 1.1 kcal/mol for the full set of 25 molecules. The scaled CM1A charges also yield the best results for the hydration of amides including the E/Z free-energy difference for N-methylacetamide in water.     

PDDG/PM3 and PDDG/MNDO: improved semiempirical methods

Two new semiempirical methods employing a Pairwise Distance Directed Gaussian modification have been developed: PDDG/PM3 and PDDG/MNDO; they are easily implemented in existing software, and yield heats of formation for compounds containing C, H, N, and O atoms with significantly improved accuracy over the standard NDDO schemes, PM5, PM3, AM1, and MNDO. The PDDG/PM3 results for heats of formation also show substantial improvement over density functional theory with large basis sets. The PDDG modifications consist of a single function, which is added to the existing pairwise core repulsion functions within PM3 and MNDO, a reparameterized semiempirical parameter set, and modified computation of the energy of formation of a gaseous atom. The PDDG addition introduces functional group information via pairwise atomic interactions using only atom-based parameters. For 622 diverse molecules containing C, H, N, and O atoms, mean absolute errors in calculated heats of formation are reduced from 4.4 to 3.2 kcal/mol and from 8.4 to 5.2 kcal/mol using the PDDG modified versions of PM3 and MNDO over the standard versions, respectively. Several specific problems are overcome, including the relative stability of hydrocarbon isomers, and energetics of small rings and molecules containing multiple heteroatoms. The internal consistency of PDDG energies is also significantly improved, enabling more reliable analysis of isomerization energies and trends across series of molecules; PDDG isomerization energies show significant improvement over B3LYP/6-31G* results. Comparison of heats of formation, ionization potentials, dipole moments, isomer, and conformer energetics, intermolecular interaction energies, activation energies, and molecular geometries from the PDDG techniques is made to experimental data and values from other semiempirical and ab initio methods.     

溴代作用对竹红菌乙素分子性质影响的量子化学研究

用半经验量子化学方法AM1和PM3对竹红菌乙素及其溴代物进行了对比计算,考察了溴代作用对竹红菌乙素分子性质的影响。两种方法所得结果均表明,溴代作用使分子的生成热、前线轨道能级及偶极矩等参数都有所降低,溴代作用也影响了竹红菌乙素分子内氢键的性质,并能使其对光的吸收产生红移。     

Comparison of semiempirical MO methods applied to large molecules

The heats of formation for 19 molecules have been calculated with PM3 and AM1 semiempirical methods. The values obtained have been compared with experimental heats of formation. With PM3 and AM1 the average differences between calculated and experimental heats of formation are 8.45 and 12.34 kcal mol^?1 respectively. There are significant differences when large molecules are considered: this suggests that the parameterization should be done including larger molecules.     

Semi-empirical molecular orbital methods including dispersion corrections for the accurate prediction of the full range of intermolecular interactions in biomolecules

Semi-empirical calculations including an empirical dispersive correction are used to calculate intermolecular interaction energies and structures for a large database containing 156 biologically relevant molecules (hydrogen-bonded DNA base pairs, interstrand base pairs, stacked base pairs and amino acid base pairs) for which MP2 and CCSD(T) complete basis set (CBS) limit estimates of the interaction energies are available. The dispersion corrected semi-empirical methods are parameterised against a small training set of 22 complexes having a range of biologically important non-covalent interactions. For the full molecule set (156 complexes), compared to the high-level ab initio database, the mean unsigned errors of the interaction energies at the corrected semi-empirical level are 1.1 (AM1-D) and 1.2 (PM3-D) kcal mol(-1), being a significant improvement over existing AM1 and PM3 methods (8.6 and 8.2 kcal mol(-1)). Importantly, the new semi-empirical methods are capable of describing the diverse range of biological interactions, most notably stacking interactions, which are poorly described by both current AM1 and PM3 methods and by many DFT functionals. The new methods require no more computer time than existing semi-empirical methods and therefore represent an important advance in the study of important biological interactions.     

Comparative parametric method 5 (PM5) study of trans-stilbene

In this work we report a comparative Austin method 1 (AM1), parametric method 3 (PM3), and parametric method 5 (PM5) studies for trans-stilbene in its ground, excited (singlet and triplet), and ionic (positive and negative polarons and bipolarons) states. We evaluated the accuracy of the recently developed PM5 method. PM5 and AM1 predict a non-planar ground and singlet states for trans-stilbene, while PM3 predicts planar ones, which is in agreement with the available experimental data. In general the PM3 and PM5 bond lengths are superior to AM1 while AM1 bond angles are superior to PM3 and PM5 when compared with available experimental data. The PM5 underestimates the cis–trans isomerization energy and and it is not a quite reliable method for the calculation of relative IP values. The presumed PM5 superior performance against AM1 and PM3 was not observed for the stilbene structures.     

Computational structural determination and energy landscape analysis of the hepatic carcinogen 2-(acetylamino)fluorene

 2-(Acetylamino)fluorene (AAF), a potent mutagen and a prototypical example of the mutagenic aromatic amines, forms covalent adducts to DNA after metabolic activation in the liver. A benchmark study of AAF is presented using a number of the most widely used molecular mechanics and semiempirical computational methods and models. The results are compared to higher-level quantum calculations and to experimentally obtained crystal structures. Hydrogen bonding between AAF molecules in the crystal phase complicates the direct comparison of gas-phase theoretical calculations with experiment, so Hartree–Fock (HF) and Becke–Perdew (BP) density functional theory (DFT) calculations are used as benchmarks for the semiempirical and molecular mechanics results. Systematic conformer searches and dihedral energy landscapes were carried out for AAF using the SYBYL and MMFF94 molecular mechanics force fields; the AM1, PM3 and MNDO semiempirical quantum mechanics methods; HF using the 3-21G*and 6-31G* basis sets; and DFT using the nonlocal BP functional and double numerical polarization basis sets. MMFF94, AM1, HF and DFT calculations all predict the same planar structures, whereas SYBYL, MNDO and PM3 all predict various nonplanar geometries. The AM1 energy landscape is in substantial agreement with HF and DFT predictions; MMFF94 is qualitatively similar to HF and DFT; and the MNDO, PM3 and SYBYL results are qualitatively different from the HF and DFT results and from each other. These results are attributed to deficiencies in MNDO, PM3 and SYBYL. The MNDO, PM3 and SYBYL models may be unreliable for compounds in which an amide group is immediately adjacent to an aromatic ring. Received: 26 May 2002 / Accepted: 12 December 2002 / Published online: 14 February 2003     

Are AM1 ligand‐protein binding enthalpies good enough for use in the rational design of new drugs?

We have examined the performance of semiempirical quantum mechanical methods in solving the problem of accurately predicting protein-ligand binding energies and geometries. Firstly, AM1 and PM3 geometries and binding enthalpies between small molecules that simulate typical ligand-protein interactions were compared with high level quantum mechanical techniques that include electronic correlation (e.g., MP2 or B3LYP). Species studied include alkanes, aromatic systems, molecules including groups with hypervalent sulfur or with donor or acceptor hydrogen bonding capability, as well as ammonium or carboxylate ions. B3LYP/6-311+G(2d,p) binding energies correlated very well with the BSSE corrected MP2/6-31G(d) values. AM1 binding enthalpies also showed good correlation with MP2 values, and their systematic deviation is acceptable when enthalpies are used for the comparison of interaction energies between ligands and a target. PM3 otherwise gave erratic energy differences in comparison to the B3LYP or MP2 approaches. As one would expect, the geometries of the binding complexes showed the known limitations of the semiempirical and DFT methods. AM1 calculations were subsequently applied to a test set consisting of "real" protein active site-ligand complexes. Preliminary results indicate that AM1 could be a valuable tool for the design of new drugs using proteins as templates. This approach also has a reasonable computational cost. The ligand-protein X-ray structures were reasonably reproduced by AM1 calculations and the corresponding AM1 binding enthalpies are in agreement with the results from the "small molecules" test set.     

PM3 study of the proton affinities of 2-, 3-, and 4-monosubstituted pyridines in the gas phase

Heats of formation (ΔH_f) and proton affinities (PA) of 2-, 3-, and 4-monosubstituted pyridines in the gas phase are calculated using the AM1 and PM3 semiempirical methods. The following substitutents are considered: NO₂, CN, CF₃, CHO, F, Cl, COCH₃, H, CH₃, OCH₃, SCH₃, NH₂, and N(CH₃)₂. The results are compared with the experimental data. Both methods reproduce the ΔH_f with comparble accuracy; the rms deviations are 4.1 (AM1) and 4.5 kcal/mol (PM3) for the free bases and 9.5 (AM1) and 9.7 kcal/mol (PM3) for their conjugated acids. The PA are systematically underestimated by both methods, but AM1 appears to be clearly better than PM3 for reproducing the experimental values. The rms deviations for AM1 and PM3 are 5.1 and 9.6 kcal/mol, respectively. This is due to a cancellation of systematic errors in the calculated ΔH_f in the AM1 case and to a summation of the errors in the PM3 case. Both methods correctly reproduce conformations of the molecules under consideration.     

Conformation and reactivity of α-oxo-ketenes: Ab initio and semiempirical (AM1, PM3) calculations

Ab initio MP2/6-31G*//MP2/6-31G* and semiempirical AM1 and PM3 calculations on a series of differently substituted α-oxo-ketenes are used to investigate E/Z-isomerism and rotational barriers in these molecules. Sterically crowded derivatives are found to exist solely as s-E conformers. The unusual stability of these derivatives thus can be attributed to their inability to adopt the s-Z conformation required for the normal α-oxo-ketene reactions. With respect to structures and energies, the PM3 method (especially in the case of highly crowded molecules) is found to be less reliable than AM1. Ab initio HF/3-21G and PM3 vibrational frequencies appear to be of sufficient accuracy for a distinction between s-Z and s-E conformers. In this respect, the AM1 method appears less reliable. © 1994 by John Wiley & Sons, Inc.     

Can semi-empirical models describe HCl dissociation in water?

We discuss the failure of commonly used AM1 and PM3 semiempirical methods to correctly describe acid dissociation. We focus our analysis on HCl because of its physicochemical importance and its relevance in atmospheric chemistry. The structure of non-dissociated and dissociated HCl – (H₂O) _n clusters is accounted for. The very bad results obtained with PM3 (and also with AM1) are related to large errors in gas-phase proton affinity of water and gas-phase acidity of HCl. Indeed, estimation of pKa values shows that neither AM1 nor PM3 are able to predict HCl dissociation in liquid water since HCl is found to be a weaker acid than H₃O⁺. We have proposed in previous works a modified PM3 approach (PM3-MAIS) adapted to intermolecular calculations. It is derived from PM3 by reparameterization of the core–core functions using ab initio data. Since parameters for H–Cl and O–Cl core–core interactions were not yet available, we have carried out the corresponding optimization. Application of the PM3-MAIS method to HCl dissociation in HCl–(H₂O) _n clusters leads to a huge improvement over PM3 results. Though the method predicts a slightly overestimated HCl acidity in water environment, the overall agreement with ab initio calculations is very satisfying and justifies efforts to develop new semiempirical methods.