An ab initio MO study on the structures and electronic states of hydrogen-bonded O3- HF and SO2-HF complexes

Ab initio molecular orbital (MO) calculations have been carried out for base-hydrogen fluoride (HF) complexes (base = O3 and SO2) in order to elucidate the structures and energetics of the complexes. The ab initio calculations were performed up to the QCISD(T)/6-311++G(d,p) level of theory. In both complexes, hydrogen-bonded structures where the hydrogen of HF orients toward one of the oxygen atoms of bases were obtained as stable forms. The calculations showed that cis and trans isomers exist in both complexes. All calculations for the SO2-HF complex indicated that the cis form is more stable in energy than the trans form. On the other hand, in O3-HF complexes, the stable structures are changed by the ab initio levels of theory used, and the energies of the cis and trans forms are close to each other. From the most sophisticated calculations (QCISD(T)/6-311++G(d,p)//QCISD/6-311+G(d) level), it was predicted that the complex formation energies for cis SO2-HF, trans SO2-HF, cis O3-HF, and trans O3-HF are 6.1, 5.7, 3.4, and 3.6 kcal/mol, respectively, indicating that the binding energy of HF to SO2 is larger than that of O3. The harmonic vibrational frequencies calculated for cis O3-HF and cis SO2-HF complexes were in good agreement with the experimental values measured by Andrews et al. Also, the calculated rotation constants for cis SO2-HF agreed with the experiment.     

Quantum mechanical map for protein-ligand binding with application to beta-trypsin/benzamidine complex

We report full ab initio Hartree-Fock calculation to compute quantum mechanical interaction energies for beta-trypsin/benzamidine binding complex. In this study, the full quantum mechanical ab initio energy calculation for the entire protein complex with 3238 atoms is made possible by using a recently developed MFCC (molecular fractionation with conjugate caps) approach in which the protein molecule is decomposed into amino acid-based fragments that are properly capped. The present MFCC ab initio calculation enables us to obtain an "interaction spectrum" that provides detailed quantitative information on protein-ligand binding at the amino acid levels. These detailed information on individual residue-ligand interaction gives a quantitative molecular insight into our understanding of protein-ligand binding and provides a guidance to rational design of potential inhibitors of protein targets.     

The structure and binding energies of the van der Waals complexes of Ar and N2 with phenol and its cation, studied by high level ab initio and density functional theory calculations

We have investigated, using both ab initio and density functional theory methods, the minimum energy structures and corresponding binding energies of the van der Waals complexes between phenol and argon or the nitrogen molecule, and the corresponding complexes involving the phenol cation. Structures were obtained at the MP2 level using a large basis, and the corresponding energies were corrected for basis set superposition error (BSSE), higher order electron correlation effects, and for basis set size. The structures of the global minima were further refined for the effects of BSSE and the corresponding binding energies were evaluated. For each neutral species, we find only a single true minimum, pi bonded for argon and OH bonded for nitrogen. For both cationic species, we find that the OH-bonded complex is preferred over other minima which we have identified as having Ar or N(2) between exogeneous atoms. The ab initio calculations are generally in excellent agreement with experimental binding energies and rotational constants. We find that the B3LYP functional is particularly poor at describing these complexes, while a density functional theory (DFT) method with an empirical correction for dispersive interactions (DFT-D) is very successful, as are some of the new functionals proposed by Zhao and Truhlar [J. Phys. Chem. A 109, 5656 (2005); J. Chem. Theory Comput. 2, 1009 (2006); Phys. Chem. Chem. Phys. 7, 2701 (2005); J. Phys. Chem. A 108, 6908 (2004)]. Both the ab initio and DFT-D methods accurately predict the intermolecular vibrational modes.     

Molecular interactions between estrogen receptor and its ligand studied by the ab initio fragment molecular orbital method

The ab initio fragment molecular orbital calculations were performed for molecular interactions of the whole estrogen receptor (ER) ligand-binding domain with a natural ligand, 17beta-estradiol (EST). The interaction energies of the ligand at the residue level were calculated using HF and MP2 methods with several basis sets. The charge-transfer (CT) interactions were also analyzed based on configuration analysis for fragment interaction. Strong electrostatic interactions were observed between the EST and surrounding charged/polarized residues, Glu353, Arg394, His524, and Thr347. Weak electrostatic and significant van der Waals dispersion interactions were observed between the EST and the many surrounding hydrophobic residues. Together with the experimental interpretations, both interactions equally contributed to the total binding energies, and it was found that the inclusion of electron correlation was essential to obtain an appropriate picture of the interaction. The strongest interaction energy was observed between Glu353 and the EST, and the CT interactions from the lone-pair orbital of the carbonyl oxygen of Glu353 to the sigma(OmicronEta) orbital of the hydroxyl group of EST were found to be important. The CT interactions from the lone-pair orbital of EST to the sigma(NuEta) of Arg394 and from the lone-pair orbital of EST to the sigma(NuEta) of His524 were also observed. These CT interactions occurred through the hydrogen-bond networks between the ER and EST. Therefore, electron donations from the ER to the EST and electron back-donations from EST to the ER were characteristic of ER-ligand binding. Our approach provides a powerful tool to understanding detailed molecular interactions at the quantum mechanical level.     

Force field parametrization and molecular dynamics simulation of flexible POSS-linked (NHC; phosphine) Ru catalytic complexes

In recent years, N-heterocyclic carbene (NHC) or phospine groups have been put forward as candidate catalysts ligands for olefin metathesis reactions to be performed using multistep methods. Some of these proposed ligands contain polyhedral oligomeric silsesquioxane (POSS) structures linked to NHC rings by means of alkyl chains. Some important properties for the prediction of catalytic activity, such as the theoretically defined buried volume, are related to the conformational characteristics of these complex ligands that can be studied through molecular dynamics simulations. However, the chemical structure of resulting catalytic complexes usually contains atoms or groups that are not included in the common forcefields used in simulations. In this work we focus on complexes formed by a catalytic metal center (Ru) with both phospine and POSS-linked NHC groups. The central part of the complexes contain atoms and groups that have bonds, bond angles, and torsional angles whose parameters have not been previously evaluated and included in existing force fields. We have performed basic ab initio quantum mechanical calculations based on the density functional theory to obtain energies for this central section. The force field parameters for bonds, bond angles, and torsional angles are then calculated from an analysis of energies calculated for the equilibrium and different locally deformed structures. Nonbonded interactions are also conveniently evaluated. From subsequent molecular dynamics simulations, we have obtained results that illustrate the conformational characteristics most closely connected with the catalytic activity.     

Semiempirical PM5 molecular orbital study on chlorophylls and bacteriochlorophylls: comparison of semiempirical, ab initio, and density functional results

The semiempirical PM5 method has been used to calculate fully optimized structures of magnesium-bacteriochlorin, magnesium-chlorin, magnesium-porphin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), and d, and bacteriochlorophylls a, b, c, d, e, f, g, and h with all homologous structures. Hartree-Fock/6-31G* ab initio and density functional B3LYP/6-31G* methods were used to optimize structures of methyl chlorophyllide a, chlorophyll c(1), and methyl bacteriochlorophyllides a and c for comparison. Spectroscopic transition energies of the chromophores and their 1:1 or 1:2 solvent complexes were calculated with the Zindo/S CIS method. The self-consistent reaction field model was used to estimate solvent shifts. The PM5 calculations predict planar structure of the porphyrin ring and central position of the four coordinated magnesium atoms in all pigments studied, in accord with the experimental, ab initio, and density functional results, a significant improvement as compared to the older semiempirical PM3 approach. Only small differences in PM5 and B3LYP/6-31G* or Hartree-Fock/6-31G* minimum energy geometries of the reference molecules were observed. Calculations show that in 1:1 solvent complexes, where the magnesium atom is five coordinated, the magnesium atom is shifted out of the plane of the porphyrin ring towards the solvent molecule, while the hexa coordinated 1:2 complexes are again planar. The PM5 method gives atomic charges that are comparable with those obtained from the Hartree-Fock/6-31G* and B3LYP/6-31G* calculations. The single point ZINDO/S CIS calculations with PM5 minimum energy structure gave excellent correlations between calculated and experimental transition energies of the chlorophylls and bacteriochlorophylls studied. Such correlations may be used for prediction of transition energies of the chromophores in protein binding sites. Calculations also predict existence of dark electronic states below the main Soret absorption band in all chromophores studied. The results suggest that the semiempirical PM5 method is a fairly reliable and computationally efficient method in predicting molecular parameters of porphyrin-like molecules.     

