scholarly journals Analysis of polarization in QM/MM modelling of biologically relevant hydrogen bonds

2008 ◽  
Vol 5 (suppl_3) ◽  
pp. 207-216 ◽  
Author(s):  
Kittusamy Senthilkumar ◽  
Jon I Mujika ◽  
Kara E Ranaghan ◽  
Frederick R Manby ◽  
Adrian J Mulholland ◽  
...  

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly important for the study of chemical reactions and systems in condensed phases. Here, we have tested the accuracy of a density functional theory-based QM/MM implementation (B3LYP/6-311+G(d,p)/CHARMM27) on a set of biologically relevant interactions by comparison with full QM calculations. Intermolecular charge transfer due to hydrogen bond formation is studied to assess the severity of spurious polarization of QM atoms by MM point charges close to the QM/MM boundary. The changes in total electron density and natural bond orbital atomic charges due to hydrogen bond formation in selected complexes obtained at the QM/MM level are compared with full QM results. It is found that charge leakage from the QM atoms to MM atomic point charges close to the QM/MM boundary is not a serious problem, at least with limited basis sets. The results are encouraging in showing that important properties of key biomolecular interactions can be treated well at the QM/MM level employing good-quality levels of QM theory.

2012 ◽  
Vol 2012 ◽  
pp. 1-12 ◽  
Author(s):  
Snehanshu Pal ◽  
T. K. Kundu

A detailed quantum chemical calculation based study of hydrogen bond formation in trimethylene glycol- (TMG-) water complex has been performed by Hatree-Fock (HF) method, second-order Møller-Plesset perturbation theory (MP2), density functional theory (DFT), and density functional theory with dispersion function (DFT-D) using 6-31++G(d,p) basis set. B3LYP DFT-D, WB97XD, M06, and M06-2X functionals are used to capture highly dispersive hydrogen bond formation. Geometrical parameters, interaction energy, deviation of potential energy curve of hydrogen-bonded O–H from that of free O–H, natural bond orbital (NBO), atom in molecule (AIM), charge transfer, and red shift are investigated. It is observed that hydrogen bond between TMG and water molecule is stronger in case of TMG acting as proton donor compared to that of water acting as proton donor, and dilute TMG solution would inhibit water cluster formation.


2010 ◽  
Vol 88 (8) ◽  
pp. 849-857 ◽  
Author(s):  
Nguyen Tien Trung ◽  
Tran Thanh Hue ◽  
Minh Tho Nguyen

