Mercury Methylation by Cobalt Corrinoids: Relativistic Effects Dictate the Reaction Mechanism

Scalar Relativistic Effects with Multiwavelets: Implementation and Benchmark

The importance of relativistic effects in quantum chemistry is widely recognized, not only for heavier elements but throughout the periodic table. At the same time, relativistic effects are strongest in the nuclear region, where the description of electrons through a linear combination of atomic orbitals becomes more challenging. Furthermore, the choice of basis sets for heavier elements is limited compared with lighter elements where precise basis sets are available. Thanks to the framework of multiresolution analysis, multiwavelets provide an appealing alternative to overcoming this challenge: they lead to robust error control and adaptive algorithms that automatically refine the basis set description until the desired precision is reached. This allows one to achieve a proper description of the nuclear region. In this work, we extended the multiwavelet-based code MRChem to the scalar zero-order regular approximation framework. We validated our implementation by comparing the total energies for a small set of elements and molecules. To confirm the validity of our implementation, we compared both against a radial numerical code for atoms and the plane-wave-based code EXCITING.

A general method for locating stationary points on the mixed-spin surface of spin-forbidden reaction with multiple spin states

Some chemical reactions proceed on multiple potential energy surfaces and are often accompanied by a change in spin multiplicity, being called spin-forbidden reactions, where the spin-orbit coupling (SOC) effects play a crucial role. In order to efficiently investigate spin-forbidden reactions with two spin states, Yang et al. [Phys. Chem. Chem. Phys. 20, 4129-4136 (2018)] proposed a two-state spin-mixing (TSSM) model, where the SOC effects between the two spin states are simulated by a geometry-independent constant. Inspired by the TSSM model, we suggest a multiple-state spin-mixing (MSSM) model in this paper for the general case with any number of spin states, and its analytic first and second derivatives have been developed for locating stationary points on the mixed-spin potential energy surface and estimating thermochemical energies. To demonstrate the performance of the MSSM model, some spin-forbidden reactions involving 5d transition elements are calculated using the density functional theory (DFT), and the results are compared with the two-component relativistic ones. It is found that MSSM DFT and two-component DFT calculations may provide very similar stationary-point information on the lowest mixed-spin/spinor energy surface, including structures, vibrational frequencies, and zero-point energies. For the reactions containing saturated 5d elements, the reaction energies by MSSM DFT and two-component DFT agree very well within 3 kcal/mol. As for the two reactions OsO+ + CH4 → OOs(CH2)+ + H2 and W + CH4 → WCH2 + H2 involving unsaturated 5d elements, MSSM DFT may also yield good reaction energies of similar accuracy but with some counterexamples. Nevertheless, the energies may be remarkably improved by a posteriori single point energy calculations using two-component DFT at the MSSM DFT optimized geometries, and the maximum error of about 1 kcal/mol is almost independent of the SOC constant used. The MSSM method as well as the developed computer program provides an effective utility for studying spin-forbidden reactions.

The inorganic chemistry of the cobalt corrinoids - an update

The inorganic chemistry of the cobalt corrinoids, derivatives of vitamin B₁₂, is reviewed, with particular emphasis on equilibrium constants for, and kinetics of, their axial ligand substitution reactions. The role the corrin ligand plays in controlling and modifying the properties of the metal ion is emphasised. Other aspects of the chemistry of these compounds, including their structure, corrinoid complexes with metals other than cobalt, the redox chemistry of the cobalt corrinoids and their chemical redox reactions, and their photochemistry are discussed. Their role as catalysts in non-biological reactions and aspects of their organometallic chemistry are briefly mentioned. Particular mention is made of the role that computational methods - and especially DFT calculations - have played in developing our understanding of the inorganic chemistry of these compounds. A brief overview of the biological chemistry of the B₁₂-dependent enzymes is also given for the reader's convenience.

Photoproduct formation in coenzyme B12-dependent CarH via a singlet pathway