Quantum chemical simulation of excited states of chlorophylls, bacteriochlorophylls and their complexes

The present review describes the use of quantum chemical methods in estimation of structures and electronic transition energies of photosynthetic pigments in vacuum, in solution and imbedded in proteins. Monomeric Mg-porphyrins, chlorophylls and bacteriochlorophylls and their solvent 1:1 and 1:2 complexes were studied. Calculations were performed for Mg-porphyrin, Mg-chlorin, Mg-bacteriochlorin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), d and bacteriochlorophylls a, b, c, d, e, f, g, h, plus several homologues. Geometries were optimised with PM3, PM3/CISD, PM5, ab initio HF (6-31G*/6-311G**) and density functional B3LYP (6-31G*/6-311G**) methods. Spectroscopic transition energies were calculated with ZINDO/S CIS, PM3 CIS, PM3 CISD, ab initio CIS, time-dependent HF and time-dependent B3LYP methods. Estimates for experimental transition energies were obtained from linear correlations of the calculated transition energies of 1:1 solvent complexes against experimentally recorded solution energies (scaling). According to the calculations in five-coordinated solvent complexes the magnesium atom lies out of the porphyrin plane, while in six-coordinated complexes the porphyrin is nearly planar. Charge densities on magnesium and nitrogen atoms were strongly dependent on the computational method deployed. Several dark states of low oscillator strength below the main Soret band were predicted for solvent complexes and chlorophylls and bacteriochlorophylls in protein environment. Such states, though not yet identified experimentally, might serve as intermediate states for excitation energy transfer in photosynthetic complexes. Q(y), Q(x) and Soret transition energies were found to depend on the orientation of the acetyl group and external pressure. A method to estimate site energies and dimeric interaction energies and to simulate absorption and CD spectra of photosynthetic complexes is described. Simulations for the light harvesting complexes Rhodospirillum molischianum, chlorosomes of Chlorobium tepidum and Chloroflexus aurantiacus, and LHC-II of Spinacia oleracea are presented as examples.     

