Automatic Determination of the Non-Covalent Stable Conformations of the NO2-Pyrene Cluster in Full Dimensionality (81D) Using the vdW-TSSCDS Approach

We present the detailed topographical characterisation (stationary points and minimum energy paths connecting them) of the full dimensional (81D) intermolecular potential energy surface associated with the non-covalent interactions between the NO2 radical and the pyrene (C16H10) molecule. The whole procedure is (quasi) fully automated. We have used our recent algorithm vdW-TSSCDS as implemented on the freely-available AutoMekin software package. To this end, a series of inexpensive classical trajectories using forces from a low-level (semi-empirical) theory are used to sample the configuration space of the system in the search for candidates to first order saddle points. These guess structures are determined by means of a graph-theory based algorithm using the concept of adjacency matrix. Low-level optimizations are followed by re-optimizations at a final high-level of theory (DFT and CCSD(T)-F12 in our case.). The resulting set of stationary points and paths connecting them constitutes the so-called reaction network. In the case of NO2-pyrene, this network exhibits four major basins which can be characterized by their point-group symmetry. A central one, of global C2 symmetry, comprises the global minimum (as well as all other permutationally related conformers) together with the corresponding C2v saddle points connecting them. This central basin is connected to three others of lower C1 symmetry. The latter can be distinguished by the projection of the position of the NO2 nitrogen atom on the pyrene plane in combination with the relative orientation of the oxygen pair pointing either inwards, outwards, upwards or downwards.

Mechanism of Prebiotic Uracil Synthesis from Urea and HC3O+ in Space

The potential energy surface for the formation of protonated uracil (UH⁺) from urea and HC₃O⁺ was explored by performing quantum chemical complete basis set-QB3 calculations. A barrierless pathway was found for the formation of UH⁺, which was estimated to occur in the interstellar medium (ISM) much faster than the timescale of chemical revolution of typical dense interstellar clouds. Investigation of further reactions of UH⁺ formed through the obtained pathway led to the conclusion that uracil could be produced on icy grain surfaces but not in the gas phase of the ISM.

New Tools for Taming Complex Reaction Networks: The Unimolecular Decomposition of Indole Revisited