The CarH photoreceptor exploits of the light-sensing ability of coenzyme B₁₂ ( adenosylcobalamin = AdoCbl) to perform its catalytic function, which includes large-scale structural changes to regulate transcription. In daylight, transcription is activated in CarH via the photo-cleavage of the Co-C_5' bond of coenzyme B₁₂. Subsequently, the photoproduct, 4',5'-anhydroadenosine (anhAdo) is formed inducing dissociation of the CarH tetramer from DNA. Several experimental studies have proposed that hydridocoblamin (HCbl) may be formed in process with anhAdo. The photolytic cleavage of the Co-C_5' bond of AdoCbl was previously investigated using photochemical techniques and the involvement of both singlet and triplet excited states were explored. Herein, QM/MM calculations were employed to probe (1) the photolytic processes which may involve singlet excited states, (2) the mechanism of anhAdo formation, and (3) whether HCbl is a viable intermediate in CarH. Time-dependent density functional theory (TD-DFT) calculations indicate that the mechanism of photodissociation of the Ado ligand involves the ligand field (LF) portion of the lowest singlet excited state (S₁) potential energy surface (PES). This is followed by deactivation to a point on the S₀ PES where the Co-C_5' bond remains broken. This species corresponds to a singlet diradical intermediate. From this point, the PES for anhAdo formation was explored, using the Co-C_5' and Co-C_4' bond distances as active coordinates, and a local minimum representing anhAdo and HCbl formation was found. The transition state (TS) for the formation of the Co-H bond of HCbl was located and its identity was confirmed by a single imaginary frequency of i1592 cm^-1. Comparisons to experimental studies and the potential role of rotation around the N-glycosidic bond of the Ado ligand were discussed.

Relativity as a Synthesis Design Principle: A Comparative Study of [3 + 2] Cycloaddition of Technetium(VII) and Rhenium(VII) Trioxo Complexes with Olefins

The difference in [3 + 2] cycloaddition reactivity between fac-[MO₃(tacn)]⁺ (M = Re, ⁹⁹Tc; tacn = 1,4,7-triazacyclononane) complexes has been reexamined with a selection of unsaturated substrates including sodium 4-vinylbenzenesulfonate, norbornene, 2-butyne, and 2-methyl-3-butyn-2-ol (2MByOH). None of the substrates was found to react with the Re cation in water at room temperature, whereas the ⁹⁹Tc reagent cleanly yielded the [3 + 2] cycloadducts. Interestingly, a bis-adduct was obtained as the sole product for 2MByOH, reflecting the high reactivity of a ⁹⁹TcO-enediolato monoadduct. On the basis of scalar relativistic and nonrelativistic density functional theory calculations of the reaction pathways, the dramatic difference in reactivity between the two metals has now been substantially attributed to differences in relativistic effects, which are much larger for the 5d metal. Furthermore, scalar-relativistic ΔG values were found to decrease along the series propene > norbornene > 2-butyne > dimethylketene, indicating major variations in the thermodynamic driving force as a function of the unsaturated substrate. The suggestion is made that scalar-relativistic effects, consisting of greater destabilization of the valence electrons of the 5d elements compared with those of the 4d elements, be viewed as a new design principle for novel ^99mTc/Re radiopharmaceuticals, as well as more generally in heavy-element coordination chemistry.

Room temperature cycloaddition reactivity of fac-[⁹⁹TcO₃(tacn)]⁺ (tacn = 1,4,7-triazacyclononane) with a variety of unsaturated substrates and the lack of such reactivity for fac-[ReO₃(tacn)]⁺ appears largely attributable to much stronger relativistic effects for Re relative to Tc, based on relativistic density functional theory calculations.

Structure determination of the HgcAB complex using metagenome sequence data: insights into microbial mercury methylation

Bacteria and archaea possessing the hgcAB gene pair methylate inorganic mercury (Hg) to form highly toxic methylmercury. HgcA consists of a corrinoid binding domain and a transmembrane domain, and HgcB is a dicluster ferredoxin. However, their detailed structure and function have not been thoroughly characterized. We modeled the HgcAB complex by combining metagenome sequence data mining, coevolution analysis, and Rosetta structure calculations. In addition, we overexpressed HgcA and HgcB in Escherichia coli, confirmed spectroscopically that they bind cobalamin and [4Fe-4S] clusters, respectively, and incorporated these cofactors into the structural model. Surprisingly, the two domains of HgcA do not interact with each other, but HgcB forms extensive contacts with both domains. The model suggests that conserved cysteines in HgcB are involved in shuttling Hg^II, methylmercury, or both. These findings refine our understanding of the mechanism of Hg methylation and expand the known repertoire of corrinoid methyltransferases in nature.

