T. Verstraelen

In this paper the authors develop a method to accurately calculate localized vibrational modes for partially optimized molecular structures or for structures containing link atoms. The method avoids artificially introduced imaginary frequencies and keeps track of the invariance under global translations and rotations. Only a subblock of the Hessian matrix has to be constructed and diagonalized, leading to a serious reduction of the computational time for the frequency analysis. The mobile block Hessian approach (MBH) proposed in this work can be regarded as an extension of the partial Hessian vibrational analysis approach proposed by Head [Int. J. Quantum Chem. 65, 827 (1997)] . Instead of giving the nonoptimized region of the system an infinite mass, it is allowed to move as a rigid body with respect to the optimized region of the system. The MBH approach is then extended to the case where several parts of the molecule can move as independent multiple rigid blocks in combination with single atoms. The merits of both models are extensively tested on ethanol and di-n-octyl-ether.

Ring opening of 1-arylmethyl-2-(cyanomethyl)aziridines with HBr afforded 3-(arylmethyl)amino-4-bromobutyronitriles via regiospecific ring opening at the unsubstituted aziridine carbon. Previous experimental and theoretical reports show treatment of the same compounds with benzyl bromide to furnish 4-amino-3-bromobutanenitriles through ring opening at the substituted aziridine carbon. To gain insights into the regioselective preference with HBr, reaction paths have been analyzed with computational methods. The effect of solvation was taken into account by the use of explicit solvent molecules. Geometries were determined at the B3LYP/6-31++G(d,p) level of theory, and a Grimme-type correction term was included for long-range dispersion interactions; relative energies were refined with the meta-hybrid MPW1B95 functional. Activation barriers confirm preference for ring opening at the unsubstituted ring carbon for HBr. HBr versus benzyl bromide ring opening was analyzed through comparison of the electronic structure of corresponding aziridinium intermediates. Although the electrostatic picture fails to explain the opposite regiospecific nature of the reaction, frontier molecular orbital analysis of LUMOs and nucleophilic Fukui functions show a clear preference of attack for the substituted aziridine carbon in the benzyl bromide case and for the unsubstituted aziridine carbon in the HBr case, successfully rationalizing the experimentally observed regioselectivity.

The neutral and anionic semiquinone radicals of the flavin adenine dinucleotide (FAD) cofactor noncovalently bound in glucose oxidase from A. niger are examined with the aid of QM/MM molecular modeling methods, enabling complete inclusion of the protein environment. Recently, the electron paramagnetic resonance (EPR) characteristics, the anisotropic g tensor and all the significant hyperfine couplings, of these flavoprotein radicals were determined at high resolution (J. Phys. Chem. B 2008, 112, 3568). A striking difference between the neutral and anionic radical forms was found to be a shift in the gy principal value. Within the QM/MM framework, geometry optimization and molecular dynamics simulations are combined with EPR property calculations, employing a recent implementation by some of the authors in the CP2K software package. In this way, spectroscopic characteristics are computed on the fly during the MD simulations of the solvated protein structure, mimicking as best as possible the experimental conditions. The general agreement between calculated and experimental EPR properties is satisfactory and on par with those calculated with other codes (Gaussian 03, ORCA). The protonation state of two histidines (His559 and His516) at the catalytic site of this flavoprotein is found to have a remarkable influence on the isotropic hyperfine coupling of one of the methyl groups on the neutral FAD radical, which is consistent with experimental findings in other flavoproteins (J. Biol. Chem. 2007, 282, 4738). Furthermore, the shift in the gy principal values between the neutral and anionic radicals is well reproduced by QM/MM simulations. Incorporation of at least the nearest protein environment of the cofactor radicals proves to be vital for a correct reproduction, indicating that this shift is a global feature of the protein rather than a local one. In addition, QM/MM techniques are used to make a prediction of relative angles between important spectroscopic principal directions, which are not readily determined by conventional EPR experiments. Significantly, the directions of the gx and the gy components of the g-tensor that lie in the plane of the isoalloxazine moiety are rotated by approximately 59° between the neutral and the anionic radicals.

TAMkin is a program for the calculation and analysis of normal modes, thermochemical properties and chemical reaction rates. At present, the output from the frequently applied software programs ADF, CHARMM, CPMD, CP2K, Gaussian, Q-Chem, and VASP can be analyzed. The normal-mode analysis can be performed using a broad variety of advanced models, including the standard full Hessian, the Mobile Block Hessian, the Partial Hessian Vibrational approach, the Vibrational Subsystem Analysis with or without mass matrix correction, the Elastic Network Model, and other combinations. TAMkin is readily extensible because of its modular structure. Chemical kinetics of unimolecular and bimolecular reactions can be analyzed in a straightforward way using conventional transition state theory, including tunneling corrections and internal rotor refinements. A sensitivity analysis can also be performed, providing important insight into the theoretical error margins on the kinetic parameters. Two extensive examples demonstrate the capabilities of TAMkin: the conformational change of the biological system adenylate kinase is studied, as well as the reaction kinetics of the addition of ethene to the ethyl radical. The important feature of batch processing large amounts of data is highlighted by performing an extended level of theory study, which TAMkin can automate significantly.