Optimization of parameters for semiempirical methods IV: extension of MNDO,AM1, and PM3 to more main group elements

The NDDO semiempirical methods MNDO, AM1, and PM3 have been extended to all the remaining non-radioactive elements of the main group, excluding the noble gases. Most of the new elements are of Groups I and II. 44 sets of parameters are presented for the following methods and elements. MNDO: Na, Mg, K, Ca, Ga, As, Se, Rb, Sr, In, Sb, Te, Cs, Ba, Tl, and Bi; AM1: Li, Be, Na, Mg, K, Ca, Ga, As, Se, Rb, Sr, In, Sn, Sb, Te, Cs, Ba, Tl, Pb, and Bi; PM3: B, Na, K, Ca, Rb, Sr, Cs, and Ba. Average errors are presented for heats of formation, molecular geometries, etc.     

A semiempirical study of the prototropic tautomerism of hypoxanthine

The total and relative energies, bond order matrices and localized MOs for the eight possible tautomers of hypoxanthine (HYP) have been calculated, with full geometry optimization, using both AM1 and MNDO methods. The AM1 relative energies show that HYP(9,1), HYP(7,1) and HYP (9,10) are the predominant species at room temperature, the two former being in larger concentration that the latter. The calculated IR spectra for these species agree well with the reported spectrum in an isolated matrix, which has been interpreted in terms of the presence of these three tautomeric forms. The MNDO method does not predict the right order, and the more stable tautomer would be HYP(9,10). The calculated structure for the HYP(9,1) species shows that the molecule is essentially planar. The bond distances compare well with those of hypoxanthine hydrochloride and guanine and also correlate well with the calculated bond orders. The proton affinities for the three more stable tautomers have also been calculated. For HYP(9,1) the prefered site of protonation is N7, whereas for HYP(7,1) the protonation occurs rather at N9. These results agree well with¹⁵N and¹³C NMR studies in DMSO.     

New Multicentre Point Charge Models for Molecular Electrostatic Potentials from Semiempirical M0-Calculations

Two quasi-multipole electrostatic models for molecular charge distributions are presented. They assign arrays of point charges to nonhydrogen atoms on the basis of hybrid orbitals or localised molecular orbitals. When used with common semiempirical MO-techniques, they reproduce natural atomic orbital derived point charge (NAO-PC) and ab initio molecular potentials well. The localised orbital technique (LMO-PC) is intuitively more attractive than the hybrid orbital-point charge (HO-PC) method, although the former is more CPU-intensive.Electronic Supplementary Material available.     

On the performance of semiempirical MO theory for dipole moments of dye molecules

The dipole moments of a set of 71 simple dye molecules calculated at the ab initio, DFT, and semiempirical levels have been compared. The DFT dipole moments are on average 16% larger than those obtained by MP2/6-31G**. AM1 and PM3 modified with an empirical correction procedure yield dipole moments essentially at the same level of accuracy as the results of non-empirical calculations. INDO/S and CNDO/S are considerably less accurate. Among different versions of spectral methods, the CISD scheme gives the best performance.Electronic Supplementary Material available.     

A molecular mechanics and semiempirical molecular orbital study on the conformation of polynorbornene chains

The conformational analysis of polynorbornene (PNB) chains was investigated with the AM1, MM2, AMBER and OPLS methods taking into consideration the possibility of binding of norbornene monomers to each other at various positions, i.e. exo–exo, exo–endo, endo–endo. The chain that is formed by connecting exo–endo positions of the monomers has lower torsional barrier energy than those formed with bonds at other positions and has more flexibility. It is determined that the thredisyndiotactic chain formed by exo–endo addition adopts a helix structure and has a coil shape. The disyndiotactic chain formed by connecting norbornene monomers in mixed type has a linear structure. It is found that the repeat unit conformations of thredisyndiotactic and disyndiotactic chains of PNB are TGTG^– and (TGTG^–)₂, respectively.     

Computational Investigations of the Chalcogen-Substituted Carboxylic Acids RC(=O)XH and RC(=X)OH and their Dimers [RC(=O)XH]2 and [RC(=X)OH]2 (X=S, Se, Te)

