Noncovalent Interactions from Models for the Moller-Plesset Adiabatic Connection.
Ontology highlight
ABSTRACT:
SUBMITTER: Daas TJ
PROVIDER: S-EPMC8280728 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
| S-EPMC9948597 | biostudies-literature
Project description:The performances of Møller-Plesset second-order perturbation theory (MP2) and density functional theory (DFT) have been assessed for the purposes of investigating the interaction between stannylenes and aromatic molecules. The complexes between SnX2 (where X = H, F, Cl, Br, and I) and benzene or pyridine are considered. Structural and energetic properties of such complexes are calculated using six MP2-type and 14 DFT methods. The assessment of the above-mentioned methods is based on the comparison of the structures and interaction energies predicted by these methods with reference computational data. A very detailed analysis of the performances of the MP2-type and DFT methods is carried out for two complexes, namely SnH2-benzene and SnH2-pyridine. Of the MP2-type methods, the reference structure of SnH2-benzene is reproduced best by SOS-MP2, whereas SCS-MP2 is capable of mimicking the reference structure of SnH2-pyridine with the greatest accuracy. The latter method performs best in predicting the interaction energy between SnH2 and benzene or pyridine. Among the DFT methods, ωB97X provides the structures and interaction energies of the SnH2-benzene and SnH2-pyridine complexes with good accuracy. However, this density functional is not as effective in reproducing the reference data for the two complexes as the best performing MP2-type methods. Next, the DFT methods are evaluated using the full test set of SnX2-benzene and SnX2-pyridine complexes. It is found that the range-separated hybrid or dispersion-corrected density functionals should be used for describing the interaction in such complexes with reasonable accuracy.
| S-EPMC4326664 | biostudies-literature
Project description:Molecular structure does not easily identify the intricate noncovalent interactions that govern many areas of biology and chemistry, including design of new materials and drugs. We develop an approach to detect noncovalent interactions in real space, based on the electron density and its derivatives. Our approach reveals the underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of nonbonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.
| S-EPMC2864795 | biostudies-literature
Project description:Approximate exchange-correlation functionals built by modeling in a nonlinear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameter-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular systems. We show that size consistency in the AC-based functionals can be restored in a very simple way at no extra computational cost. Results on a large set of benchmark molecular interaction energies show that functionals based on the interaction strength interpolation approximations are significantly more accurate than second-order perturbation theory.
| S-EPMC5994725 | biostudies-literature
Project description:Insulin degrading enzyme (IDE), a metalloprotease that degrades amyloid-β (Aβ) peptides and insulin, is associated with Alzheimer's disease and diabetes. The mechanism of IDE catalyzed degrading of Aβ peptides, which is of fundamental importance in the design of therapeutic methods for Alzheimer's disease, has not been fully understood. In this work, combined quantum mechanics and molecular mechanics (QM/MM) style Møller-Plesset second order perturbation theory (MP2) geometry optimization calculations are performed to investigate the catalytic mechanism of the Aβ40 Phe19-Phe20 peptide bond cleavage by human IDE. The analyses using QM/MM MP2 optimization suggest that a neutral water molecule is at the active site of the enzyme-substrate (ES) complex. The water molecule is in hydrogen bonding with the nearby anionic Glu111 of IDE but not directly bound to the catalytic Zn ion. This is confirmed by QM/MM DFTB3 molecular dynamics simulation. Our studies also reveal that the hydrolysis of the Aβ40 Phe19-Phe20 peptide bond by IDE consists of four key steps. The neutral water is first activated by moving toward and binding to the Zn ion. A gem-diol intermediate is then formed by the activated neutral water molecule attacking the C atom of the Phe19-Phe20 peptide bond. The next is the protonation of the N atom of Phe19-Phe20 peptide bond to form an intermediate with an elongated C-N bond. The final step is the breaking of the Phe19-Phe20 C-N bond. The final step is the rate-determining step with a calculated Gibbs free energy of activation of 17.34 kcal/mol, in good agreement with the experimental value 16.7 kcal/mol. This mechanism provides the basis for the design of biochemical methods to modulate the activity of IDE in humans.
| S-EPMC6670292 | biostudies-literature
Project description:The construction of density-functional approximations is explored by modeling the adiabatic connection locally, using energy densities defined in terms of the electrostatic potential of the exchange-correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to assess the accuracy of a range of local interpolation models against accurate exchange-correlation energy densities. The importance of the strictly correlated electrons (SCE) functional describing the strong coupling limit is emphasized, enabling the corresponding interpolated functionals to treat strong correlation effects. In addition to exploring the performance of such models numerically for the helium and beryllium isoelectronic series and the dissociation of the hydrogen molecule, an approximate analytic model is presented for the initial slope of the local adiabatic connection. Comparisons are made with approaches based on global models, and prospects for future approximations based on the local adiabatic connection are discussed.
| S-EPMC4910137 | biostudies-other
Project description:The electron-nucleus hyperfine coupling constant is a challenging property for density functional methods. For accurate results, hybrid functionals with a large amount of exact exchange are often needed and there is no clear "one-for-all" functional which describes the hyperfine coupling interaction for a large set of nuclei. To alleviate this unfavorable situation, we apply the adiabatic connection random phase approximation (RPA) in its post-Kohn-Sham fashion to this property as a first test. For simplicity, only the Fermi-contact and spin-dipole terms are calculated within the nonrelativistic and the scalar-relativistic exact two-component framework. This requires to solve a single coupled-perturbed Kohn-Sham equation to evaluate the relaxed density matrix, which comes with a modest increase in computational demands. RPA performs remarkably well and substantially improves upon its Kohn-Sham (KS) starting point while also reducing the dependence on the KS reference. For main-group systems, RPA outperforms global, range-separated, and local hybrid functionals─at similar computational costs. For transition-metal compounds and lanthanide complexes, a similar performance as for hybrid functionals is observed. In contrast, related post-Hartree-Fock methods such as Møller-Plesset perturbation theory or CC2 perform worse than semilocal density functionals.
| S-EPMC11372758 | biostudies-literature
Project description:There is a biomedical need to develop molecular recognition systems that selectively target the interfaces of protein and lipid aggregates in biomembranes. This is an extremely challenging problem in supramolecular chemistry because the biological membrane is a complex dynamic assembly of multifarious molecular components with local inhomogeneity. Two simplifying concepts are presented as a framework for basing molecular design strategies. The first generalization is that association of two binding partners in a biomembrane will be dominated by one type of non-covalent interaction which is referred to as the keystone interaction. Structural mutations in membrane proteins that alter the strength of this keystone interaction will likely have a major effect on biological activity and often will be associated with disease. The second generalization is to view the structure of a cell membrane as three spatial regions, that is, the polar membrane surface, the midpolar interfacial region and the non-polar membrane interior. Each region has a distinct dielectric, and the dominating keystone interaction between binding partners will be different. At the highly polar membrane surface, the keystone interactions between charged binding partners are ion-ion and ion-dipole interactions; whereas, ion-dipole and ionic hydrogen bonding are very influential at the mid-polar interfacial region. In the non-polar membrane interior, van der Waals forces and neutral hydrogen bonding are the keystone interactions that often drive molecular association. Selected examples of lipid and transmembrane protein association systems are described to illustrate how the association thermodynamics and kinetics are dominated by these keystone noncovalent interactions.
| S-EPMC2857542 | biostudies-literature
Project description:Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The ACn method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, ACn is more efficient than existing ab initio multireference dynamic correlation methods. ACn affords excellent results for singlet-triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.
| S-EPMC9150121 | biostudies-literature
Project description:One of the most familiar carbon-centered noncovalent interactions (NCIs) involving an antibonding π*-orbital situated at the Bürgi-Dunitz angle from the electron donor, mostly lone pairs of electrons, is known as n → π* interactions, and if it involves a σ* orbital in a linear fashion, then it is known as the carbon bond. These NCIs can be intra- or inter-molecular and are usually weak in strength but have a paramount effect on the structure and function of small-molecular crystals and proteins. Surprisingly, the experimental evidence of such interactions in the solution phase is scarce. It is even difficult to determine the interaction energy in the solution. Using NMR spectroscopy aided with molecular dynamics (MD) simulation and high-level quantum mechanical calculations, herein we provide the experimental evidence of intermolecular carbon-centered NCIs in solution. The challenge was to find appropriate heterodimers that could sustain room temperature thermal energy and collisions from the solvent molecules. However, after several trial model compounds, the pyridine-N-oxide:dimethyltetracyanocyclopropane (PNO-DMTCCP) complex was found to be a good candidate for the investigation. NBO analyses show that the PNO:DMTCCP complex is stabilized mainly by intermolecular n → π* interaction when a weaker carbon bond gives extra stability to the complex. From the NMR study, it is observed that the NCIs between DMTCCP and PNO are enthalpy driven with an enthalpy change of -28.12 kJ mol-1 and dimerization energy of ∼-38 kJ mol-1 is comparable to the binding energies of a conventional hydrogen-bonded dimer. This study opens up a new strategy to investigate weak intermolecular interactions such as n → π* interaction and carbon bonds in the solution phase.
| S-EPMC9749111 | biostudies-literature