Molecular second-quantized Hamiltonian: Electron correlation and non-adiabatic coupling treated on an equal footing.

We introduce a new theoretical and computational framework for treating molecular quantum mechanics without the Born-Oppenheimer approximation. The molecular wavefunction is represented in a tensor-product space of electronic and vibrational basis functions, with electronic basis chosen to reproduce the mean-field electronic structure at all geometries. We show how to transform the Hamiltonian to a fully second-quantized form with creation/annihilation operators for electronic and vibrational quantum particles, paving the way for polynomial-scaling approximations to the tensor-product space formalism. In addition, we make a proof-of-principle application of the new Ansatz to the vibronic spectrum of C2.

The Molpro quantum chemistry package.

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Novel long-chain neurotoxins from Bungarus candidus distinguish the two binding sites in muscle-type nicotinic acetylcholine receptors.

αδ-Bungarotoxins, a novel group of long-chain α-neurotoxins, manifest different affinity to two agonist/competitive antagonist binding sites of muscle-type nicotinic acetylcholine receptors (nAChRs), being more active at the interface of α-δ subunits. Three isoforms (αδ-BgTx-1-3) were identified in Malayan Krait (Bungarus candidus) from Thailand by genomic DNA analysis; two of them (αδ-BgTx-1 and 2) were isolated from its venom. The toxins comprise 73 amino acid residues and 5 disulfide bridges, being homologous to α-bungarotoxin (α-BgTx), a classical blocker of muscle-type and neuronal α7, α8, and α9α10 nAChRs. The toxicity of αδ-BgTx-1 (LD50 = 0.17-0.28 µg/g mouse, i.p. injection) is essentially as high as that of α-BgTx. In the chick biventer cervicis nerve-muscle preparation, αδ-BgTx-1 completely abolished acetylcholine response, but in contrast with the block by α-BgTx, acetylcholine response was fully reversible by washing. αδ-BgTxs, similar to α-BgTx, bind with high affinity to α7 and muscle-type nAChRs. However, the major difference of αδ-BgTxs from α-BgTx and other naturally occurring α-neurotoxins is that αδ-BgTxs discriminate the two binding sites in the Torpedo californica and mouse muscle nAChRs showing up to two orders of magnitude higher affinity for the α-δ site as compared with α-ε or α-γ binding site interfaces. Molecular modeling and analysis of the literature provided possible explanations for these differences in binding mode; one of the probable reasons being the lower content of positively charged residues in αδ-BgTxs. Thus, αδ-BgTxs are new tools for studies on nAChRs.

Direct quantum dynamics using variational Gaussian wavepackets and Gaussian process regression.

The method of direct variational quantum nuclear dynamics in a basis of Gaussian wavepackets, combined with the potential energy surfaces fitted on-the-fly using Gaussian process regression, is described together with its implementation. Enabling exact and efficient analytic evaluation of Hamiltonian matrix elements, this approach allows for black-box quantum dynamics of multidimensional anharmonic molecular systems. Example calculations of intra-molecular proton transfer on the electronic ground state of salicylaldimine are provided, and future algorithmic improvements as well as the potential for multiple-state non-adiabatic dynamics are discussed.

Electron and nuclear dynamics following ionisation of modified bismethylene-adamantane.

We have simulated the coupled electron and nuclear dynamics using the Ehrenfest method upon valence ionisation of modified bismethylene-adamantane (BMA) molecules where there is an electron transfer between the two π bonds. We have shown that the nuclear motion significantly affects the electron dynamics after a few fs when the electronic states involved are close in energy. We have also demonstrated how the non-stationary electronic wave packet determines the nuclear motion, more precisely the asymmetric stretching of the two π bonds, illustrating "charge-directed reactivity". Taking into account the nuclear wave packet width results in the dephasing of electron dynamics with a half-life of 8 fs; this eventually leads to the equal delocalisation of the hole density over the two methylene groups and thus symmetric bond lengths.

A complete description of tunnelling using direct quantum dynamics simulation: Salicylaldimine proton transfer.

We demonstrate here conclusively that the variational multiconfiguration Gaussian (vMCG) method converges to the grid based full quantum dynamics multiconfiguration time-dependent Hartree result for a tunnelling problem in many dimensions, using the intramolecular proton transfer in salicylaldimine as a model system. The 13-dimensional model potential energy surface was obtained from Hartree Fock energies with the 6-31G* basis set and the expectation value of the flux operator along the transition mode was used as a benchmark characteristic. As well as showing excellent convergence of the vMCG method on the model surface using a local harmonic approximation and a moderate number of basis functions, we show that the direct dynamics version of the vMCG also performs very well, usually needs the same number of Gaussians to converge, and converges to exact results if those are obtained on an accurately fitted surface. Finally, we make an important observation that the width of the Gaussian basis functions must be chosen very carefully to obtain accurate results with the use of the frozen-width approximation.

The cinchona primary amine-catalyzed asymmetric epoxidation and hydroperoxidation of α,ß-unsaturated carbonyl compounds with hydrogen peroxide.

