Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases

2013 ◽  
Vol 91 (7) ◽  
pp. 559-572 ◽  
Author(s):  
Jennifer L. Kellie ◽  
Stacey D. Wetmore
Keyword(s):  
Active Sites ◽  
Noncovalent Interactions ◽  
Single Point ◽  
Enzymatic Reaction ◽  
Basis Set ◽  
Basis Sets ◽  
Dna Glycosylases ◽  
Dft Methods ◽  
Primary Focus ◽  
High Level

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.

The X40x10 Halogen Bonding Benchmark Revisited: Surprising Importance of (n-1)d Subvalence Correlation

2018 ◽  
Author(s):  
Manoj Kumar Kesharwani ◽  
Nitai Sylvetsky ◽  
Debashree Manna ◽  
Jan M.L. Martin
Keyword(s):  
Noncovalent Interactions ◽  
Halogen Bonding ◽  
Basis Set ◽  
Basis Sets ◽  
Coupled Cluster ◽  
Dissociation Energies ◽  
Explicitly Correlated ◽  
Two Factors ◽  
High Level ◽  
Cluster Methods

<p>We have re-evaluated the X40x10 benchmark for halogen bonding using conventional and explicitly correlated coupled cluster methods. For the aromatic dimers at small separation, improved CCSD(T)–MP2 “high-level corrections” (HLCs) cause substantial reductions in the dissociation energy. For the bromine and iodine species, (n-1)d subvalence correlation increases dissociation energies, and turns out to be more important for noncovalent interactions than is generally realized; ; (n-1)sp subvalence correlation is much less important. The (n-1)d subvalence term is dominated by core-valence correlation; with the smaller cc-pVDZ-F12-PP and cc-pVTZ-F12-PP basis sets, basis set convergence for the core-core contribution becomes sufficiently erratic that it may compromise results overall. The two factors conspire to generate discrepancies of up to 0.9 kcal/mol (0.16 kcal/mol RMS) between the original X40x10 data and the present revision.</p>

The X40x10 Halogen Bonding Benchmark Revisited: Surprising Importance of (n-1)d Subvalence Correlation

2018 ◽  
Author(s):  
Manoj Kumar Kesharwani ◽  
Nitai Sylvetsky ◽  
Debashree Manna ◽  
Jan M.L. Martin
Keyword(s):  
Noncovalent Interactions ◽  
Halogen Bonding ◽  
Basis Set ◽  
Basis Sets ◽  
Coupled Cluster ◽  
Dissociation Energies ◽  
Explicitly Correlated ◽  
Two Factors ◽  
High Level ◽  
Cluster Methods

<p>We have re-evaluated the X40x10 benchmark for halogen bonding using conventional and explicitly correlated coupled cluster methods. For the aromatic dimers at small separation, improved CCSD(T)–MP2 “high-level corrections” (HLCs) cause substantial reductions in the dissociation energy. For the bromine and iodine species, (n-1)d subvalence correlation increases dissociation energies, and turns out to be more important for noncovalent interactions than is generally realized; ; (n-1)sp subvalence correlation is much less important. The (n-1)d subvalence term is dominated by core-valence correlation; with the smaller cc-pVDZ-F12-PP and cc-pVTZ-F12-PP basis sets, basis set convergence for the core-core contribution becomes sufficiently erratic that it may compromise results overall. The two factors conspire to generate discrepancies of up to 0.9 kcal/mol (0.16 kcal/mol RMS) between the original X40x10 data and the present revision.</p>

scholarly journals The S66 noncovalent interactions benchmark reconsidered using explicitly correlated methods near the basis set limit

2017 ◽  
Author(s):  
Manoj Kumar Kesharwani ◽  
Amir Karton ◽  
Nitai Sylvetsky ◽  
Jan M. L. Martin
Keyword(s):  
Noncovalent Interactions ◽  
Computational Cost ◽  
The Other ◽  
Basis Set ◽  
Basis Sets ◽  
Other Hand ◽  
Explicitly Correlated ◽  
Explicitly Correlated Methods ◽  
The One ◽  
High Level

<p>The S66 benchmark for noncovalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 “high-level corrections” are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66x8 problem set in particular, and for accurate calculations on noncovalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, while half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.</p>

scholarly journals The S66 noncovalent interactions benchmark reconsidered using explicitly correlated methods near the basis set limit

2017 ◽  
Author(s):  
Manoj Kumar Kesharwani ◽  
Amir Karton ◽  
Nitai Sylvetsky ◽  
Jan M. L. Martin
Keyword(s):  
Noncovalent Interactions ◽  
Computational Cost ◽  
The Other ◽  
Basis Set ◽  
Basis Sets ◽  
Final Version ◽  
Explicitly Correlated ◽  
Explicitly Correlated Methods ◽  
The One ◽  
High Level