Connor J. Cooper et al. expressed HgcA and HgcB in Escherichia coli and modeled the structure of the HgcAB complex by combining metagenome sequence data, coevolution analysis, and ab initio structure calculations. This study provides insights into the biochemical mechanism of mercury (Hg) methylation.

Spin-orbit coupling and vibronic transitions of Ce(C3H4) and Ce(C3H6) formed by the Ce reaction with propene: Mass-analyzed threshold ionization and relativistic quantum computation

A Ce atom reaction with propene is carried out in a pulsed laser vaporization molecule beam source. Several Ce-hydrocarbon species formed by the C-H and C-C bond activation of propene are observed by time-of-flight mass spectrometry, and Ce(C₃H_n) (n = 4 and 6) are characterized by mass-analyzed threshold ionization (MATI) spectroscopy and density functional theory, multiconfiguration, and relativistic quantum chemical calculations. The MATI spectrum of each species consists of two vibronic band systems, each with several vibronic bands. Ce(C₃H₆) is identified as an inserted species with Ce inserting into an allylic C-H bond of propene and Ce(C₃H₄) as a metallocycle through 1,2-vinylic dehydrogenation. Both species have a C_s structure with the Ce 4f¹6s¹ ground valence electron configuration in the neutral molecule and the Ce 4f¹ configuration in the singly charged ion. The two vibronic band systems observed for each species are attributed to the ionization of two pairs of the lowest spin-orbit coupled states with each pair being nearly degenerate.

Evidence of the Elusive Gold-Induced Non-classical Hydrogen Bonding in Aqueous Environments

The ability of a gold ion to act as a proton acceptor in hydrogen bonding continues to remain an open question. Heavy-atom effects and secondary competitive interactions in gold complexes make it challenging to precisely establish the identity of gold-ion-induced hydrogen bonding via experimental techniques. In such situations, computational chemistry may play an important role. Herein we have performed Born-Oppenheimer molecular dynamics simulations to study the behavior of [Au(CH₃)₂)]^- in bulk and interfacial aqueous environments. The simulation results suggest that the [Au(CH₃)₂)]^- complex forms one and two gold-ion-induced hydrogen bonds with the water molecules in interfacial and bulk environments, respectively. The calculated probabilities of key hydrogen-bonded configurations of [Au(CH₃)₂)]^-, combined distribution functions, and diffusion coefficients further support this unusual hydrogen-bonding interaction. In summary, the present results suggest that gold-ion-induced hydrogen bonding in an actual solvent environment may be feasible. These results will improve our understanding about the role of weak interactions in transition metal complexes and, thus, will have implications in catalysis and supramolecular assemblies.

Spin-Forbidden Reactivity of Transition Metal Oxo Species: Exploring the Potential Energy Surfaces

Spin-forbidden reactions are frequently encountered when transition metal oxo species are involved, particularly in oxygen transfer reactivity. The computational study of such reactions is challenging, because reactants and products are located on different spin potential energy surfaces (PESs). One possible approach to describe these reactions is the so-called minimum energy crossing point (MECP) between the diabatic reactants and products PESs. Alternatively, inclusion of spin-orbit coupling (SOC) effects allows to locate a saddle point on a single adiabatic PES (TS SOC). The TS SOC approach is rarely applied because of its high computational cost. Recently evidence for a TS SOC impact on significantly lowering the activation barrier in dioxygen addition to a carbene-gold(I)-hydride complex reaction (Chem. Sci. 2016, 7, 7034-7039) or even on predicting a qualitatively different reaction mechanism in mercury methylation by cobalt corrinoid (Angew. Chem. Int. Ed. 2016, 55, 11503-11506) has been put forward. Using MECP and TS SOC approaches a systematic analysis is provided here of three prototypical transition metal oxo spin-forbidden processes to investigate their implications on reactivity. Cycloaddition of ethylene to chromyl chloride (CrO₂ Cl₂ +C₂ H₄ ), iron oxide cation insertion into the hydrogen molecule (FeO⁺ +H₂ ) and H-abstraction from toluene by a Mn^V -oxo-porphyrin cation (MnOP(H₂ O)⁺ +C₆ H₅ CH₃ ) are case studies. For all these processes the MECP and TS SOC results are compared, which show that the spin-forbidden reactivity of transition metal oxo species can be safely described by a MECP approach, at least for the first-row transition metals investigated here, where the spin-orbit coupling is relatively weak. However, for the Mn-oxo reactivity, the MECP and TS SOC have been found to be crucial for a correct description of the reaction mechanism. In particular, the TS SOC approach allows to straightforwardly explore detailed features of the adiabatic potential energy surface which in principle could affect the overall reaction rate in cases where the involved diabatic PESs are tricky.