The level of detail attained in the computational description of reaction mechanisms can be vastly improved through tools for automated chemical space exploration, particularly for systems of small to medium size. Under this approach, the unimolecular decomposition landscape for indole was explored through the automated reaction mechanism discovery program AutoMeKin. Nevertheless, the sheer complexity of the obtained mechanisms might be a hindrance regarding their chemical interpretation. In this spirit, the new Python library amk-tools has been designed to read and manipulate complex reaction networks, greatly simplifying their overall analysis. The package provides interactive dashboards featuring visualizations of the network, the three-dimensional (3D) molecular structures and vibrational normal modes of all chemical species, and the corresponding energy profiles for selected pathways. The combination of the joined mechanism generation and postprocessing workflow with the rich chemistry of indole decomposition enabled us to find new details of the reaction (obtained at the CCSD(T)/aug-cc-pVTZ//M06-2X/MG3S level of theory) that were not reported before: (i) 16 pathways leading to the formation of HCN and NH₃ (via amino radical); (ii) a barrierless reaction between methylene radical and phenyl isocyanide, which might be an operative mechanism under the conditions of the interstellar medium; and (iii) reaction channels leading to both hydrogen cyanide and hydrogen isocyanide, of potential astrochemical interest as the computed HNC/HCN ratios greatly exceed the calculated equilibrium value at very low temperatures. The reported reaction networks can be very valuable to supplement databases of kinetic data, which is of remarkable interest for pyrolysis and astrochemical studies.

Electron and ion spectroscopy of the cyclo-alanine-alanine dipeptide

The VUV photoionisation and photofragmentation of cyclo-alanine-alanine (cAA) has been studiedin a joint experimental and theoretical work. The photoelectron spectrum and the photoelectron-photoion coincidence (PEPICO) measurements, which enable a control...

Exploring the Effects of Methylation on the CID of Protonated Lysine: A Combined Experimental and Computational Approach

We report the results of experiments, simulations, and DFT calculations that focus on describing the reaction dynamics observed within the collision-induced dissociation of l-lysine-H⁺ and its side-chain methylated analogues, N_ε-methyl-l-lysine-H⁺ (Me₁-lysine-H⁺), N_ε,N_ε-dimethyl-l-lysine-H⁺ (Me₂-lysine-H⁺), and N_ε,N_ε,N_ε-trimethyl-l-lysine-H⁺ (Me₃-lysine-H⁺). The major pathways observed in the experimental measurements were m/z 130 and 84, with the former dominant at low collision energies and the latter at intermediate to high collision energies. The m/z 130 peak corresponds to loss of N(CH₃)_nH_3-n, while m/z 84 has the additional loss of H₂CO₂ likely in the form of H₂O + CO. Within the time frame of the direct dynamics simulations, m/z 130 and 101 were the most populous peaks, with the latter identified as an intermediate to m/z 84. The simulations allowed for the determination of several reaction pathways that result in these products. A graph theory analysis enabled the elucidation of the significant structures that compose each peak. Methylation results in the preferential loss of the side-chain amide group and a reduction of cyclic structures within the m/z 84 peak population in simulations.

AutoMeKin2021: An open-source program for automated reaction discovery

AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin.

Current Challenges and Recent Developments in Mass Spectrometry-Based Metabolomics

High-resolution mass spectrometry (MS) has advanced the study of metabolism in living systems by allowing many metabolites to be measured in a single experiment. Although improvements in mass detector sensitivity have facilitated the detection of greater numbers of analytes, compound identification strategies, feature reduction software, and data sharing have not kept up with the influx of MS data. Here, we discuss the ongoing challenges with MS-based metabolomics, including de novo metabolite identification from mass spectra, differentiation of metabolites from environmental contamination, chromatographic separation of isomers, and incomplete MS databases. Because of their popularity and sensitive detection of small molecules, this review focuses on the challenges of liquid chromatography-mass spectrometry-based methods. We then highlight important instrumentational, experimental, and computational tools that have been created to address these challenges and how they have enabled the advancement of metabolomics research.

A map of mass spectrometry-based in silico fragmentation prediction and compound identification in metabolomics

Metabolomics, the comprehensive study of the metabolome, and lipidomics-the large-scale study of pathways and networks of cellular lipids-are major driving forces in enabling personalized medicine. Complicated and error-prone data analysis still remains a bottleneck, however, especially for identifying novel metabolites. Comparing experimental mass spectra to curated databases containing reference spectra has been the gold standard for identification of compounds, but constructing such databases is a costly and time-demanding task. Many software applications try to circumvent this process by utilizing cutting-edge advances in computational methods-including quantum chemistry and machine learning-and simulate mass spectra by performing theoretical, so called in silico fragmentations of compounds. Other solutions concentrate directly on experimental spectra and try to identify structural properties by investigating reoccurring patterns and the relationships between them. The considerable progress made in the field allows recent approaches to provide valuable clues to expedite annotation of experimental mass spectra. This review sheds light on individual strengths and weaknesses of these tools, and attempts to evaluate them-especially in view of lipidomics, when considering complex mixtures found in biological samples as well as mass spectrometer inter-instrument variability.

Excited State Deactivation Mechanism in Protonated Uracil: New Insights from Theoretical Studies

We have conducted here a theoretical exploration, discussing the distinct excited state lifetimes reported experimentally for the two lowest lying protonated isomers of uracil. In this regard, the first-principal computational levels as well as the nonadiabatic surface hopping dynamics have been employed. It has been revealed that relaxation of the ¹ππ* state of enol-enol form (EE⁺) to the ground is barrier-free via out-of-plane coordinates, resulting in an ultrashort S₁ lifetime of this species. For the second most stable isomer (EK⁺), however, a significant barrier predicted in the CASPT2 S₁ potential energy profile along the twisting coordinate has been proposed to explain the relevant long lifetime reported experimentally.

Insights into Cisplatin Binding to Uracil and Thiouracils from IRMPD Spectroscopy and Tandem Mass Spectrometry

The monofunctional primary complexes cis-[PtCl(NH₃)₂(L)]⁺, formed by the reaction of cisplatin, a major chemotherapeutic agent, with four nucleobases L, i.e., uracil (U), 2-thiouracil (2SU), 4-thiouracil (4SU), and 2,4-dithiouracil (24dSU), have been studied by a combination of infrared multiple photon dissociation (IRMPD) action spectroscopy in both the fingerprint (900-1900 cm^-1) and the N-H/O-H stretching (3000-3800 cm^-1) ranges, energy-resolved collision-induced dissociation (CID) mass spectrometry, and density functional calculations at the B3LYP/LACVP/6-311G** level. On the basis of the comparison across the experimental features and the linear IR spectra of conceivable structures, the cisplatin residue is found to promote a monodentate interaction preferentially with the O4(S4) atoms of the canonical forms of U, 4SU, and 24dSU and to the S2 atom of 2SU, yielding the most stable structures. Additional absorptions reveal the presence of minor, alternative tautomers in the sampled ion populations of 2SU and 24dSU, underlying the ability of cisplatin to increase the prospect of (therapeutically beneficial) nucleic acid strand disorder. Implication of these evidence may provide insights into drug mechanism and design.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.

Collisional dynamics simulations revealing fragmentation properties of Zn(ii)-bound poly-peptide

Collisional simulations show how peptide fragmentation is modified by the presence of Zn(ii).

A Trajectory-Based Method to Explore Reaction Mechanisms

The tsscds method, recently developed in our group, discovers chemical reaction mechanisms with minimal human intervention. It employs accelerated molecular dynamics, spectral graph theory, statistical rate theory and stochastic simulations to uncover chemical reaction paths and to solve the kinetics at the experimental conditions. In the present review, its application to solve mechanistic/kinetics problems in different research areas will be presented. Examples will be given of reactions involved in photodissociation dynamics, mass spectrometry, combustion chemistry and organometallic catalysis. Some planned improvements will also be described.

Chemical dynamics simulations of CID of peptide ions: comparisons between TIK(H+)2 and TLK(H+)2 fragmentation dynamics, and with thermal simulations

Gas phase unimolecular fragmentation of the two model doubly protonated tripeptides threonine-isoleucine-lysine (TIK) and threonine-leucine-lysine (TLK) is studied using chemical dynamics simulations. Attention is focused on different aspects of collision induced dissociation (CID): fragmentation pathways, energy transfer, theoretical mass spectra, fragmentation mechanisms, and the possibility of distinguishing isoleucine (I) and leucine (L). Furthermore, discussion is given regarding the differences between single collision CID activation, which results from a localized impact between the ions and a colliding molecule N₂, and previous thermal activation simulation results; Z. Homayoon, S. Pratihar, E. Dratz, R. Snider, R. Spezia, G. L. Barnes, V. Macaluso, A. Martin-Somer and W. L. Hase, J. Phys. Chem. A, 2016, 120, 8211-8227. Upon thermal activation unimolecular fragmentation is statistical and in accord with RRKM unimolecular rate theory. Simulations show that in collisional activation some non-statistical fragmentation occurs, including shattering, which is not present when the ions dissociate statistically. Products formed by non-statistical shattering mechanisms may be related to characteristic mass spectrometry peaks which distinguish the two isomers I and L.

Threshold for shattering fragmentation in collision-induced dissociation of the doubly protonated tripeptide TIK(H+)2

Determination of shattering threshold for unimolecular dissociation of a model tripeptide.

An automated method to find reaction mechanisms and solve the kinetics in organometallic catalysis

A novel computational method is proposed in this work for use in discovering reaction mechanisms and solving the kinetics of transition metal-catalyzed reactions. The method does not rely on either chemical intuition or assumed a priori mechanisms, and it works in a fully automated fashion. Its core is a procedure, recently developed by one of the authors, that combines accelerated direct dynamics with an efficient geometry-based post-processing algorithm to find transition states (Martinez-Nunez, E., J. Comput. Chem.2015, 36, 222-234). In the present work, several auxiliary tools have been added to deal with the specific features of transition metal catalytic reactions. As a test case, we chose the cobalt-catalyzed hydroformylation of ethylene because of its well-established mechanism, and the fact that it has already been used in previous automated computational studies. Besides the generally accepted mechanism of Heck and Breslow, several side reactions, such as hydrogenation of the alkene, emerged from our calculations. Additionally, the calculated rate law for the hydroformylation reaction agrees reasonably well with those obtained in previous experimental and theoretical studies.

Gas-phase interactions of organotin compounds with cysteine

The gas-phase interactions of cysteine with di-organotin and tri-organotin compounds have been studied by mass spectrometry experiments and quantum calculations. Positive-ion electrospray spectra show that the interaction of di- and tri-organotins with cysteine results in the formation of [(R)₂ Sn(Cys-H)]⁺ and [(R)₃ Sn(Cys)]⁺ ions, respectively. MS/MS spectra of [(R)₂ Sn(Cys-H)]⁺ complexes are characterized by numerous fragmentation processes, notably associated with elimination of NH₃ and (C,H₂ ,O₂ ). Several dissociation routes are characteristic of each given organic species. Upon collision, both the [(R)₃ Sn(Gly)]⁺ and [(R)₃ Sn(Cys)]⁺ complexes are associated with elimination of the intact amino acid, leading to the formation of [(R)₃ Sn]⁺ cation. But for the latter complex, two additional fragmentation processes are observed, associated with the elimination of NH₃ and C₃ H₄ O₂ S. Calculations indicate that the interaction between organotins and cysteine is predominantly electrostatic but also exhibits a considerable covalent character, which is slightly more pronounced in tri-organotin complexes. A preferred bidentate interaction of the type -η² -S-NH₂ , with sulfur and the amino group, is observed. As for the [(R)₃ Sn(Cys)]⁺ complexes, their stability is due to the combination of the hydrogen bond taking place between the amino group and the sulfur lone pair and the interaction between the carboxylic oxygen atom and the metal. Copyright © 2016 John Wiley & Sons, Ltd.