Semiempirical and ab initio theoretical methods have been used to investigate molecular structures of the chalcogen-substituted carboxylic acid isomers RC(=O)XH (chalcogenol acid) and RC(=X)OH (chalcogenon acid). A recent experimental report suggests that the chalcogenon isomers, although less stable at room temperature, predominate at low temperature in polar solvents and that there is only a small barrier to isomerization between the isomers. Theoretical calculations have been used to locate minimum energy structures of chalcogen-substituted carboxylic acid isomers and to calculate energy differences between pairs of isomers. Carboxylic acids are well known to dimerize, especially in the gas phase and in non-polar solvents. We have, therefore, also calculated energies of dimerization of the chalcogen-substituted acids by optimizing the geometries of the symmetric dimers. We note that the PM3 level of theory is only qualitatively correct for sulfur- and selenium-containing species but fails even qualitatively for the tellurium-containing compounds. Ab initio results confirm the experimental observations and provide good estimates of both isomerization and dimerization energies. We conclude that for many functional groups with tautomers RC(=X)YH and RC(=Y)XH, the more acidic tautomer is the one with the acid proton on the smaller, more electronegative atom, although in many cases this may not be the more stable tautomer.Electronic Supplementary Material available.     

DFT study of the gas phase proton transfer in guanine assisted by water, methanol, and hydrogen peroxide

A computational study of hydrogen-bonded complexes between the oxo-/hydroxo-amino N7/9H tautomers of guanine and water, methanol, and hydrogen peroxide has been performed at the B3LYP/6-31+G(d) level of theory. The mechanisms of the water-, methanol-, and hydrogen peroxide-assisted proton transfers in guanine were studied and compared with the intramolecular proton transfer in guanine in the gas phase. It was found that the assisted proton transfers pass through about three times lower energy barriers than those found for isolated guanine tautomers. Figure DFT study of the gas phase proton transfer in guanine assisted by water, methanol and hydrogen peroxide     

DFT: a dynamic study of the interaction of ethanol and methanol with platinum

The interaction between alcohol molecules and platinum (Pt) was studied using molecular dynamics (MD; Born-Oppenheimer method). Alcohol molecules like ethanol and methanol present a similar molecular structure, with a methyl group (CH₃) at one end and a fragment of hydroxyl (OH) at the other. This fact generates two orientations that are considered in the interaction with Pt. The MD calculation results for these two orientations indicate a preferential orientation due to energy interactions. A plausible reaction mechanism that takes into account the interaction between Pt and alcohol is presented. The charge transference obtained from the Pt–alcohol interaction was also analyzed. The energy for the two orientations was calculated by indicating the preferential orientation. The methyl and hydroxyl groups are involved in heterolytic breakage of hydrogen bonds, joined to a carbon atom in the former and to an oxygen atom in the latter; however, the methyl group reaction seems to be the most important.     

DFT studies of the ionization of alpha and beta glycopyranosyl donors

Current attempts at mimicking the transition states (TSs) of glycosyl processing enzymes (GPEs) that proceed through TSs with a high degree of oxacarbenium ion formation suffer from a paucity of data about the conformations of such oxacarbenium ions. Because TSs are maxima, the current models based on minimized structures may need some refinement. As part of studies directed at optimizing chemical glycosylation the ionization of 3,4,6-tri-O-acetyl-alpha/beta-D-glucopyranosyl chlorides and triflates, 2,3,4,6-tetra-O-methyl-alpha/beta-D-glucopyranosyl fluorides, chlorides and triflates, 2,3,4,6-tetra-O-methyl-alpha/beta-D-mannopyranosyl fluorides, 2,3-di-O-methyl 4,6-O-benzylidene alpha/beta-D-mannopyranosyl triflates and 2,3-di-O-methyl 4,6-O-benzylidene alpha/beta-D-glucopyranosyl triflates was studied by a prototypic density functional theory (DFT) procedure. In all cases, the alpha-anomers ionized smoothly to 4H3 half chair conformations or adjacent envelopes. By contrast, all beta-anomers exhibited an abrupt conformational change before ionization was complete. The nature of the conformations sampled depends on both the leaving group and the protecting group. The methods presented can be readily adapted to the study of any GPE or chemical glycosylation and provide a method for initial evaluation of plausible TSs, which in turn can be used in mimetic design.     