Quantum Molecular Dynamics, Topological, Group Theoretical and Graph Theoretical Studies of Protein-Protein Interactions

BACKGROUND

Protein-protein interactions (PPIs) are becoming increasingly important as PPIs form the basis of multiple aggregation-related diseases such as cancer, Creutzfeldt-Jakob, and Alzheimer's diseases. This mini-review presents hybrid quantum molecular dynamics, quantum chemical, topological, group theoretical, graph theoretical, and docking studies of PPIs. We also show how these theoretical studies facilitate the discovery of some PPI inhibitors of therapeutic importance.

OBJECTIVE

The objective of this review is to present hybrid quantum molecular dynamics, quantum chemical, topological, group theoretical, graph theoretical, and docking studies of PPIs. We also show how these theoretical studies enable the discovery of some PPI inhibitors of therapeutic importance.

METHODS

This article presents a detailed survey of hybrid quantum dynamics that combines classical and quantum MD for PPIs. The article also surveys various developments pertinent to topological, graph theoretical, group theoretical and docking studies of PPIs and highlight how the methods facilitate the discovery of some PPI inhibitors of therapeutic importance.

RESULTS

It is shown that it is important to include higher-level quantum chemical computations for accurate computations of free energies and electrostatics of PPIs and Drugs with PPIs, and thus techniques that combine classical MD tools with quantum MD are preferred choices. Topological, graph theoretical and group theoretical techniques are shown to be important in studying large network of PPIs comprised of over 100,000 proteins where quantum chemical and other techniques are not feasible. Hence, multiple techniques are needed for PPIs.

CONCLUSION

Drug discovery and our understanding of complex PPIs require multifaceted techniques that involve several disciplines such as quantum chemistry, topology, graph theory, knot theory and group theory, thus demonstrating a compelling need for a multi-disciplinary approach to the problem.

Kinetics of Enzymatic Mercury Methylation at Nanomolar Concentrations Catalyzed by HgcAB

Methylmercury (MeHg) is a potent bioaccumulative neurotoxin that is produced by certain anaerobic bacteria and archaea. Mercury (Hg) methylation has been linked to the gene pair hgcAB, which encodes a membrane-associated corrinoid protein and a ferredoxin. Although microbial Hg methylation has been characterized in vivo, the cellular biochemistry and the specific roles of the gene products HgcA and HgcB in Hg methylation are not well understood. Here, we report the kinetics of Hg methylation in cell lysates of Desulfovibrio desulfuricans ND132 at nanomolar Hg concentrations. The enzymatic Hg methylation mediated by HgcAB is highly oxygen sensitive, irreversible, and follows Michaelis-Menten kinetics, with an apparent K_m of 3.2 nM and V _max of 19.7 fmol · min^-1 · mg^-1 total protein for the substrate Hg(II). Although the abundance of HgcAB in the cell lysates is extremely low, Hg(II) was quantitatively converted to MeHg at subnanomolar substrate concentrations. Interestingly, increasing thiol/Hg(II) ratios did not impact Hg methylation rates, which suggests that HgcAB-mediated Hg methylation effectively competes with cellular thiols for Hg(II), consistent with the low apparent K_m Supplementation of 5-methyltetrahydrofolate or pyruvate did not enhance MeHg production, while both ATP and a nonhydrolyzable ATP analog decreased Hg methylation rates in cell lysates under the experimental conditions. These studies provide insights into the biomolecular processes associated with Hg methylation in anaerobic bacteria.IMPORTANCE The concentration of Hg in the biosphere has increased dramatically over the last century as a result of industrial activities. The microbial conversion of inorganic Hg to MeHg is a global public health concern due to bioaccumulation and biomagnification of MeHg in food webs. Exposure to neurotoxic MeHg through the consumption of fish represents a significant risk to human health and can result in neuropathies and developmental disorders. Anaerobic microbial communities in sediments and periphyton biofilms have been identified as sources of MeHg in aquatic systems, but the associated biomolecular mechanisms are not fully understood. In the present study, we investigate the biochemical mechanisms and kinetics of MeHg formation by HgcAB in sulfate-reducing bacteria. These findings advance our understanding of microbial MeHg production and may help inform strategies to limit the formation of MeHg in the environment.

