The 2-pyridone/2-hydroxypyridine tautomerism in gas phase: An excited state roaming reaction.

Multiconfigurational methods (CASSCF and CASPT2) were employed to gain a new understanding of the mechanism of the gas-phase phototautomerization of 2-pyridone/2-hydroxypyridine. Potential energy curves and crossing points of the low-lying excited states were analyzed. The results show that the tautomerization only occurs from 2-pyridone to 2-hydroxypyridine after electronic excitation to the S1 (ππ*) state. From this state, the system would be able to reach a conical intersection between S1 and the dissociative S2 (πσ*) due to vibrational effects. Then, it can evolve to the hydroxy form in its ground state by reaching an intersection seam between the S0 and the πσ* states. For this to happen, a roaming process responsible for the hydrogen atom migration would be required; otherwise, the system would revert to the 2-pyridone tautomer. The unfeasibility of the reverse process after optical excitation from the lactim to the lactam form is explained by the great amount of energy needed to reach the conical intersection between the ππ* and πσ* states. These findings would provide new insights into the understanding of the photophysics and photochemistry of a primordial heterocycle, considered a prebiotic model known to be found in interstellar clouds.

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis.

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations.

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

Open questions on toxic heavy metals Cd, Hg and Pb binding small components of DNA and nucleobases. Are there any predictable trends?

In this perspective article, we provide a bibliographic compilation of experimental and theoretical work on Cd, Hg, and Pb, and analyze in detail the bonding of M2+ and CH3M+ (M = Zn, Cd, Hg, Pb) with urea and thiourea as suitable models for larger biochemical bases. Through the use of DFT calculations, we have found that although in principle binding energies decrease according to ionic size (Zn2+ > Cd2+ > Pb2+), Hg2+ largely breaks the trend. Through the use of EDA (Energy Decomposition Analysis) it is possible to explain this behavior, which is essentially due to the strong contribution of polarization to the binding. This conclusion is ratified by the NEDA (Natural Energy Decomposition Analysis) formalism, showing that the charge transfer term is very large in all cases, but particularly in the case of the mercury-thiourea system. The general trends observed for the interactions with CH3M+ monocations show however CH3Hg+ binding energies systematically smaller than the CH3Zn+ ones, likely because the relativistic contraction of the Hg orbitals is very much attenuated by the attachment to the methyl group. Finally, we have investigated the gas-phase reactivity between EtHg+ and uracil to compare it with that exhibited by CH3Hg+ and n-ButHg+ previously described in the literature. This comparison gathers new information that highlights the importance of the length of the alkyl chain attached to the metal on the mechanisms of these reactions. For methyl mercury, only the alkyl transfer process is allowed; for butyl mercury, protonation is clearly favored, and for ethyl mercury, both paths are competitive experimentally.

Characterization of cisplatin/membrane interactions by QM/MM energy decomposition analysis.

We extend for the first time a quantum mechanical energy decomposition analysis scheme based on deformation electron densities to a hybrid electrostatic embedding quantum mechanics/molecular mechanics framework. The implemented approach is applied to characterize the interactions between cisplatin and a dioleyl-phosphatidylcholine membrane, which play a key role in the permeation mechanism of the drug inside the cells. The interaction energy decomposition into electrostatic, induction, dispersion and Pauli repulsion contributions is performed for ensembles of geometries to account for conformational sampling. It is evidenced that the electrostatic and repulsive components are predominant in both polar and non-polar regions of the bilayer.

Tracking the Transition from Pericyclic to Pseudopericyclic Reaction Mechanisms Using Multicenter Electron Delocalization Analysis: The [1,3] Sigmatropic Rearrangement.

Herein, the power of multicenter electron delocalization analysis to elucidate the intricacies of concerted reaction mechanisms is brought to light by tracking the transition of [1,3] sigmatropic rearrangements from the high-barrier pericyclic mechanism in 1-butene to the barrierless pseudopericyclic mechanism in 1,2-diamino-1-nitrosooxyethane. This transition has been progressively achieved by substituting the migrating group, changing the donor and acceptor atoms, and functionalizing the alkene unit with weak and strong electron-donating and electron-withdrawing groups. Fourteen [1,3] sigmatropic reactions with electronic energy barriers ranging from 1 to 89 kcal/mol have been investigated. A very good correlation has been found between the barrier and the four-center electron delocalization at the transition state, the latter calculated for the atoms involved in the four-centered ring adduct formed along the reaction path. Surprisingly, the barrier has been found to be independent of the bond strength between the migrating group and the donor atom so that only the changes induced in the multicenter bonding control the kinetics of the reaction. Additional insights into the effect of atom substitution and group functionalization have also been extracted from the analysis of the multicenter electron delocalization profiles along the reaction path and qualitatively supported by the topological analysis of the electron density.