<p>[Final version of record available at http://dx.doi.org/10.1071/CH17588]</p><p>The S66 benchmark for noncovalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 “high-level corrections” are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66x8 problem set in particular, and for accurate calculations on noncovalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, while half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.<br></p>

scholarly journals The S66 noncovalent interactions benchmark reconsidered using explicitly correlated methods near the basis set limit

2017 ◽  
Author(s):  
Manoj Kumar Kesharwani ◽  
Amir Karton ◽  
Nitai Sylvetsky ◽  
Jan M. L. Martin
Keyword(s):  
Noncovalent Interactions ◽  
Computational Cost ◽  
The Other ◽  
Basis Set ◽  
Basis Sets ◽  
Other Hand ◽  
Explicitly Correlated ◽  
Explicitly Correlated Methods ◽  
The One ◽  
High Level

<p>The S66 benchmark for noncovalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 “high-level corrections” are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66x8 problem set in particular, and for accurate calculations on noncovalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, while half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.<br></p>

A Systematic DFT Study of Two Elementary Steps Involving Hydrogen Peroxide and the Hydroxyl Radical in the Fenton Reaction

2018 ◽  
Author(s):  
Danilo Carmona ◽  
David Contreras ◽  
Oscar A. Douglas-Gallardo ◽  
Stefan Vogt-Geisse ◽  
Pablo Jaque ◽  
...  
Keyword(s):  
Hydrogen Peroxide ◽  
Hydroxyl Radical ◽  
Ab Initio ◽  
Fenton Reaction ◽  
Basis Set ◽  
Basis Sets ◽  
Dft Methods ◽  
Elementary Steps ◽  
Oxygen Oxygen ◽  
Complete Basis Set

The Fenton reaction plays a central role in many chemical and biological processes and has various applications as e.g. water remediation. The reaction consists of the iron-catalyzed homolytic cleavage of the oxygen-oxygen bond in the hydrogen peroxide molecule and the reduction of the hydroxyl radical. Here, we study these two elementary steps with high-level ab-initio calculations at the complete basis set limit and address the performance of different DFT methods following a specific classification based on the Jacob´s ladder in combination with various Pople's basis sets. Ab-initio calculations at the complete basis set limit are in agreement to experimental reference data and identified a significant contribution of the electron correlation energy to the bond dissociation energy (BDE) of the oxygen-oxygen bond in hydrogen peroxide and the electron affinity (EA) of the hydroxyl radical. The studied DFT methods were able to reproduce the ab-initio reference values, although no functional was particularly better for both reactions. The inclusion of HF exchange in the DFT functionals lead in most cases to larger deviations, which might be related to the poor description of the two reactions by the HF method. Considering the computational cost, DFT methods provide better BDE and EA values than HF and post--HF methods with an almost MP2 or CCSD level of accuracy. However, no systematic general prediction of the error based on the employed functional could be established and no systematic improvement with increasing the size in the Pople's basis set was found, although for BDE values certain systematic basis set dependence was observed. Moreover, the quality of the hydrogen peroxide, hydroxyl radical and hydroxyl anion structures obtained from these functionals was compared to experimental reference data. In general, bond lengths were well reproduced and the error in the angles were between one and two degrees with some systematic trend with the basis sets. From our results we conclude that DFT methods present a computationally less expensive alternative to describe the two elementary steps of the Fenton reaction. However, choice of approximated functionals and basis sets must be carefully done and the provided benchmark allows a systematic validation of the electronic structure method to be employed

DFT Benchmark Study of the O--O Bond Dissociation Energy in Peroxides Validated with High-Level Ab-Initio Calculations

2019 ◽  
Author(s):  
Danilo Carmona ◽  
Pablo Jaque ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Reference Value ◽  
Basis Set ◽  
Basis Sets ◽  
Mean Deviation ◽  
Bond Dissociation ◽  
Dissociation Energies ◽  
The Mean ◽  
Experimental Values ◽  
High Level ◽  
Almost All

<div><div><div><p>Peroxides play a central role in many chemical and biological pro- cesses such as the Fenton reaction. The relevance of these compounds lies in the low stability of the O–O bond which upon dissociation results in radical species able to initiate various chemical or biological processes. In this work, a set of 64 DFT functional-basis set combinations has been validated in terms of their capability to describe bond dissociation energies (BDE) for the O–O bond in a database of 14 ROOH peroxides for which experimental values ofBDE are available. Moreover, the electronic contributions to the BDE were obtained for four of the peroxides and the anion H2O2− at the CBS limit at CCSD(T) level with Dunning’s basis sets up to triple–ζ quality provid- ing a reference value for the hydrogen peroxide anion as a model. Almost all the functionals considered here yielded mean absolute deviations around 5.0 kcal mol−1. The smallest values were observed for the ωB97 family and the Minnesota M11 functional with a marked basis set dependence. Despite the mean deviation, order relations among BDE experimental values of peroxides were also considered. The ωB97 family was able to reproduce the relations correctly whereas other functionals presented a marked dependence on the chemical nature of the R group. Interestingly, M11 functional did not show a very good agreement with the established order despite its good performance in the mean error. The obtained results support the use of similar validation strategies for proper prediction of BDE or other molecular properties by DF Tmethods in subsequent related studies.</p></div></div></div>

scholarly journals The S66 Non-Covalent Interactions Benchmark Reconsidered Using Explicitly Correlated Methods Near the Basis Set Limit

2018 ◽  
Vol 71 (4) ◽  
pp. 238 ◽  
Author(s):  
Manoj K. Kesharwani ◽  
Amir Karton ◽  
Nitai Sylvetsky ◽  
Jan M. L. Martin
Keyword(s):  
Computational Cost ◽  
The Other ◽  
Basis Set ◽  
Basis Sets ◽  
Explicitly Correlated ◽  
Non Covalent Interactions ◽  
Explicitly Correlated Methods ◽  
The One ◽  
High Level ◽  
Covalent Interactions

The S66 benchmark for non-covalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 ‘high-level corrections’ are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66×8 problem set in particular, and for accurate calculations on non-covalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, whereas half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.

Theoretical study on NO3-initiated oxidation of acenaphthene in the atmosphere

2008 ◽  
Vol 86 (2) ◽  
pp. 129-137 ◽  
Author(s):  
Xiaohui Qu ◽  
Qingzhu Zhang ◽  
Wenxing Wang
Keyword(s):  
Product Information ◽  
Environmental Concern ◽  
Single Point ◽  
Basis Set ◽  
Orbital Theory ◽  
Night Time ◽  
Point Energy ◽  
Secondary Reactions ◽  
Favorable Reaction ◽  
High Level

Acenaphthene is widespread and toxic, and thus of substantial environmental concern. The reaction with NO3 radicals is an important atmospheric loss process of acenaphthene at night time. In this work, the mechanism for the NO3-initiated atmospheric oxidation reaction of acenaphthene has been studied using high level molecular orbital theory. Geometries of all the related species have been optimized at the MPWB1K level with the 6–31G(d,p) basis set. The single-point energy calculations have been carried out at the MPWB1K/6–311+G(3df,2p) level. The possible secondary reactions were also studied. Several energetically favorable reaction pathways were revealed for the first time.Key words: acenaphthene, NO3 radicals, reaction mechanism, product information, oxidation degradation.

scholarly journals Basis Set Effects in the Description of the Cl-O Bond in ClO and XClO/ClOX Isomers (X = H, O, and Cl) Using DFT and CCSD(T) Methods

2019 ◽  
Vol 2019 ◽  
pp. 1-23 ◽  
Author(s):  
Kenneth Irving ◽  
Martina Kieninger ◽  
Oscar N. Ventura
Keyword(s):  
Density Functional ◽  
Basis Set ◽  
Basis Sets ◽  
Enthalpies Of Formation ◽  
Isodesmic Reactions ◽  
Dft Methods ◽  
Functional Methods ◽  
The One ◽  
Isomeric Structures ◽  
Different Sources

The performance of a group of density functional methods of progressive complexity for the description of the ClO bond in a series of chlorine oxides was investigated. The simplest ClO radical species and the two isomeric structures XClO/ClOX for each X = H, Cl, and O were studied using the PW91, TPSS, B3LYP, PBE0, M06, M06-2X, BMK, and B2PLYP functionals. Geometry optimizations and reaction enthalpies and enthalpies of formation for each species were calculated using Pople basis sets and the (aug)-cc-pVnZ Dunning sets, with n = D, T, Q, 5, and 6. For the calculation of enthalpies of formation, atomization and isodesmic reactions were employed. Both the precision of the methods with respect to the increase of the basis sets, as well as their accuracy, were gauged by comparing the results with the more accurate CCSD(T) calculations, performed using the same basis sets as for the DFT methods. The results obtained employing composite chemical methods (G4, CBS-QB3, and W1BD) were also used for the comparisons, as well as the experimental results when they are available. The results obtained show that error compensation is the key for successful description of molecular properties (geometries and energies) by carefully selecting the method and basis sets. In general, expansion of the one-electron basis set to the limit of completeness does not improve results at the DFT level, but just the opposite. The enthalpies of formation calculated at the CCSD(T)/aug-cc-pV6Z for the species considered are generally in agreement with experimental determinations and the most accurate theoretical values. Different sources of error in the calculations are discussed in detail.

Sign in / Sign up

Export Citation Format

Share Document