Basis set convergence studies of Hartree–Fock calculations of molecular properties within the resolution of the identity approximation

Artemyev, A.I.; Bibikov, A. V.; Zayets, Valentin A.; Bodrenko, Igor

doi:10.1063/1.1947193

Cited by 26 publications

(12 citation statements)

References 22 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Dimer calculations were performed with the ALGOQMT software, 18 which exploits the very efficient resolution of the identity approximate integral technique. Molecular geometry was optimized by the PC GAMESS ͑Ref.…”

Section: Methodsmentioning

confidence: 99%

Ab initio study of the effects of orientation and corrugation for H2 adsorbed on polycyclic aromatic hydrocarbons

Donchev¹

2007

The Journal of Chemical Physics

View full text Add to dashboard Cite

This paper reports state-of-the-art ab initio calculations at the second order of Moller-Plesset perturbation theory of molecular hydrogen interaction with polycyclic aromatic hydrocarbons (PAHs) up to coronene (C(24)H(12)). The effects of both H(2) orientation with respect to the PAH plane and of PAH corrugation have been carefully investigated. It was found that the energetic preference for the perpendicular H(2) orientation over the parallel one decreases rapidly as the PAH grows. The detailed study of coronene-H(2) potential energy surface has made it possible to estimate from first principles the graphite-H(2) binding energy. This ab initio estimate is shown to be in fair agreement with experiment. As revealed by the example of benzene-H(2) dimer, the electrostatic energy component plays an important role in PAH-H(2) interactions, that stems largely from the charge penetration effect and, therefore, cannot be reproduced by the simple quadrupole-quadrupole interaction model.

show abstract

Section: Methodsmentioning

confidence: 99%

Ab initio study of the effects of orientation and corrugation for H2 adsorbed on polycyclic aromatic hydrocarbons

Donchev¹

2007

The Journal of Chemical Physics

View full text Add to dashboard Cite

show abstract

“…The results are especially promising because further progress is possible and its direction is clear: extremely accurate ab initio calculations can be performed with advanced methods, e.g., ref. 23, followed by further refinement of the functional forms and parameterization of the QMPFF interaction potentials to track the ab initio data.…”

Section: Discussionmentioning

confidence: 99%

Water properties from first principles: Simulations by a general-purpose quantum mechanical polarizable force field

Donchev¹,

Галкин²,

Illarionov³

et al. 2006

Proc. Natl. Acad. Sci. U.S.A.

View full text Add to dashboard Cite

We have recently introduced a quantum mechanical polarizable force field (QMPFF) fitted solely to high-level quantum mechanical data for simulations of biomolecular systems. Here, we present an improved form of the force field, QMPFF2, and apply it to simulations of liquid water. The results of the simulations show excellent agreement with a variety of experimental thermodynamic and structural data, as good or better than that provided by specialized water potentials. In particular, QMPFF2 is the only ab initio force field to accurately reproduce the anomalous temperature dependence of water density to our knowledge. The ability of the same force field to successfully simulate the properties of both organic molecules and water suggests it will be useful for simulations of proteins and protein-ligand interactions in the aqueous environment.G eneral-purpose force fields, from Levitt's early protein potential (1) to modern models such as CHARMM, OPLS-AA, MMFF, and AMBER (2-5), which approximate molecular potentials by simple analytical formulas, are in wide use for computational studies of biological systems ranging from the simplest molecular clusters to large complexes involving proteins. In the latter case, the investigations encounter serious computational problems, primarily related to proper conformational sampling and adequate treatment of the long-range intermolecular interactions; however, with advancements in simulation methodologies and the increase in computer speed these difficulties are alleviated so the accuracy of the underlying models becomes the dominant factor.Protein and protein-ligand interactions usually take place in an aqueous environment, which contributes critically to their energetics, e.g., by hydrogen bonding and the hydrophobic effect. Hence, a force field should accurately reproduce the properties of both organic compounds and water if it is to be used for precise calculations of protein-ligand binding, as required for example in drug-design applications. Moreover, the quality of the applications of a force field to water can be considered as a criterion for the accuracy of the approach as a whole. Hence, it is disconcerting that no general-purpose force field has previously succeeded in accurately describing key properties of liquid water.On the other hand, impressive progress has been made in theoretical studies using specialized water potentials. Many of these potentials are empirical, i.e., they have been fitted to experimental data on the thermodynamics and kinetics of liquid water and in some cases ice. The most advanced of these models, such as the pairwise additive TIP5P (6) and polarizable (7-9) potentials, generally provide an accurate description of the most important properties of water and͞or ice. However, no one model is yet able to reproduce in detail the diversity of thermodynamic and kinetic experimental data on both gas and condensed phases under a range of conditions. Moreover, these empirical water potentials cannot be transferred to more general molecular systems ...