Spectroscopy and structures of copper complexes with ethylenediamine and methyl-substituted derivatives

Copper complexes of ethylenediamine (en), N-methylethylenediamine (meen), N,N-dimethylethylenediamine (dmen), N,N,N'-trimethylethylenediamine (tren), and N,N,N',N'-tetramethylethylenediamine (tmen) are synthesized in laser-vaporization supersonic molecular beams and studied by pulsed-field ionization zero electron kinetic energy (ZEKE) and photoionization efficiency spectroscopies and second-order Moller-Plesset perturbation theory. Precise ionization energies and vibrational frequencies of Cu-en, -meen, and -dmen are measured from the ZEKE spectra, and ionization thresholds of Cu-tren and -tmen are estimated from the photoionization efficiency spectra. The measured vibrational modes span a frequency range of 35-1646 cm(-1) and include metal-ligand stretch and bend, hydrogen-bond stretch, and ligand-based torsion. A number of low-energy structures with Cu binding to one or two nitrogen atoms are predicted for each complex by the ab initio calculations. The combination of the spectroscopic measurements and ab initio calculations has identified a hydrogen-bond-stabilized monodentate structure for the Cu-en complex and bidentate cyclic structures for the methyl-substituted derivatives. The change of the Cu binding from the monodentate to the bidentate mode arises from the competition between copper coordination and hydrogen bonding.     

