Natural bond orbital analysis in the ONETEP code: Applications to large protein systems

First principles electronic structure calculations are typically performed in terms of molecular orbitals (or bands), providing a straightforward theoretical avenue for approximations of increasing sophistication, but do not usually provide any qualitative chemical information about the system. We can derive such information via post‐processing using natural bond orbital (NBO) analysis, which produces a chemical picture of bonding in terms of localized Lewis‐type bond and lone pair orbitals that we can use to understand molecular structure and interactions. We present NBO analysis of large‐scale calculations with the ONETEP linear‐scaling density functional theory package, which we have interfaced with the NBO 5 analysis program. In ONETEP calculations involving thousands of atoms, one is typically interested in particular regions of a nanosystem whilst accounting for long‐range electronic effects from the entire system. We show that by transforming the Non‐orthogonal Generalized Wannier Functions of ONETEP to natural atomic orbitals, NBO analysis can be performed within a localized region in such a way that ensures the results are identical to an analysis on the full system. We demonstrate the capabilities of this approach by performing illustrative studies of large proteins—namely, investigating changes in charge transfer between the heme group of myoglobin and its ligands with increasing system size and between a protein and its explicit solvent, estimating the contribution of electronic delocalization to the stabilization of hydrogen bonds in the binding pocket of a drug‐receptor complex, and observing, in situ, the n → π* hyperconjugative interactions between carbonyl groups that stabilize protein backbones. © 2012 Wiley Periodicals, Inc.     

Molecular docking,experimental FT-IR spectra,UV–Vis spectra,vibrational analysis,electronic properties,Fukui function analysis of a potential bioactive agent – Proflavine

Proflavine, having the molecular formula C₁₃H₁₁N₃, is a well-known urinary antiseptic and anticancer medication (3,6-diaminoacridine). In this communication, Quantum chemical computations of Proflavine's geometry have been performed and examined in the ground state. The optimized structure and wavenumbers of the molecule's vibrational bands were investigated using the DFT/B3LYP method and 6–311G (d, p) as the basis set. The calculated vibrational frequencies are compared to experimental IR spectra. The link between thermodynamic characteristics and temperature has been studied. The computed IR frequencies correlate well with the experiments, as indicated by the correlation factor (R² = 0.99). The UV spectra of the title molecule are calculated by using Time Dependent Density Functional Theory (TD-DFT). The molecule's interactions with other species were described using an analysis of a HOMO (Highest Occupied Molecular Orbital) and LUMO (Lowest Unoccupied Molecular Orbital). Natural bond orbitals (NBO) analysis was used to investigate intramolecular and intermolecular hydrogen bonding and their second-order stabilization energies and conjugative and hyperconjugative interactions. By computing the first hyperpolarizability, nonlinear optical (NLO) analysis was utilized to explore the molecule's nonlinear optical properties.     

One-pot synthesis,X-ray crystal structure,and identification of potential molecules against COVID-19 main protease through structure-guided modeling and simulation approach

Although antimicrobial resistance before the Covid-19 pandemic is a top priority for global public health, research is already ongoing on novel organic compounds with antimicrobial and antiviral properties in changing medical environments in connection with Covid 19. Thanks to the Biginelli reaction, which allows the synthesis of pyrimidine compounds, blockers of calcium channels, antibodies, antiviral, antimicrobial, anti-inflammatory, or antioxidant therapeutic compounds were investigated. In this paper, we aim to present Biginelli's synthesis, its therapeutic properties, and the structural–functional relationship in the test compounds that allows the synthesis of antimicrobial compounds. Both the DFT and TD-DFT computations of spectral data, molecular orbitals (HOMO, LUMO) analysis, and electrostatic potential (MEP) surfaces are carried out as an add-on to synthetic research. Hirshfeld surface analysis was also used to segregate the different intermolecular hydrogen bonds involved in the molecular packing strength. Natural Bond Orbital (NBO) investigation endorses the existence of intermolecular interactions mediated by lone pair, bonding, and anti-bonding orbitals. The dipole moment, linear polarizability, and first hyperpolarizabilities have been explored as molecular parameters. All findings based on DFT exhibit the best consistency with experimental findings, implying that synthesized molecules are highly stable. To better understand the binding mechanism of the SARS-CoV-2 M^pro, we performed molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations.     