Solution structure of molybdic acid from Raman spectroscopy and DFT analysis

Protonation of produces the well-characterized polymolybdates, but at concentrations below 10⁻³ M the dominant species is monomeric molybdic acid, H₂MoO₄. It is likely to be the species adsorbed on manganese oxide, a process thought to control levels in the ocean, because of the strong proton dependence of adsorption. The molecular structure of H₂MoO₄ is elusive, since it occurs only in dilute solutions. Using 244 nm laser excitation, near resonance with O → Mo charge-transfer electronic transitions of H₂MoO₄, we have detected a 919 cm⁻¹ Raman band assignable to ν_sMoO. Using DFT, we have computed geometries and vibrational modes for the various structures consistent with the H₂MoO₄ formula. We tested the computations on a series of Mo(VI) oxo complexes with known vibrational frequencies, at several levels of theory. Best agreement with experimental values, at reasonable computational cost, was obtained with the B3LYP functional, employing a LANL2DZ ECP basis set for Mo and the 6-311+G(2df,p) basis set for O and H. Among the possible H₂MoO₄ structures only those based on the MoO₃ unit, with one, two or three coordinated water molecules, gave a scaled frequency for ν_sMoO that was within two standard deviations of 919 cm⁻¹. Best agreement was obtained for MoO₃(H₂O)₃. The MoO₂ and MoO structures gave frequencies that were too high. The Mo(OH)₆ structure could be excluded, because its vibrational frequencies shift down strongly upon H/D exchange, whereas the 919 cm⁻¹H₂MoO₄ band shifts up 1 cm⁻¹ in D₂O.     

The metal bonding domain of the antitumor drug Fe(II)-bleomycin: a DFT investigation

The geometric and electronic structure of ferrous complexes of bleomycin (Fe(II)BLM) has been investigated by means of density functional theory (DFT) calculations. The active site of this antitumor drug is a highly distorted octahedral complex, with the coordination sphere completed by the five known endogenous ligands, including pyrimidine, imidazole, deprotonated amide, and secondary and primary amines. We have addressed the controversial issue of the nature of the sixth axial ligand, which we have identified as the oxygen of the carbamoyl group. Our conclusions are further validated by a comparison with structural data derived from NMR experiments. Moreover, because of the high sensitivity of structural data on the pH of the environment, we have investigated the effect of a different protonation state of the histidine amide on the geometric structure of the Fe(II)BLM complex. The extensive model of the active site of bleomycin considered in this work allows us to check the limitations of previous investigations based on simplified models.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.An erratum to this article can be found at     

Synthesis, DFT calculations, linear and nonlinear optical properties of binuclear phthalocyanine gallium chloride

The axially substituted binuclear GaCl/GaCl phthalocyanine 1 with an unsymmetrical pattern of substitution has been prepared and its nonlinear optical (NLO) properties determined. The resulting binuclear complex retains approximately the same transition energies of monomeric (RO)₈PcGaCl as far as the linear optical spectrum is concerned, although 1 has a double concentration of central atoms per molecule and an enlarged conjugated ligand. The lack of significant spectral shifts in passing from mononuclear to binuclear complexes has been rationalized theoretically by means of density functional theory calculations. The purpose of the present study is to determine whether binuclearity affects the optical limiting behavior of 1 with respect to monomeric (RO)₈PcGaCl in the NLO regime determined by nanosecond laser pulses. Figure Bis axially substituted binuclear phthalocyanine: synthesis, DFT calculations and NLO properties Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday.     

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

Investigation of the inhibition pathway of glucosamine synthase by N3-(4-methoxyfumaroyl)-L-2,3-diaminopropanoic acid by semiempirical quantum mechanical and molecular mechanics methods