Approximate additivity of anion-pi interactions: an ab initio study on anion-pi, anion-pi(2) and anion-pi(3) complexes

We have studied the additivity of the anion-pi interaction using high level ab initio calculations. We have optimized chloride and bromide complexes with one, two and three aromatic units (such as trifluoro-s-triazine and s-triazine). We have analyzed the interaction using the atoms in molecules theory and studied the charge transfer using several methods for deriving atomic charges. The results revealed additivities of both the geometries and the binding energies. We have also proposed a neutral receptor for chloride based on multiple anion-pi interactions. Finally, we have simulated solvent effects within the self-consistent reaction field model.     

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.     

Theoretical analysis of binding specificity of influenza viral hemagglutinin to avian and human receptors based on the fragment molecular orbital method

The hemagglutinin (HA) protein of the influenza virus binds to the host cell receptor in the early stage of viral infection. A change in binding specificity from avian 2-3 to human 2-6 receptor is essential for optimal human-to-human transmission and pandemics. Therefore, it is important to reveal the key factors governing the binding affinity of HA-receptor complex at the molecular level for the understanding and prediction of influenza pandemics. In this work, on the basis of ab initio fragment molecular orbital (FMO) method, we have carried out the interaction energy analysis of HA-receptor complexes to quantitatively elucidate the binding specificity of HAs to avian and human receptors. To discuss the binding property of influenza HA comprehensively, a number of HAs from human H1, swine H1, avian H3 and avian H5 viruses were analyzed. We performed detailed investigations about the interaction patterns of complexes of various HAs and receptor analogues, and revealed that intra-molecular interactions between conserved residues in HA play an important role for HA-receptor binding. These results may provide a hint to understand the role of conserved acidic residues at the receptor binding site which are destabilized by the electrostatic repulsion with sialic acid. The calculated binding energies and interaction patterns between receptor and HAs are consistent with the binding specificities of each HA and thus explain the receptor binding mechanism. The calculated results in the present analysis have provided a number of viewpoints regarding the models for the HA-receptor binding specificity associated with mutated residues. Examples include the role of Glu190 and Gln226 for the binding specificity of H5 HA. Since H5 HA has not yet been adapted to human receptor and the mechanism of the specificity change is unknown, this result is helpful for the prediction of the change in receptor specificity associated with forthcoming possible pandemics.     

Can MM/GBSA calculations be sped up by system truncation?

We have studied whether calculations of the binding free energy of small ligands to a protein by the MM/GBSA approach (molecular mechanics combined with generalized Born and surface area solvation) can be sped up by including only a restricted number of atoms close to the ligand. If the protein is truncated before the molecular dynamics (MD) simulations, quite large changes are observed for the calculated binding energies, for example, 4 kJ/mol average difference for a radius of 19 Å for the binding of nine phenol derivatives to ferritin. The results are improved if no atoms are fixed in the simulations, with average and maximum errors of 2 and 3 kJ/mol at 19 Å and 3 and 6 kJ/mol at 7 Å. Similar results are obtained for two additional proteins, p38α MAP kinase and factor Xa. On the other hand, if energies are calculated on snapshots that are truncated after the MD simulation, all residues more than 8.5 Å from the ligand can be omitted without changing the energies by more than 1 kJ/mol on average (maximum error 1.4 kJ/mol). At the molecular mechanics level, the gain in computer time for such an approach is small. However, it shows what size of system should be used if the energies instead are calculated with a more demanding method, for example, quantum‐mechanics. © 2017 Wiley Periodicals, Inc.     

Comparative study of imidazole hydration: Ab initio and electrostatic calculations vs. Cambridge structural database analysis

