We present an overview of recent developments in the computational chemistry environment ChemShell for multiscale QM/MM modelling of biomolecular and materials catalysts, together with a survey of recent catalysis applications using ChemShell.

High-resolution crystal structures of enzymes in relevant redox states have transformed our understanding of enzyme catalysis. Recent developments have demonstrated that X-rays can be used, via the generation of solvated electrons, to drive reactions in crystals at cryogenic temperatures (100 K) to generate `structural movies' of enzyme reactions. However, a serious limitation at these temperatures is that protein conformational motion can be significantly supressed. Here, the recently developed MSOX (multiple serial structures from one crystal) approach has been applied to nitrite-bound copper nitrite reductase at room temperature and at 190 K, close to the glass transition. During both series of multiple structures, nitrite was initially observed in a `top-hat' geometry, which was rapidly transformed to a `side-on' configuration before conversion to side-on NO, followed by dissociation of NO and substitution by water to reform the resting state. Density functional theory calculations indicate that the top-hat orientation corresponds to the oxidized type 2 copper site, while the side-on orientation is consistent with the reduced state. It is demonstrated that substrate-to-product conversion within the crystal occurs at a lower radiation dose at 190 K, allowing more of the enzyme catalytic cycle to be captured at high resolution than in the previous 100 K experiment. At room temperature the reaction was very rapid, but it remained possible to generate and characterize several structural states. These experiments open up the possibility of obtaining MSOX structural movies at multiple temperatures (MSOX-VT), providing an unparallelled level of structural information during catalysis for redox enzymes.

Interactive molecular dynamics simulation in virtual reality (iMD-VR) is emerging as a promising technique in molecular science. Here, we demonstrate its use in a range of fifteen applications in materials science and heterogeneous catalysis. In this work, the iMD-VR package Narupa is used with the MD package, DL_POLY [1]. We show how iMD-VR can be used to: (i) investigate the mechanism of lithium fast ion conduction by directing the formation of defects showing that vacancy transport is favoured over interstitialcy mechanisms, and (ii) guide a molecule through a zeolite pore to explore diffusion within zeolites, examining in detail the motion of methyl n-hexanoate in H-ZSM-5 zeolite and identifying bottlenecks restricting diffusion. iMD-VR allows users to manipulate these systems intuitively, to drive changes in them and observe the resulting changes in structure and dynamics. We make these simulations available, as a resource for both teaching and research. All simulation files, with videos, can be found online (https://doi.org/10.5281/zenodo.8252314) and are provided as open-source material.

Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms. This article is part of a discussion meeting issue ‘Supercomputing simulations of advanced materials’.

Microbial nitrite reductases are denitrifying enzymes that are a major component of the global nitrogen cycle. Multiple structures measured from one crystal (MSOX data) of copper nitrite reductase at 240 K, together with molecular-dynamics simulations, have revealed protein dynamics at the type 2 copper site that are significant for its catalytic properties and for the entry and exit of solvent or ligands to and from the active site. Molecular-dynamics simulations were performed using different protonation states of the key catalytic residues (Asp_CATand His_CAT) involved in the nitrite-reduction mechanism of this enzyme. Taken together, the crystal structures and simulations show that the Asp_CATprotonation state strongly influences the active-site solvent accessibility, while the dynamics of the active-site `capping residue' (Ile_CAT), a determinant of ligand binding, are influenced both by temperature and by the protonation state of Asp_CAT. A previously unobserved conformation of Ile_CATis seen in the elevated temperature series compared with 100 K structures. DFT calculations also show that the loss of a bound water ligand at the active site during the MSOX series is consistent with reduction of the type 2 Cu atom.