Environmental Mercury Chemistry - In Silico

Mercury (Hg) is a global environmental contaminant. Major anthropogenic sources of Hg emission include gold mining and the burning of fossil fuels. Once deposited in aquatic environments, Hg can undergo redox reactions, form complexes with ligands, and adsorb onto particles. It can also be methylated by microorganisms. Mercury, especially its methylated form methylmercury, can be taken up by organisms, where it bioaccumulates and biomagnifies in the food chain, leading to detrimental effects on ecosystem and human health. In support of the recently enforced Minamata Convention on Mercury, a legally binding international convention aimed at reducing the anthropogenic emission of-and human exposure to-Hg, its global biogeochemical cycle must be understood. Thus, a detailed understanding of the molecular-level interactions of Hg is crucial. The ongoing rapid development of hardware and methods has brought computational chemistry to a point that it can usefully inform environmental science. This is particularly true for Hg, which is difficult to handle experimentally due to its ultratrace concentrations in the environment and its toxicity. The current account provides a synopsis of the application of computational chemistry to filling several major knowledge gaps in environmental Hg chemistry that have not been adequately addressed experimentally. Environmental Hg chemistry requires defining the factors that determine the relative affinities of different ligands for Hg species, as they are critical for understanding its speciation, transformation and bioaccumulation in the environment. Formation constants and the nature of bonding have been determined computationally for environmentally relevant Hg(II) complexes such as chlorides, hydroxides, sulfides and selenides, in various physical phases. Quantum chemistry has been used to determine the driving forces behind the speciation of Hg with hydrochalcogenide and halide ligands. Of particular importance is the detailed characterization of solvation effects. Indeed, the aqueous phase reverses trends in affinities found computationally in the gas phase. Computation has also been used to investigate complexes of methylmercury with (seleno)amino acids, providing a molecular-level understanding of the toxicological antagonism between Hg and selenium (Se). Furthermore, evidence is emerging that ice surfaces play an important role in Hg transport and transformation in polar and alpine regions. Therefore, the diffusion of Hg and its ions through an idealized ice surface has been characterized. Microorganisms are major players in environmental mercury cycling. Some methylate inorganic Hg species, whereas others demethylate methylmercury. Quantum chemistry has been used to investigate catalytic mechanisms of enzymatic Hg methylation and demethylation. The complex interplay between the myriad chemical reactions and transport properties both in and outside microbial cells determines net biogeochemical cycling. Prospects for scaling up molecular work to obtain a mechanistic understanding of Hg cycling with comprehensive multiscale biogeochemical modeling are also discussed.

Relativistic Effects on a Metal-Metal Bond: Osmium Corrole Dimers

A series of metal-metal bonded osmium corrole dimers, {Os[T pXPC]}₂, were synthesized in reasonably good yields (35-46%) via the interaction of the corresponding free-base meso-tris( p-X-phenyl)corroles (H₃[T pXPC], X = CF₃, H, CH₃, and OCH₃), Os₃(CO)₁₂, and potassium carbonate in 1,2,4-trichlorobenzene under an inert atmosphere at 180 °C over several hours. The complexes are only the second class of Os corroles reported to date (the first being Os^VIN corroles) and also the second class of metal-metal bonded metallocorrole dimers (the other being Ru corrole dimers). Comparison of the X-ray structures, redox potentials, and optical spectra of analogous Ru and Os corrole dimers, along with scalar-relativistic DFT calculations, has provided an experimentally calibrated account of relativistic effects in these complexes. Three of the Os corrole dimers (X = CF₃, H, and OCH₃) were analyzed with single-crystal X-ray diffraction analysis, revealing inversion-related corrole rings with eclipsed Os-N bonds and Os-Os distances of ∼2.24 Å that are ∼0.06 Å longer than the Ru-Ru distances in the analogous Ru corrole dimers. Interestingly, a comparison of scalar-relativistic and nonrelativistic DFT calculations indicates that this difference in metal-metal bond distance does not, in fact, reflect a differential relativistic effect. For a given corrole ligand, the Ru and Os corrole dimers exhibit nearly identical oxidation potentials but dramatically different reduction potentials, with the Os values ∼0.5 V lower relative to Ru, suggesting that whereas oxidation occurs in a ligand-centered manner, reduction is substantially metal-centered, which indeed was confirmed by scalar-relativistic calculations. The calculations further indicate that approximately a third of the ∼0.5 V difference in reduction potentials can be ascribed to relativity. The somewhat muted value of this relativistic effect appears to be related to the finding that reduction of an Os corrole dimer is not exclusively metal-based but that a significant amount of spin density is delocalized over to the corrole ligand; in contrast, reduction of an Ru corrole dimer occurs exclusively on the Ru-Ru linkage. For isoelectronic complexes, the Ru and Os corrole dimers exhibit substantially different UV-vis spectra. A key difference is a strong near-UV feature of the Os series, which in energy terms is blue-shifted by ∼0.55 V relative to the analogous feature of the Ru series. TDDFT calculations suggest that this difference may be related to higher-energy Os(5d)-based LUMOs in the Os case relative to analogous MOs for Ru.