We studied geometries and energies of complexes between water and neutral or protonated imidazole by ab initio molecular orbital calculations using the 4-31G basis set with and without the counterpoise correction. Positions of hydration sites and relative binding energies could be also estimated by using the electrostatic field map of imidazole as calculated by our bond increment method. The reliability of the calculations is confirmed by comparing the geometries of the imidazole-water complex to the experimental ones from the Cambridge Structural Database. These were obtained by X-ray diffraction studies on crystals with water bound to a molecule containing the imidazole fragment.     

NDDO/MC: A new semiempirical SCFMO method for transition metal complexes

A new semiempirical SCF MO procedure available for prediction of the transition-metal complexes' binding energy and molecular geometry is developed. The features of the method are (i) an explicit account of the orthogonality of the basis set; (ii) use of a new formula for the resonance integral; and (iii) an effective account of the Coulomb correlation of electrons in the calculation of the two-electron integrals based on approach of a model Coulomb hole function. The parameterization for H, C, N, O, Co, and Ni atoms is presented. The results of NDDO /MC (NDDO for Metal Compounds) calculations of molecular geometries and binding energies for a number of organic compounds and more than 30 cobalt and nickel complexes are compared with the available experimental and ab initio data. The average absolute errors for the binding energies of organic molecules and metal complexes are 8.8 and 5.0 kcal/mol, respectively. © 1992 John Wiley & Sons, Inc.     

Representation of Zn(II) complexes in polarizable molecular mechanics. Further refinements of the electrostatic and short-range contributions. Comparisons with parallel ab initio computations

We present refinements of the SIBFA molecular mechanics procedure to represent the intermolecular interaction energies of Zn(II). The two first-order contributions, electrostatic (E(MTP)), and short-range repulsion (E(rep)), are refined following the recent developments due to Piquemal et al. (Piquemal et al. J Phys Chem A 2003, 107, 9800; and Piquemal et al., submitted). Thus, E(MTP) is augmented with a penetration component, E(pen), which accounts for the effects of reduction in electronic density of a given molecular fragment sensed by another interacting fragment upon mutual overlap. E(pen) is fit in a limited number of selected Zn(II)-mono-ligated complexes so that the sum of E(MTP) and E(pen) reproduces the Coulomb contribution E(c) from an ab initio Hartree-Fock energy decomposition procedure. Denoting by S, the overlap matrix between localized orbitals on the interacting monomers, and by R, the distance between their centroids, E(rep) is expressed by a S(2)/R term now augmented with an S(2)/R(2) one. It is calibrated in selected monoligated Zn(II) complexes to fit the corresponding exchange repulsion E(exch) from ab initio energy decomposition, and no longer as previously the difference between (E(c) + E(exch)) and E(MTP). Along with the reformulation of the first-order contributions, a limited recalibration of the second-order contributions was carried out. As in our original formulation (Gresh, J Comput Chem 1995, 16, 856), the Zn(II) parameters for each energy contribution were calibrated to reproduce the radial behavior of its ab initio HF counterpart in monoligated complexes with N, O, and S ligands. The SIBFA procedure was subsequently validated by comparisons with parallel ab initio computations on several Zn(II) polyligated complexes, including binuclear Zn(II) complexes as in models for the Gal4 and beta-lactamase metalloproteins. The largest relative error with respect to the RVS computations is 3%, and the ordering in relative energies of competing structures reproduced even though the absolute numerical values of the ab initio interaction energies can be as large as 1220 kcal/mol. A term-to-term identification of the SIBFA contributions to their ab initio counterparts remained possible even for the largest sized complexes.     

The binding energies of NO-Rg (Rg = He, Ne, Ar) determined by velocity map imaging