Using cinchona alkaloid-derived primary amines as catalysts and aqueous hydrogen peroxide as the oxidant, we have developed highly enantioselective Weitz-Scheffer-type epoxidation and hydroperoxidation reactions of α,ß-unsaturated carbonyl compounds (up to 99.5:0.5 er). In this article, we present our full studies on this family of reactions, employing acyclic enones, 5-15-membered cyclic enones, and α-branched enals as substrates. In addition to an expanded scope, synthetic applications of the products are presented. We also report detailed mechanistic investigations of the catalytic intermediates, structure-activity relationships of the cinchona amine catalyst, and rationalization of the absolute stereoselectivity by NMR spectroscopic studies and DFT calculations.

A microiterative intrinsic reaction coordinate method for large QM/MM systems.

Intrinsic reaction coordinate (IRC) computations are a valuable tool in theoretical studies of chemical reactions, but they can usually not be applied in their current form to handle large systems commonly described by quantum mechanics/molecular mechanics (QM/MM) methods. We report on a development that tackles this problem by using a strategy analogous to microiterative transition state optimization. In this approach, the IRC equations only govern the motion of a core region that contains at least the atoms directly involved in the reaction, while the remaining degrees of freedom are relaxed after each IRC step. This strategy can be used together with any existing IRC procedure. The present implementation covers the stabilized Euler, local quadratic approximation, and Hessian predictor-corrector algorithms for IRC calculations. As proof of principle, we perform tests at the QM level on small gas-phase systems and validate the results by comparisons with standard IRC procedures. The broad applicability of the method is demonstrated by IRC computations for two enzymatic reactions using standard QM/MM setups.

Theoretical rotation-vibration spectrum of thioformaldehyde.

We present a variational calculation of the first comprehensive T = 300 K rovibrational line list for thioformaldehyde, H2CS. It covers 41,809 rovibrational levels for states up to J(max) = 30 with vibrational band origins up to 5000 cm(-1) and provides the energies and line intensities for 547,926 transitions from the ground vibrational state to these levels. It is based on our previously reported accurate ab initio potential energy surface and a newly calculated ab initio dipole moment surface. Minor empirical adjustments are made to the ab initio equilibrium geometry to reduce systematic errors in the predicted intra-band rotational energy levels. The rovibrational energy levels and transition intensities are computed variationally by using the methods implemented in the computer program TROVE. Transition wavelengths and intensities are found to be in excellent agreement with the available experimental data. The present calculations correctly reproduce the observed resonance effects, such as intensity borrowing, thus reflecting the high accuracy of the underlying ab initio surfaces. We report a detailed analysis of several vibrational bands, especially those complicated by strong Coriolis coupling, to facilitate future laboratory assignments.

Quantum mechanics/molecular mechanics dual Hamiltonian free energy perturbation.

The dual Hamiltonian free energy perturbation (DH-FEP) method is designed for accurate and efficient evaluation of the free energy profile of chemical reactions in quantum mechanical/molecular mechanical (QM/MM) calculations. In contrast to existing QM/MM FEP variants, the QM region is not kept frozen during sampling, but all degrees of freedom except for the reaction coordinate are sampled. In the DH-FEP scheme, the sampling is done by semiempirical QM/MM molecular dynamics (MD), while the perturbation energy differences are evaluated from high-level QM/MM single-point calculations at regular intervals, skipping a pre-defined number of MD sampling steps. After validating our method using an analytic model potential with an exactly known solution, we report a QM/MM DH-FEP study of the enzymatic reaction catalyzed by chorismate mutase. We suggest guidelines for QM/MM DH-FEP calculations and default values for the required computational parameters. In the case of chorismate mutase, we apply the DH-FEP approach in combination with a single one-dimensional reaction coordinate and with a two-dimensional collective coordinate (two individual distances), with superior results for the latter choice.

Quantum mechanical/molecular mechanical study on the mechanism of the enzymatic Baeyer-Villiger reaction.

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP(+), while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme-reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP(+) cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP(+) cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step.

Retaining glycosyltransferase mechanism studied by QM/MM methods: lipopolysaccharyl-α-1,4-galactosyltransferase C transfers α-galactose via an oxocarbenium ion-like transition state.

Glycosyltransferases (GTs) catalyze the highly specific biosynthesis of glycosidic bonds and, as such, are important both as drug targets and for biotechnological purposes. Despite their broad interest, fundamental questions about their reaction mechanism remain to be answered, especially for those GTs that transfer the sugar with net retention of the configuration at the anomeric carbon (retaining glycosyltransferases, ret-GTs). In the present work, we focus on the reaction catalyzed by lipopolysaccharyl-α-1,4-galactosyltransferase C (LgtC) from Neisseria meningitides. We study and compare the different proposed mechanisms (S(N)i, S(N)i-like, and double displacement mechanism via a covalent glycosyl-enzyme intermediate, CGE) by using density functional theory (DFT) and quantum mechanics/molecular mechanics (QM/MM) calculations on the full enzyme. We characterize a dissociative single-displacement (S(N)i) mechanism consistent with the experimental data, in which the acceptor substrate attacks on the side of the UDP leaving group that acts as a catalytic base. We identify several key interactions that help this front-side attack by stabilizing the transition state. Among them, Gln189, the putative nucleophile in a double displacement mechanism, is shown to favor the charge development at the anomeric center by about 2 kcal/mol, compatible with experimental mutagenesis data. We predict that using 3-deoxylactose as acceptor would result in a reduction of k(cat) to 0.6-3% of that for the unmodified substrates. The reactions of the Q189A and Q189E mutants have also been investigated. For Q189E, there is a change in mechanism since a CGE can be formed which, however, is not able to evolve to products. The current findings are discussed in the light of the available experimental data and compared with those for other ret-GTs.