The hydrogen-bonded interactions in the simple (HNZ)2 dimers, with Z = O and S, were investigated using quantum chemical calculations with the second-order Møller–Plesset perturbation (MP2), coupled-cluster with single, double (CCSD), and triple excitations (CCSD(T)) methods in conjunction with the 6-311++G(2d,2p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets. Six-membered cyclic structures were found to be stable complexes for the dimers (HNO)2, (HNS)2, and (HNO–HNS). The pair (HNS)2 has the largest complexation energy (–11 kJ/mol), and (HNO)2 the smallest one (–9 kJ/mol). A bond length contraction and a frequency blue shift of the N–H bond simultaneously occur upon hydrogen bond formation of the N–H···S type, which has rarely been observed before. The stronger the intramolecular hyperconjugation and the lower the polarization of the X–H bond involved as proton donor in the hydrogen bond, the more predominant is the formation of a blue-shifting hydrogen bond.


2021 ◽  
Author(s):  
Qian Tang ◽  
Ting Huang ◽  
Ruisi Huang ◽  
Hongyu Cao ◽  
Lihao Wang ◽  
...  

Abstract The hydrogen bond formation with formic acid would affect the complementary pair of bases between uracil and adenine, but the binding modes and spectral properties of hydrogen bonds are still obscure. Density functional theory and time-dependent density functional theory were applied to investigate the intermolecular hydrogen bonds between uracil and formic acid. The reduced density gradient (RDG), bond lengths and vibration absorption frequencies revealed that the most probable uracil-formic acid (U-FA) interaction mode formed in the position c of FA and the site 1 of U, that is, the mode 1c. The theoretical parameters in excited state complexes manifested that the variety of hydrogen bond configurations led to different degrees of strengthening or weakening of molecular interaction. In the implicit solvent (water), the formations of O-H∙∙∙O in the uracil-formic acid complexes were promoted obviously. These theoretical studies would positively affect the researches of life science and medicinal chemistry.


2015 ◽  
Vol 51 (1) ◽  
pp. 213-219 ◽  
Author(s):  
Farhoush Kiani ◽  
Mehran Abbaszadeh ◽  
Mohammad Pousti ◽  
Fardad Koohyar

In the present work, acid dissociation constant (pKa) values of muscimol derivatives were calculated using the Density Functional Theory (DFT) method. In this regard, free energy values of neutral, protonated and deprotonated species of muscimol were calculated in water at the B3LYP/6-31G(d) basis sets. The hydrogen bond formation of all species had been analyzed using the Tomasi's method. It was revealed that the theoretically calculated pKa values were in a good agreement with the existing experimental pKa values, which were determined from capillary electrophoresis, potentiometric titration and UV-visible spectrophotometric measurements.


2013 ◽  
Vol 2013 ◽  
pp. 1-16 ◽  
Author(s):  
Snehanshu Pal ◽  
T. K. Kundu

A detailed theoretical study of hydrogen-bond formation in different polyethylene glycol + water complex and dipropylene glycol + water have been performed by Hartree Fock (HF) method, second-order Møller-Plesset perturbation theory (MP2), and density functional theory (DFT) using 6-31++G(d,p) basis set. B3LYP DFT-D, WB97XD, M06, and M06-2X functionals have been used to describe highly dispersive hydrogen-bond formation appropriately. Geometrical parameters, interaction energies, deformation energies, deviation of potential energy curves of hydrogen bonded O–H from that of free O–H, frontier orbitals, and charge transfer have been studied to analyze stability and nature of hydrogen bond formation of various glycol and water complexes. It is found that WB97XD is best among all the applied DFT functionals to describe hydrogen bond interaction, and intermolecular hydrogen bonds have higher covalent character and accordingly higher strength when glycol acts as proton donor for glycol + 1 water complex system.


2021 ◽  
pp. 3323-3335
Author(s):  
Mustafa M. Kadhim ◽  
Layla A. Al. Juber ◽  
Ahmed S. M. Al-Janabi

    Seven Zn-dithiocarbamate complexes were suggested as corrosion inhibitors. Density functional theory (DFT) was used to predict the ability of inhibition. Room temperature conditions were applied to suggest the optimization of complexes, physical properties, and corrosion parameters. In addition, the HOMO, LUMO, dipole moment, energy gap, and other parameters were used to compare the inhibitors efficiency. Gaussian 09 software with LanL2DZ basis set was used. Total electron density (TED) and electrostatic surface potential (ESP) were utilized to show the sites of adsorption according to electron density.


2019 ◽  
Author(s):  
Kamal Batra ◽  
Stefan Zahn ◽  
Thomas Heine

<p>We thoroughly benchmark time-dependent density- functional theory for the predictive calculation of UV/Vis spectra of porphyrin derivatives. With the aim to provide an approach that is computationally feasible for large-scale applications such as biological systems or molecular framework materials, albeit performing with high accuracy for the Q-bands, we compare the results given by various computational protocols, including basis sets, density-functionals (including gradient corrected local functionals, hybrids, double hybrids and range-separated functionals), and various variants of time-dependent density-functional theory, including the simplified Tamm-Dancoff approximation. An excellent choice for these calculations is the range-separated functional CAM-B3LYP in combination with the simplified Tamm-Dancoff approximation and a basis set of double-ζ quality def2-SVP (mean absolute error [MAE] of ~0.05 eV). This is not surpassed by more expensive approaches, not even by double hybrid functionals, and solely systematic excitation energy scaling slightly improves the results (MAE ~0.04 eV). </p>


2021 ◽  
Vol 03 (02) ◽  
pp. 090-096
Author(s):  
Yusuke Ishigaki ◽  
Kota Asai ◽  
Takuya Shimajiri ◽  
Tomoyuki Akutagawa ◽  
Takanori Fukushima ◽  
...  

The crystal structures of a series of tetracyanonaphthoquinodimethanes fused with a selenadiazole or thiadiazole ring revealed that their molecular packing is determined mainly by two intermolecular interactions: chalcogen bond (ChB) and weak hydrogen bond (WHB). ChB between Se and a cyano group dictates the packing of selenadiazole derivatives, whereas the S-based ChB is much weaker and competes with WHB in thiadiazole analogues. This difference can be explained by different electrostatic potentials as revealed by density functional theory calculations. A proper molecular design that weakens WHB can change the contribution of ChB in determining the crystal packing of thiadiazole derivatives.


Sign in / Sign up

Export Citation Format

Share Document