Integrated Modeling Program, Applied Chemical Theory (IMPACT).

We provide an overview of the IMPACT molecular mechanics program with an emphasis on recent developments and a description of its current functionality. With respect to core molecular mechanics technologies we include a status report for the fixed charge and polarizable force fields that can be used with the program and illustrate how the force fields, when used together with new atom typing and parameter assignment modules, have greatly expanded the coverage of organic compounds and medicinally relevant ligands. As we discuss in this review, explicit solvent simulations have been used to guide our design of implicit solvent models based on the generalized Born framework and a novel nonpolar estimator that have recently been incorporated into the program. With IMPACT it is possible to use several different advanced conformational sampling algorithms based on combining features of molecular dynamics and Monte Carlo simulations. The program includes two specialized molecular mechanics modules: Glide, a high-throughput docking program, and QSite, a mixed quantum mechanics/molecular mechanics module. These modules employ the IMPACT infrastructure as a starting point for the construction of the protein model and assignment of molecular mechanics parameters, but have then been developed to meet specialized objectives with respect to sampling and the energy function.

Pseudospectral Local Second-Order Møller-Plesset Methods for Computation of Hydrogen Bonding Energies of Molecular Pairs.

We present a methodology for computing the binding energy of molecular dimers based on extrapolation of pseudospectral local second-order Moller-Plesset (MP2), or PS-LMP2, energies to the basis set limit. The extrapolation protocol is based on carrying out PS-LMP2 calculations with the Dunning cc-pVTZ (-f) and cc-pVQZ (-g) basis sets and then using a simple two-parameter function to compute the final basis set limit results. The function is parametrized to ultralarge basis set MP2 calculations for 5 molecular pairs taken from the literature and then tested by calculating results for a set of formamide dimers for which such calculations have also been carried out. The results agree to within ca. 0.2 kcal/mol with the conventional MP2 large basis set calculations. A specialized, but relatively simple, protocol is described for eliminating noise due to overcompleteness of the basis set. Timing results are presented for the LMP2 calculations, and comparisons are made with the LMP2 methodology of the QChem program. CPU time required by each of the methods scales as N(3), where N is the number of the basis functions, with the PS-LMP2 approach displaying a 2- to 3-fold advantage in the prefactor. We also discuss one set of test cases for which the PS-LMP2 results disagree with those obtained from an alternative type of MP2 calculation, N-methyl acetamide (NMA) dimers, and show that the results for liquid-state simulations using polarizable parameters derived by fitting to the PS-LMP2 binding energies appear to produce better results when compared with experimental data. The convergence issues associated with the alternative MP2 formulation remain to be investigated.

A Polarizable Force Field and Continuum Solvation Methodology for Modeling of Protein-Ligand Interactions.

A polarizable force field, and associated continuum solvation model, have been developed for the explicit purpose of computing and studying the energetics and structural features of protein binding to the wide range of ligands with potential for medicinal applications. Parameters for the polarizable force field (PFF) are derived from gas-phase ab initio calculations and then utilized for applications in which the protein binding to ligands occurs in aqueous solvents, wherein the charge distributions of proteins and ligands can be dramatically altered. The continuum solvation model is based on a self-consistent reaction field description of solvation, incorporating an analytical gradient, that allows energy minimizations (and, potentially, molecular dynamics simulations) of protein/ligand systems in continuum solvent. This technology includes a nonpolar model describing the cost of cavity formation, and van der Waals interactions, between the continuum solvent and protein/ligand solutes. Tests of the structural accuracy and computational stability of the methodology, and timings for energy minimizations of proteins and protein/ligand systems in the condensed phase, are reported. In addition, the derivation of polarizability, electrostatic, exchange repulsion, and torsion parameters from ab initio data is described, along with the use of experimental solvation energies for determining parameters for the solvation model.