We report velocity map imaging measurements of the binding energies, D(0), of NO-Rg (Rg = He, Ne, Ar) complexes. The X state binding energies determined are 3.0 ± 1.8, 28.6 ± 1.7, and 93.5 ± 0.9 cm(-1) for NO-He, -Ne, and -Ar, respectively. These values compare reasonably well with ab initio calculations. Because the ?-X transitions were unable to be observed for NO-He and NO-Ne, values for the binding energies in the ? state of these complexes have not been determined. Based on our X state value and the reported ?-X origin band position, the ? state binding energy for NO-Ar was determined to be 50.6 ± 0.9 cm(-1).     

Vibrational and rotational cooling of NO+ in collisions with He

A quantum mechanical investigation of the vibrational and rotational deactivation of NO(+) in collisions with He atoms in the cold and ultracold regime is presented. Ab initio potential energy calculations are carried out at BCCD(T) level and a new global 3D potential energy surface (PES) is obtained by fitting ab initio points within the reproducing kernel Hilbert space method. As a first test of this PES the bound state energies of the (3)He-NO(+) and (4)He-NO(+) complexes are calculated and compared to previous rigid rotor calculations. The efficiency of the vibrational and the rotational cooling of this molecular ion using a buffer gas of helium is then investigated by performing close coupling scattering calculations for collision energy ranging from 10(-6) to 2000 cm(-1). The calculations are performed for the two isotopes (3)He and (4)He and the results are compared to the available experimental data.     

Classical trajectory study of collisions of Ar with alkanethiolate self-assembled monolayers: potential-energy surface effects on dynamics

We have investigated collisions between Ar and alkanethiolate self-assembled monolayers (SAMs) using classical trajectory calculations with several potential-energy surfaces. The legitimacy of the potential-energy surfaces is established through comparison with molecular-beam data and ab initio calculations. Potential-energy surfaces used in previous work overestimate the binding of Ar to the SAM, leading to larger energy transfer than found in the experiments. New calculations, based on empirical force fields that better reproduce ab initio calculations, exhibit improved agreement with the experiments. In particular, polar-angle-dependent average energies calculated with explicit-atom potential-energy surfaces are in excellent agreement with the experiments. Polar- and azimuthal-angle-dependent product translational energies are examined to gain deeper insight into the dynamics of Ar+SAM collisions.     

Sexual attraction in the silkworm moth. Nature of binding of bombykol in pheromone binding protein--an ab initio study

An analysis of the crystal structure of [BmPBP...bombykol] complex identified nine amino acid residues involved in a variety of intermolecular interactions binding the ligand. Using simple model fragments as the representatives of the residues, the interaction energies of their complexes with bombykol were calculated using high-level ab initio methods. The results were discussed in terms of the method and basis set dependence and were further corrected to account for their pair nonadditivities. This enabled us to describe quantitatively the nature and origin of the binding forces in terms of contribution of the individual amino acids and individual types of interaction to the overall stability. All of these interactions are well defined and cannot be considered as nonspecific hydrophobic interactions, one of the major conclusions of this work.     

Ab initio study of phosphorescent emitters based on rare-earth complexes with organic ligands for organic electroluminescent devices

An ab initio approach is developed for calculation of low-lying excited states in Ln(3+) complexes with organic ligands. The energies of the ground and excited states are calculated using the XMCQDPT2/CASSCF approximation; the 4f electrons of the Ln(3+) ion are included in the core, and the effects of the core electrons are described by scalar quasirelativistic 4f-in-core pseudopotentials. The geometries of the complexes in the ground and triplet excited states are fully optimized at the CASSCF level, and the resulting excited states have been found to be localized on one of the ligands. The efficiency of ligand-to-lanthanide energy transfer is assessed based on the relative energies of the triplet excited states localized on the organic ligands with respect to the receiving and emitting levels of the Ln(3+) ion. It is shown that ligand relaxation in the excited state should be properly taken into account in order to adequately describe energy transfer in the complexes. It is demonstrated that the efficiency of antenna ligands for lanthanide complexes used as phosphorescent emitters in organic light-emitting devices can be reasonably predicted using the procedure suggested in this work. Hence, the best antenna ligands can be selected in silico based on theoretical calculations of ligand-localized excited energy levels.