The general effective fragment potential (EFP) method provides model potentials for any molecule that is derived from first principles, with no empirically fitted parameters. The EFP method has been interfaced with most currently used ab initio single-reference and multireference quantum mechanics (QM) methods, ranging from Hartree-Fock and coupled cluster theory to multireference perturbation theory. The most recent innovations in the EFP model have been to make the computationally expensive charge transfer term much more efficient and to interface the general EFP dispersion and exchange repulsion interactions with QM methods. Following a summary of the method and its implementation in generally available computer programs, these most recent new developments are discussed.

Conspectus Natural metalloenzymes are often the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions. However, metalloenzymes are occasionally surprising in their selection of catalytic metals, and in their responses to metal substitution. Indeed, from the isolated standpoint of producing the best catalyst, a chemist designing from first-principles would likely choose a different metal. For example, some enzymes employ a redox active metal where a simple Lewis acid is needed. Such are several hydrolases. In other cases, substitution of a non-native metal leads to radical improvements in reactivity. For example, histone deacetylase 8 naturally operates with Zn2+ in the active site but becomes much more active with Fe2+. For β-lactamases, the replacement of the native Zn2+ with Ni2+ was suggested to lead to higher activity as predicted computationally. There are also intriguing cases, such as Fe2+- and Mn2+-dependent ribonucleotide reductases and W4+- and Mo4+-dependent DMSO reductases, where organisms manage to circumvent the scarcity of one metal (e.g., Fe2+) by creating protein structures that utilize another metal (e.g., Mn2+) for the catalysis of the same reaction. Naturally, even though both metal forms are active, one of the metals is preferred in every-day life, and the other metal variant remains dormant until an emergency strikes in the cell. These examples lead to certain questions. When are catalytic metals selected purely for electronic or structural reasons, implying that enzymatic catalysis is optimized to its maximum? When are metal selections a manifestation of competing evolutionary pressures, where choices are dictated not just by catalytic efficiency but also by other factors in the cell? In other words, how can enzymes be improved as catalysts merely through the use of common biological building blocks available to cells? Addressing these questions is highly relevant to the enzyme design community, where the goal is to prepare maximally efficient quasi-natural enzymes for the catalysis of reactions that interest humankind. Due to competing evolutionary pressures, many natural enzymes may not have evolved to be ideal catalysts and can be improved for the isolated purpose of catalysis in vitro when the competing factors are removed. The goal of this Account is not to cover all the possible stories but rather to highlight how variable enzymatic catalysis can be. We want to bring up possible factors affecting the evolution of enzyme structure, and the large- and intermediate-scale structural and electronic effects that metals can induce in the protein, and most importantly, the opportunities for optimization of these enzymes for catalysis in vitro.

The accurate representation of nitrogen-containing heterocycles is essential for modeling biological systems. In this study, the general effective fragment potential (EFP2) method is used to model dimers of benzene and pyridine, complexes for which high-level theoretical data -including large basis spin-component-scaled second-order perturbation theory (SCS-MP2), symmetry-adapted perturbation theory (SAPT), and coupled cluster with singles, doubles, and perturbative triples (CCSD(T))-are available. An extensive comparison of potential energy curves and components of the interaction energy is presented for sandwich, T-shaped, parallel displaced, and hydrogen-bonded structures of these dimers. EFP2 and CCSD(T) potential energy curves for the sandwich, T-shaped, and hydrogen-bonded dimers have an average root-mean-square deviation (RMSD) of 0.49 kcal/mol; EFP2 and SCS-MP2 curves for the parallel displaced dimers have an average RMSD of 0.52 kcal/mol. Additionally, results are presented from an EFP2 Monte Carlo/simulated annealing (MC/SA) computation to sample the potential energy surface of the benzene-pyridine and pyridine dimers.

A method for calculating the dispersion energy between molecules modeled with the general effective fragment potential (EFP2) method and those modeled using a full quantum mechanics (QM) method, e.g., Hartree-Fock (HF) or second-order perturbation theory, is presented. C 6 dispersion coefficients are calculated for pairs of orbitals using dynamic polarizabilities from the EFP2 portion, and dipole integrals and orbital energies from the QM portion of the system. Dividing by the sixth power of the distance between localized molecular orbital centroids yields the first term in the commonly employed London series expansion. A C 8 term is estimated from the C 6 term to achieve closer agreement with symmetry adapted perturbation theory values. Two damping functions for the dispersion energy are evaluated. By using terms that are already computed during an ordinary HF or EFP2 calculation, the new method enables accurate and extremely rapid evaluation of the dispersioninteraction between EFP2 and QM molecules. The dispersion interaction between quantum mechanics and effective fragment potential molecules Quentin A. Smith, 1 Klaus Ruedenberg, 1 Mark S. Gordon, 1,a) and Lyudmila V. Slipchenko 2,a) A method for calculating the dispersion energy between molecules modeled with the general effective fragment potential (EFP2) method and those modeled using a full quantum mechanics (QM) method, e.g., Hartree-Fock (HF) or second-order perturbation theory, is presented. C 6 dispersion coefficients are calculated for pairs of orbitals using dynamic polarizabilities from the EFP2 portion, and dipole integrals and orbital energies from the QM portion of the system. Dividing by the sixth power of the distance between localized molecular orbital centroids yields the first term in the commonly employed London series expansion. A C 8 term is estimated from the C 6 term to achieve closer agreement with symmetry adapted perturbation theory values. Two damping functions for the dispersion energy are evaluated. By using terms that are already computed during an ordinary HF or EFP2 calculation, the new method enables accurate and extremely rapid evaluation of the dispersion interaction between EFP2 and QM molecules.