Intra and intermolecular hydrogen bonding in formohydroxamic acid

The presence of hydrogen bonding interactions in several tautomeric forms of formohydroxamic acid (FHA) and 1:1 association among the tautomeric forms and water‐coordinated tautomeric forms of FHA is explored theoretically. Out of the seven equilibrium structures, four tautomeric forms have been selected for aggregation with single water molecule and dimer formation. Fifteen aggregates of FHA with H₂O have been optimized at MP2/AUG‐cc‐PVDZ level and analyzed for intramolecular and intermolecular H‐bond interactions. Twenty‐seven dimers of the four tautomeric forms have been obtained at MP2/6‐31+G* level. The stabilization energies associated with dimerization and adduct formation with water are the result of H‐bond interactions and range from very weak to medium. The atomic charges and NBO analysis indicate that the electrostatic and the charge transfer are the important components favoring H‐bond formation. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Theoretical Study and NBO Analysis of the Decomposition Kinetics of Oxetane, 2-Methyloxetane and 2,2-Dimethyloxetane

A theoretical study of the thermal decomposition kinetics of oxetane (1), 2-methyloxetane (2), and 2,2-dimethyloxetane (3) has been carried out at the B3LYP/6-311+G**, B3PW91/6-311+G**, and MPW1PW91/6-311+G** levels of theory. The MPW1PW91/6-311+G** method was found to give a reasonable good agreement with the experimental kinetics and thermodynamic parameters. The decomposition reaction of compounds 1～3 yields formaldehyde and the corresponding substituted olefin. Based on the optimized ground state geometries using MPW1PW91/6-311+G** method, the natural bond orbital (NBO) analysis of donor-acceptor (bond-antibond) interactions revealed that the stabilization energies associated with the electronic delocalization from σC3-C4 bonding to σ*O1-C2 antibonding orbitals decrease from compounds 1 to 3. The σC3-C4→σO1-C2 resonance energies for compounds 1～3 are 2.63, 2.59 and 2.45 kcal mol-1, respectively. Further, the results showed that the energy gaps between σC3-C4 bonding and σ*O1-C2 antibonding orbitals decrease from compounds 1 to 3. Also, the decomposition process in these compounds are controlled by σ→σ* resonance energies. Moreover, the obtained order of energy barriers could be explained by the number of electron-releasing methyl groups substituted to the Csp3 atom (which is attached to oxygen atom). NBO analysis shows that the occupancies of σCsp3-O bonds decrease for compounds 1～3 as 3<2<1, and those of σCsp3-O bonds increase in the opposite order (3 > 2 > 1). This fact illustrates a comparatively easier thermal decomposition of the sCsp3-O bond in compound 3 compared to compound 2, and in compound 2 compared to compound 1. NBO results indicate that these reactions are occurring through a concerted and asynchronous four-membered cyclic transition state type of mechanism.     

OCS,CO2,N2O与烯、炔烃之间π…π作用的电子密度拓扑研究

采用MP2/aug-cc-pVDZ方法对氧硫化碳(OCS)、二氧化碳(CO2)、一氧化二氮(N2O)与乙烯(C2H4)、乙炔(C2H2)、2-丁炔(C4H6)之间形成的平行构型复合物中的分子间相互作用进行了理论研究.复合物的相互作用能按照B…C2H4B…C2H2>B…C4H6(B=OCS,CO2,N2O)的顺序依次减小.采用电子密度拓扑分析理论方法,讨论了复合物中π…π作用的成键特性.电子密度拓扑分析表明复合物中形成了弱的分子间相互作用,且以静电作用为主;π电子密度分子图与全电子密度分子图中键径方向是一致的,说明π…π作用在本文所讨论的体系中起着很重要的作用.NBO分析表明净电荷迁移从电子给体C2H4,C2H2,C4H6到电子受体OCS,CO2,N2O,迁移数按照B…C2H4相似文献   

Halogen Interactions in Halogenated Oxindoles: Crystallographic and Computational Investigations of Intermolecular Interactions