Glucosamine synthase (E.C. 2.6.1.16) is a promising target in antifungal drug design. It has been reported that its potent inhibitor, N³-(4-methoxyfu-maroyl)-L-2,3-diaminopropanoic acid (FMDP), inactivates the enzyme by the Michael addition of the S-H group to the FMDP molecule followed by cyclisation reactions. In this study we have investigated, by means of semiempirical MNDO, PM3 and molecular mechanics methods, the energetics and kinetic possibility of the formation of various stereoisomers of the products of cyclisation of the Michael addition products detected experimentally. It was found that the substituted 1,4-thiazin-3-one can be formed in one step under alkaline conditions; the stereoisomers of this compound predicted to be the most stable on the basis of theoretical calculations are also the dominant ones in reality.Abbreviations FMDP N3-(4-methoxyfumaroyl)-L-2,3-diaminopropanoic acid - MNDO Modified Neglect of Diatomic Overlap - AM1 Advanced Method 1 - PM3 (full abbreviation: MNDOPM3) Modified Neglect of Diatomic Overlap, Parametrised Method 3 Correspondence to: A. Liwo     

The Interplay of Angle Strain and Aromaticity: Molecular and Electronic Structures of [0n]Paracyclophanes

The belt-like polyphenylenes, [0_n]paracyclophanes, (n = 5 and 6), have been investigated using semi-empirical, ab initio and DFT methods. The molecular structure, rotational barrier on twisting a single phenyl ring and the aromatic character within each ring as well as in the whole molecule have been evaluated. [0₅]Paracyclophane is predicted to have a quinonoid structure. In contrast, the equatorial pentaphenyl fragment found in C₇₀ as well as the hexagons of the less strained [0₆]paracyclophane have benzenoid character. Approximate band structures have been derived for larger cycles of [0_n] paracyclophanes.Electronic Supplementary Material available.     

Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements

Several modifications that have been made to the NDDO core-core interaction term and to the method of parameter optimization are described. These changes have resulted in a more complete parameter optimization, called PM6, which has, in turn, allowed 70 elements to be parameterized. The average unsigned error (AUE) between calculated and reference heats of formation for 4,492 species was 8.0 kcal mol⁻¹. For the subset of 1,373 compounds involving only the elements H, C, N, O, F, P, S, Cl, and Br, the PM6 AUE was 4.4 kcal mol⁻¹. The equivalent AUE for other methods were: RM1: 5.0, B3LYP 6–31G*: 5.2, PM5: 5.7, PM3: 6.3, HF 6–31G*: 7.4, and AM1: 10.0 kcal mol⁻¹. Several long-standing faults in AM1 and PM3 have been corrected and significant improvements have been made in the prediction of geometries. Figure Calculated structure of the complex ion [Ta₆Cl₁₂]²⁺ (footnote): Reference value in parenthesis Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Adsorption of CO on Cu (110) and (100) surfaces using COSMO-based DFT

Density functional theory (DFT) combined with the conductor-like solvent model (COSMO) can provide valuable atomistic level insights into CO adsorption on Cu surface interactions in liquid paraffin. The objective of this research was to investigate the solvent effect of liquid paraffin. It was found that both structural parameters and relative energies are very sensitive to the COSMO solvent model. Solvent effects can improve the stability of CO adsorption on Cu (110) and (100) surfaces and the extent of CO activation.     

DFT study of the intrinsic conformations of [2Fe-2S-4(SCH3)] clusters and their influence on exchange coupling

We have determined the equilibrium conformations of the diiron(III) cluster [2Fe-2S-4(SCH₃)]²⁻ using density functional theory. The conformers have dihedral Fe-Fe-S-C angles of ∼0° and ±120°. The relative energies of the conformers can be accurately parameterized with a small number of side-chain repulsion parameters. Of the 17 conformers identified on the basis of the ideal values for the dihedrals, 10 conformers are stable in both the ferromagnetic and broken symmetry state for the cluster. The exchange coupling constants for the seven energetically lowest conformers are predicted to belong to a narrow range, 150 cm⁻¹ ? J ? 178 cm⁻¹. The cluster conformers found in proteins do not coincide with any of the intrinsic ones, due to distortion of one of the dihedral angles under the influence of the protein scaffold.