show abstract

“…In particular, we find that the QMPFF3 parameters underestimate the barrier for coplanar positioning of the amidine group (terminal NH 2 -C-NH 2 ) with the adjacent aromatic ring by Ϸ1-1.5 kcal/mol in comparison with the ab initio quantum value. [The latter was calculated by using our in-house AlgoQMT software (43,44) at the MP2/TZ* level, the same level as in QM calculations used for the QMPFF3 parameterization in ref. 36.]…”

Section: Resultsmentioning

confidence: 99%

Application of a polarizable force field to calculations of relative protein–ligand binding affinities

Khoruzhii¹,

Donchev²,

Галкин³

et al. 2008

Proc. Natl. Acad. Sci. U.S.A.

View full text Add to dashboard Cite

An explicitly polarizable force field based exclusively on quantum data is applied to calculations of relative binding affinities of ligands to proteins. Five ligands, differing by replacement of an atom or functional group, in complexes with three serine proteases-trypsin, thrombin, and urokinase-type plasminogen activator-with available experimental binding data are used as test systems. A special protocol of thermodynamic integration was developed and used to provide sufficiently low levels of systematic error along with high numerical efficiency and statistical stability. The calculated results are in excellent quantitative (rmsd ‫؍‬ 1.0 kcal/mol) and qualitative (R 2 ‫؍‬ 0.90) agreement with experimental data. The potential of the methodology to explain the observed differences in the ligand affinities is also demonstrated. molecular dynamics simulation ͉ serine protease ͉ drug design T he ability to accurately calculate the binding affinity, or equivalently binding free energy, of a ligand for a protein would be highly useful in the field of drug design for lead selection and optimization. Although screening and docking methods (1, 2) have been successful in filtering large chemical databases, they cannot provide definitive calculations of binding energy because of the simplified scoring functions used and the restricted number of states tested. Hence, more accurate calculations have generally been based on molecular mechanics models (3, 4). In these methods the required properties of statistical ensembles are determined by molecular dynamical or Monte Carlo simulations of systems by using a physically grounded description of interactions between particles. The theoretical thermodynamic foundation of such methods is clear, simple, and well established (5).Despite these advantages, as well as initial optimism and long development, only a limited number of successful simulations have been published during the past decade (e.g., refs. 6-14; for review of early results, see refs. 3 and 4). In part, this is explained by many methodological difficulties that must be overcome, especially in accurate and efficient description of long-range interactions and adequate sampling of the conformational space. Major efforts devoted to solving these problems have resulted both in partial success (15, 16) and the recognition of some principal limitations. Notably, several theoretically based techniques have been developed, such as umbrella sampling, the concept of potential of mean force, and artificial restraining potentials, which restrict or decompose the conformational space and simplify adequate sampling (e.g., 11, 17-20).The other principal problem has been the quality of the model potentials or force fields (FFs) used to describe atom-atom interactions. The most questionable point in this respect is the role of nonadditive effects, particularly electronic polarizability. Widely used FFs such as MMFF (21), AMBER (22), OPLS (23), CHARMM (24), and GROMOS (25) are not explicitly polarizable but rather include polarizabil...

show abstract

Basis set convergence studies of Hartree–Fock calculations of molecular properties within the resolution of the identity approximation

Cited by 26 publications

References 22 publications

Ab initio study of the effects of orientation and corrugation for H2 adsorbed on polycyclic aromatic hydrocarbons

Ab initio study of the effects of orientation and corrugation for H2 adsorbed on polycyclic aromatic hydrocarbons

Water properties from first principles: Simulations by a general-purpose quantum mechanical polarizable force field

Application of a polarizable force field to calculations of relative protein–ligand binding affinities

Contact Info

Product

Resources

About