Microbial Mercury Methylation in Aquatic Environments: A Critical Review of Published Field and Laboratory Studies

Methylmercury (MeHg) is an environmental contaminant of concern because it biomagnifies in aquatic food webs and poses a health hazard to aquatic biota, piscivorous wildlife and humans. The dominant source of MeHg to freshwater systems is the methylation of inorganic Hg (IHg) by anaerobic microorganisms; and it is widely agreed that in situ rates of Hg methylation depend on two general factors: the activity of Hg methylators and their uptake of IHg. A large body of research has focused on the biogeochemical processes that regulate these two factors in nature; and studies conducted within the past ten years have made substantial progress in identifying the genetic basis for intracellular methylation and defining the processes that govern the cellular uptake of IHg. Current evidence indicates that all Hg methylating anaerobes possess the gene pair hgcAB that encodes proteins essential for Hg methylation. These genes are found in a large variety of anaerobes, including iron reducers and methanogens; but sulfate reduction is the metabolic process most often reported to show strong links to MeHg production. The uptake of Hg substrate prior to methylation may occur by passive or active transport, or by a combination of both. Competitive inhibition of Hg uptake by Zn speaks in favor of active transport and suggests that essential metal transporters are involved. Shortly after its formation, MeHg is typically released from cells, but the efflux mechanisms are unknown. Although methylation facilitates Hg depuration from the cell, evidence suggests that the hgcAB genes are not induced or favored by Hg contamination. Instead, high MeHg production can be linked to high Hg bioavailability as a result of the formation of Hg(SH)₂, HgS nanoparticles, and Hg-thiol complexes. It is also possible that sulfidic conditions require strong essential metal uptake systems that inadvertently bring Hg into the cytoplasm of Hg methylating microbes. In comparison with freshwaters, Hg methylation in open ocean waters appears less restricted to anoxic environments. It does seem to occur mainly in oxygen deficient zones (ODZs), and possibly within anaerobic microzones of settling organic matter, but MeHg (CH₃Hg⁺) and Me₂Hg ((CH₃)₂Hg) have been shown to form also in surface water samples from the euphotic zone. Future studies may disclose whether several different pathways lead to Hg methylation in marine waters and explain why Me₂Hg is a significant Hg species in oceans but seemingly not in most freshwaters.

Automated reaction path searches for spin-forbidden reactions

Many catalytic and biomolecular reactions containing transition metals involve changes in the electronic spin state. These processes are referred to as "spin-forbidden" reactions within nonrelativistic quantum mechanics framework. To understand detailed reaction mechanisms of spin-forbidden reactions, one must characterize reaction pathways on potential energy surfaces with different spin states and then identify crossing points. Here we propose a practical computational scheme, where only the lowest mixed-spin eigenstate obtained from the diagonalization of the spin-coupled Hamiltonian matrix is used in reaction path search calculations. We applied this method to the ^6,4 FeO⁺  + H₂ → ^6,4 Fe⁺  + H₂ O, ^6,4 FeO⁺  + CH₄ → ^6,4 Fe⁺  + CH₃ OH, and ⁷ Mn⁺  + OCS → ⁵ MnS⁺  + CO reactions, for which crossings between the different spin states are known to play essential roles in the overall reaction kinetics. © 2018 Wiley Periodicals, Inc.