The phenoxyl radical-water complex--a matrix isolation and computational study.

The phenoxyl radical 1 was generated in high yields by flash vacuum pyrolysis of allyl phenyl ether 2 with subsequent trapping of the products in argon at 3 K. In water-doped argon matrices, an OH···O complex between 1 and water is formed that could be characterized by IR spectroscopy. Several isotopomers of the complex were generated, and the IR spectra compared to results of density functional theory calculations. Other dimers between 1 and water were not found under these conditions. QM/MM calculations in simulated argon matrices reveal that an OH···π complex is unstable even at a time scale of picoseconds. This finding has implications on the related interaction between the tyrosyl radical and the water in biological systems.

Ultrafast Photoinduced Dynamics of 1,3-Cyclohexadiene Using XMS-CASPT2 Surface Hopping.

A full-dimensional simulation of the photodissociation of 1,3-cyclohexadiene in the manifold of three electronic states was performed via nonadiabatic surface hopping dynamics using extended multistate complete active space second-order perturbation (XMS-CASPT2) electronic structure theory with fully analytic nonadiabatic couplings. With the 47 ± 8% product quantum yield calculated from the 136 trajectories, generally 400 fs-long, and an estimated excited lifetime of 89 ± 9 fs, our calculations provide a detailed description of the nonadiabatic deactivation mechanism, showing the existence of an extended conical intersection seam along the reaction coordinate. The nature of the preferred reaction pathways on the ground state is discussed and extensive comparison to the previously published full dimensional dynamics calculations is provided.

Toward Accurate QM/MM Reaction Barriers with Large QM Regions Using Domain Based Pair Natural Orbital Coupled Cluster Theory.

The hydroxylation reaction catalyzed by p-hydroxybenzoate hydroxylase and the Baeyer-Villiger reaction catalyzed by cyclohexanone monooxygenase are investigated by means of quantum mechanical/molecular mechanical (QM/MM) calculations at different levels of QM theory. The geometries of the stationary points along the reaction profile are obtained from QM/MM geometry optimizations, in which the QM region is treated by density functional theory (DFT). Relative energies are determined from single-point QM/MM calculations using the domain-based local pair natural orbital coupled cluster DLPNO-CCSD(T) method as QM component. The results are compared with single-point DFT/MM energies obtained using popular density functionals and with available experimental and computational data. It is found that the choice of the QM method strongly affects the computed energy profiles for these reactions. Different density functionals provide qualitatively different energy barriers (variations of the order of 10 kcal/mol in both reactions), thus limiting the confidence in DFT/MM computational predictions of energy profiles. On the other hand, the use of the DLPNO-CCSD(T) method in conjunction with large QM regions and basis sets makes it possible to achieve high accuracy. A critical discussion of all the technical aspects of the calculations is given with the aim of aiding computational chemists in the application of the DLPNO-CCSD(T) methodology in QM/MM calculations.

Quantum mechanical/molecular mechanical study on the enantioselectivity of the enzymatic Baeyer-Villiger reaction of 4-hydroxycyclohexanone.

We report a combined quantum mechanical/molecular mechanical (QM/MM) study of the effect of mutations of the Phe434 residue in the active site of cyclohexanone monooxygenase (CHMO) on its enantioselectivity toward 4-hydroxycyclohexanone. In terms of our previously established model of the enzymatic Baeyer-Villiger reaction, enantioselectivity is governed by the preference toward the equatorial ((S)-selectivity) or axial ((R)-selectivity) conformation of the substituent at the C4 carbon atom of the cyclohexanone ring in the Criegee intermediate and the subsequent rate-limiting transition state for migration (TS2). We assess the enantiopreference by locating all relevant TS2 structures at the QM/MM level. In the wild-type enzyme we find that the axial conformation is energetically slightly more stable, thus leading to a small excess of (R)-product. In the Phe434Ser mutant, there is a hydrogen bond between the serine side chain and the equatorial substrate hydroxyl group that is retained during the whole reaction, and hence there is pronounced reverse (S)-enantioselectivity. Another mutant, Phe434Ile, is shown to preserve and enhance the (R)-selectivity. All these findings are in accordance with experiment. The QM/MM calculations allow us to explain the effect of point mutations on CHMO enantioselectivity for the first time at the molecular level by an analysis of the specific interactions between substrate and active-site environment in the TS2 structures that satisfy the basic stereoelectronic requirement of anti-periplanarity for the migrating σ-bond.