X-ray structural determinations and computational studies were used to investigate halogen interactions in two halogenated oxindoles. Comparative analyses of the interaction energy and the interaction properties were carried out for Br···Br, C-H···Br, C-H···O and N-H···O interactions. Employing Møller–Plesset second-order perturbation theory (MP2) and density functional theory (DFT), the basis set superposition error (BSSE) corrected interaction energy (E_int(BSSE)) was determined using a supramolecular approach. The E_int(BSSE) results were compared with interaction energies obtained by Quantum Theory of Atoms in Molecules (QTAIM)-based methods. Reduced Density Gradient (RDG), QTAIM and Natural bond orbital (NBO) calculations provided insight into possible pathways for the intermolecular interactions examined. Comparative analysis employing the electron density at the bond critical points (BCP) and molecular electrostatic potential (MEP) showed that the interaction energies and the relative orientations of the monomers in the dimers may in part be understood in light of charge redistribution in these two compounds.     

Analysis of vibrational spectra of L-alanylglycine based on density functional theory calculations

FT Raman and IR spectra of the crystallized biologically active molecule, L-alanylglycine (L-Ala-Gly) have been recorded and analyzed. The equilibrium geometry, bonding features and harmonic vibrational frequencies of L-Ala-Gly have been investigated with the help of B3LYP density functional theory (DFT) method. The calculated molecular geometry has been compared with the experimental data. The assignments of the vibrational spectra have been carried out with the help of normal coordinate analysis (NCA) following the scaled quantum mechanical force field methodology (SQMFF). The optimized geometry shows the non-planarity of the peptide group of the molecule. Potential energy surface (PES) scan studies has also been carried out by ab initio calculations with B3LYP/6-311+G** basis set. The red shifting of NH3+ stretching wavenumber indicates the formation of N-H...O hydrogen bonding. The change in electron density (ED) in the sigma* antibonding orbitals and E2 energies have been calculated by natural bond orbital analysis (NBO) using DFT method. The NBO analysis confirms the occurrence of strong intermolecular hydrogen bonding in the molecule.     

Computation and analysis of 19F substituent chemical shifts of some bridgehead‐substituted polycyclic alkyl fluorides

The ¹⁹F NMR shieldings for several remotely substituted rigid polycyclic alkyl fluorides with common sets of substituents covering a wide range of electronic effects were calculated using the DFT‐GIAO theoretical model. The level of theory, B3LYP/6–311+G(2d,p), was chosen based on trial calculations which gave good agreement with experimental values where known. The optimized geometries were used to obtain various molecular parameters (fluorine natural charges, electron occupancies on fluorine of lone pairs and of the C? F bond, and hybridization states) by means of natural bond orbital (NBO) analysis which could help in understanding electronic transmission mechanisms underlying ¹⁹F substituent chemical shifts (SCS) in these systems. Linear regression analysis was employed to explore the relationship between the calculated ¹⁹F SCS and polar substituent constants and also the NBO derived molecular parameters. The ¹⁹F SCS are best described by an electronegativity parameter. The most pertinent molecular parameters appear to be the occupation number of the NBO p‐type fluorine lone pair and the occupation number of the C? F antibonding orbital. This trend suggests that in these types of rigid saturated systems hyperconjugative interactions play a key role in determining the ¹⁹F SCS. Electrostatic field effects appear to be relatively unimportant. Copyright © 2003 John Wiley & Sons, Ltd.     

Strongly Bound π-Hole Tetrel Bonded Complexes between H2SiO and Substituted Pyridines. Influence of Substituents

Ab initio calculation at the MP2/aug-cc-pVTZ level has been performed on the π-hole based N^…Si tetrel bonded complexes between substituted pyridines and H₂SiO. The primary aim of the study is to find out the effect of substitution on the strength and nature of this tetrel bond, and its similarity/difference with the N^…C tetrel bond. Correlation between the strength of the N^…Si bond and several molecular properties of the Lewis acid (H₂SiO) and base (pyridines) are explored. The properties of the tetrel bond are analyzed using AIM, NBO, and symmetry-adapted perturbation theory calculations. The complexes are characterized with short N^…Si intermolecular distances and high binding energies ranging between −142.72 and −115.37 kJ/mol. The high value of deformation energy indicates significant geometrical distortion of the monomer units. The AIM and NBO analysis reveal significant coordinate covalent bond character of the N⋅⋅⋅Si π-hole bond. Sharp differences are also noticed in the orbital interactions present in the N⋅⋅⋅Si and N⋅⋅⋅C tetrel bonds.     

Substituent effects on the electronic characteristics of pentacene derivatives for organic electronic devices: dioxolane-substituted pentacene derivatives with triisopropylsilylethynyl functional groups

The intramolecular electronic structures and intermolecular electronic interactions of 6,13-bis(triisopropylsilylethynyl)pentacene (TIPS pentacene), 6,14-bis-(triisopropylsilylethynyl)-1,3,9,11-tetraoxa-dicyclopenta[b,m]-pentacene (TP-5 pentacene), and 2,2,10,10-tetraethyl-6,14-bis-(triisopropylsilylethynyl)-1,3,9,11-tetraoxa-dicyclopenta[b,m]pentacene (EtTP-5 pentacene) have been investigated by the combination of gas-phase and solid-phase photoelectron spectroscopy measurements. Further insight has been provided by electrochemical measurements in solution, and the principles that emerge are supported by electronic structure calculations. The measurements show that the energies of electron transfer such as the reorganization energies, ionization energies, charge-injection barriers, polarization energies, and HOMO-LUMO energy gaps are strongly dependent on the particular functionalization of the pentacene core. The ionization energy trends as a function of the substitution observed for molecules in the gas phase are not reproduced in measurements of the molecules in the condensed phase due to polarization effects in the solid. The electronic behavior of these materials is impacted less by the direct substituent electronic effects on the individual molecules than by the indirect consequences of substituent effects on the intermolecular interactions. The ionization energies as a function of film thickness give information on the relative electrical conductivity of the films, and all three molecules show different material behavior. The stronger intermolecular interactions in TP-5 pentacene films lead to better charge transfer properties versus those in TIPS pentacene films, and EtTP-5 pentacene films have very weak intermolecular interactions and the poorest charge transfer properties of these molecules.     

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.     

Effect of hyperconjugation on ionization energies of hydroxyalkyl radicals

On the basis of electronic structure calculations and molecular orbital analysis, we offer a physical explanation of the observed large decrease (0.9 eV) in ionization energies (IE) in going from hydroxymethyl to hydroxyethyl radical. The effect is attributed to hyperconjugative interactions between the sigma CH orbitals of the methyl group in hydroxyethyl, the singly occupied p orbital of carbon, and the lone pair p orbital of oxygen. Analyses of vertical and adiabatic IEs and hyperconjugation energies computed by the natural bond orbital (NBO) procedure reveal that the decrease is due to the destabilization of the singly occupied molecular orbital in hydroxyethyl radical as well as structural relaxation of the cation maximizing the hyperconjugative interactions. The stabilization is achieved due to the contraction of the CO and CC bonds, whereas large changes in torsional angles bear little effect on the total hyperconjugation energies and, consequently, IEs.     

分子内结构环境对解离能的影响

为了理解化学键的这一结构效应, 本文对具有相同化学键而分子内结构环境不同的系列分子进行了计算研究, 讨论了化学键结构环境对解离能的影响.     

Intramolecular hydrogen bonding in structural conformers of 2‐amino methylene malonaldehyde: AIM and NBO studies

The molecular structure and intramolecular hydrogen bond energies of 44 conformers of 2‐Amino methylene malonaldehyde were investigated at MP2 and B3LYP levels of theory using the standard 6‐311++G** basis set and AIM and NBO analysis. The calculated geometrical parameters and conformational analysis in gas phase show that the closed ring via intramolecular hydrogen bonded conformers of this compound are more stable than the other ones. Hydrogen bond energies for H‐bonded conformers were obtained from the related rotamers method (RRM) and Schuster method, and also the nature of H‐bonding of them has been investigated by means of the Bader theory of atoms in molecules, which is based on topological properties of the electron density. Delocalization effects can be identified from the presence of off diagonal elements of the Fock matrix in the NBO basis. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Optical, redox, and NLO properties of tricyanovinyl oligothiophenes: comparisons between symmetric and asymmetric substitution patterns

A series of tricyanovinyl (TCV)-substituted oligothiophenes was synthesized and investigated with a number of physical methods including UV/Vis, IR, and Raman spectroscopy, nonlinear optical (NLO) measurements, X-ray diffraction, and cyclic voltammetry. Mono- or disubstituted oligomers were prepared by the reaction of tetracyanoethylene with mono- or dilithiated oligomers. The comparative effects of the symmetric and asymmetric substitutions in the electronic and molecular properties have been addressed. These oligomers display dramatic reductions in both their optical and electrochemical band gaps in comparison with unsubstituted molecules. The analysis of the electronic properties of the molecules was assisted by density functional theory calculations, which are in excellent agreement with the experimental data. TCV substitution influences the energies of the frontier orbitals, especially with respect to the stabilization of LUMO orbitals. X-ray structural characterization of a monosubstituted oligomer exhibits pi-stacking with favorable intermolecular interactions. NLO results agree with the role of the intramolecular charge-transfer feature in the asymmetric samples. These results furthermore exalt the role of conformational flexibility in the disubstituted compounds and reveal an unexpected nonlinear optical activity for symmetric molecules. Regarding the electronic structure, the interpretation of the vibrational data reflects the balanced interplay between aromatic and quinoid forms, finely tuned by the chain length and substitution pattern. The electronic and structural properties are consistent with the semiconducting properties exhibited by these materials in thin film transistors (TFTs).     

A highly distorted hexacoordinated silver(I) complex: synthesis,crystal structure,and DFT studies

A new highly distorted hexacoordinated silver(I) complex [AgL₂NO₃] with 2-(bis(methylthio)methylene)-1-phenylbutane-1,3-dione (L) as ligand is synthesized and characterized using elemental analysis, FTIR, NMR, and X-ray single-crystal structure analysis. The ligand (L) and the nitrate group act as bidentate ligands. The geometry around the silver ion has an intermediate configuration between a trigonal prism (TP) and an octahedron (OCT). Continuous shape measure (CShM) analysis indicated a closer configuration to TP than OCT. Experimentally and theoretically, the Ag–S bonds are shorter than any of the Ag–O bonds, indicating a stronger interaction between Ag⁺ (soft metal) and S-atom as a softer site than oxygen. Natural bond orbital (NBO) analyses showed higher interaction energies between the S-atom lone pairs and the Ag–antibonding NBO (8.61–31.39 kcal/mol) than LP(O)→Ag (3.48–11.46 kcal/mol). The acceptor antibonding NBO of the Ag atom has mainly s-orbital character. The Ag atom has a natural charge of +0.7579 e at the experimental structure, suggesting that negative charge was transferred from the ligand (0.0666 e) and nitrate (0.1090 e) to the Ag ion. Using Hirshfeld surface analysis, the important intermolecular interactions between molecular units within the crystal lattice of the ligand and its Ag-complex were analyzed and compared.     

Effect of green solvents physical,chemical, biological and bonding nature on 5-acetyl-thiophene-2-carboxylic acid by DFT and TD-DFT approach – An antiviral agent

The density functional theory (DFT) is applied to 5-acetyl-thiophene-2-carboxylic acid (5AT2CA). To determine the optimal structure and different physical, chemical, and biological characteristics, the B3LYP technique and 6–311++G(d,p) basis set are employed. The binding energies and ellipticity are determined using the atoms in molecules theory (AIM). NBO analysis is used to study the exchanges between the contributor and receiver by energies. By utilizing HOMO-LUMO values and important electronic parameters, stability is determined. The molecular electrostatic potential (MEP) and Fukui function from Mulliken charges are used to determine the reactive regions of the molecule. The TD-DFT technique obtains the electronic transition using UV–Vis spectrum with various solvents. NLO studies were carried out on this molecule. Studies of temperature effect on 5AT2CA are done by thermodynamic parameters. Drug-likeness and molecular docking tests are used to evaluate the bioactivity and antiviral properties.