Elucidating Heavy Atom Tunneling Kinetics in the Cope Rearrangement of Semibullvalene.

In this work, we characterize the temperature dependence of kinetic properties in heavy atom tunneling reactions by means of molecular dynamics simulations, including nuclear quantum effects (NQEs) via Path Integral theory. To this end, we consider the prototypical Cope rearrangement of semibullvalene. The reaction was studied in the 25-300K temperature range observing that the inclusion of NQEs modifies the temperature behavior of both free energy barriers and dynamical recrossing factors with respect to classical dynamics. Notably, while in classical simulations the activation free energy shows a very weak temperature dependence, it becomes strongly dependent on temperature when NQEs are included. This temperature behavior shows a transition from a regime where the quantum effects are limited and can mainly be traced back to zero point energy, to a low temperature regime where tunneling plays a dominant role. In this regime, the free energy curve tunnels below the potential energy barrier along the reaction coordinate,  allowing much faster reaction rates. Finally, the temperature dependence of the rate constants obtained from molecular dynamics simulations was compared with available experimental data and with semi-classical transition state theory calculations, showing comparable behaviors and similar transition temperatures from  thermal to (deep) tunneling regime.

C5 Pentacle Structures: A Localization-Delocalization Matrices Approach.

This article explores the possible presence of a pentacle valence bond structure in C 5 ${_5 }$ cyclic molecules. At this end, we have used quantum chemistry tools to elucidate the possible arrangement and the nature of chemical bonds within linear, cyclic, and three-dimensional structures only formed by five carbon atoms. While the linear structure is clearly the most stable one, local minima were obtained for both bi- and three-dimensional structures. Using the localization-delocalization matrices approach, we characterize both the minimum linear structure and the cyclic ones. Interestingly, the linear structure is a combination of ionic and covalent bonds, albeit the four distances are almost identical, when using Density Functional Theory. For cyclic C 5 ${_5 }$ , the pentacle bonding arrangement emerges as a significant Lewis structure, indicative of an unusual formal configuration characterized by five intersecting C-C bonds. Our calculations show that this pentacle arrangement in cyclic C 5 ${_5 }$ scheme is also present in the more known cyclo-pentadienyl molecule.

Piezoelectric and microfluidic tuning of an infrared cavity for vibrational polariton studies.

We developed a microfluidic system for vibrational polariton studies, which consists of two microfluidic chips: one for solution mixing and another for tuning an infrared cavity made of a pair of gold mirrors and a PDMS (polydimethylsiloxane) spacer. We show that the cavity of the system can be accurately tuned with either piezoelectric actuators or microflow-induced pressure to result in resonant coupling between a cavity mode and a variational mode of the solution molecules. Acrylonitrile solutions were chosen to prove the concept of vabriational strong coupling (VSC) of a CN stretching mode with light inside the cavity. We also show that the Rabi splitting energy is linearly proportional to the square root of molecular concentration, thereby proving the relevance and reliability of the system for VSC studies.

Quasi-classical simulations of resonance Raman spectra based on path integral linearization.

Based on a linearization approximation coupled with path integral formalism, we propose a method derived from the propagation of quasi-classical trajectories to simulate resonance Raman spectra. This method is based on ground state sampling followed by an ensemble of trajectories on the mean surface between the ground and excited states. The method was tested on three models and compared to a quantum mechanics solution based on a sum-over-states approach: harmonic and anharmonic oscillators and the HOCl molecule (hypochlorous acid). The method proposed is able to correctly characterize resonance Raman scattering and enhancement, including the description of overtones and combination bands. The absorption spectrum is obtained at the same time, and the vibrational fine structure can be reproduced for long excited state relaxation times. The method can also be applied to dissociating excited states (as is the case for HOCl).

Environmental and Nuclear Quantum Effects on Double Proton Transfer in the Guanine-Cytosine Base Pair.

In the present letter, we investigate the double proton transfer (DPT) tautomerization process in guanine-cytosine (GC) DNA base pairs. In particular, we study the influence of the biological environment on the mechanism, the kinetics and thermodynamics of such DPT. To this end, we present a molecular dynamics (MD) study in the tight-binding density functional theory framework, and compare the reactivity of the isolated GC dimer with that of the same dimer embedded in a small DNA structure. The impact of nuclear quantum effects (NQEs) is also evaluated using Path Integral based MD. Results show that in the isolated dimer, the DPT occurs via a concerted mechanism, while in the model biological environment, it turns into a stepwise process going through an intermediate structure. One of the water molecules in the vicinity of the proton transfer sites plays an important role as it changes H-bond pattern during the DPT reaction. The inclusion of NQEs has the effect of speeding up the tautomeric-to-canonical reaction, reflecting the destabilization of both the tautomeric and intermediate forms.

Light-induced infrared difference spectroscopy on three different forms of orange carotenoid protein: focus on carotenoid vibrations.

Orange carotenoid protein (OCP) is a photoactive carotenoprotein involved in photoprotection of cyanobacteria, which uses a keto-catorenoid as a chromophore. When it absorbs blue-green light, it converts from an inactive OCPO orange form to an activated OCPR red form, the latter being able to bind the light-harvesting complexes facilitating thermal dissipation of the excess of absorbed light energy. Several research groups have focused their attention on the photoactivation mechanism, characterized by several steps, involving both carotenoid photophysics and protein conformational changes. Among the used techniques, time-resolved IR spectroscopy have the advantage of providing simultaneously information on both the chromophore and the protein, giving thereby the possibility to explore links between carotenoid dynamics and protein dynamics, leading to a better understanding of the mechanism. However, an appropriate interpretation of data requires previous assignment of marker IR bands, for both the carotenoid and the protein. To date, some assignments have concerned specific α-helices of the OCP backbone, but no specific marker band for the carotenoid was identified on solid ground. This paper provides evidence for the assignment of putative marker bands for three carotenoids bound in three different OCPs: 3'-hydroxyechineone (3'-hECN), echinenone (ECN), canthaxanthin (CAN). Light-induced FTIR difference spectra were recorded in H2O and D2O and compared with spectra of isolated carotenoids. The use of DFT calculations allowed to propose a description for the vibrations responsible of several IR bands. Interestingly, most bands are located at the same wavenumber for the three kinds of OCPs suggesting that the conformation of the three carotenoids is the same in the red and in the orange form. These results are discussed in the framework of recent time-resolved IR studies on OCP.

Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations.

In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.

Challenges in the Development of Intravenous Neurokinin-1 Receptor Antagonists: Results of a Safety and Pharmacokinetics Dose-Finding, Phase 1 Study of Intravenous Fosnetupitant.

Oral NEPA is the fixed-combination antiemetic comprising netupitant (neurokinin-1 receptor antagonist [NK1 RA]) and palonosetron (5-hydroxytryptamine-3 receptor antagonist [5-HT3 RA]). Intravenous (IV) NEPA, containing fosnetupitant, a water-soluble N-phosphoryloxymethyl prodrug of netupitant, has been developed. Fosnetupitant does not require excipients or solubility enhancers often used to increase IV NK1 RA water solubility, preventing the occurrence of hypersensitivity and infusion-site reactions associated with these products. In this phase 1 study, subjects received a 30-minute placebo or fosnetupitant (17.6-353 mg) infusion and an oral NEPA or placebo capsule, with 2-sequence crossover treatment for fosnetupitant 118- to 353-mg dose cohorts. IV fosnetupitant safety and pharmacokinetics were evaluated, and its equivalence to an oral netupitant 300-mg dose was defined. Overall, 158 healthy volunteers were enrolled. All adverse events (AEs) were mild or moderate in intensity. Doppler-identified infusion-site asymptomatic thrombosis occurred in 5.4% (fosnetupitant) and 1.2% (oral NEPA) of subjects. The frequency or number of treatment-related AEs did not increase with ascending fosnetupitant doses. The most common treatment-related AEs were headache (fosnetupitant, 8.1%; oral NEPA, 12.7%) and constipation (fosnetupitant, 1.4%; oral NEPA, 7.5%). A fosnetupitant 235-mg dose was equivalent, in terms of netupitant exposure, to 300-mg oral netupitant. The safety profile of a single fosnetupitant 235-mg infusion was similar to that of single-dose oral NEPA.

Efficient and Accurate Description of Diels-Alder Reactions Using Density Functional Theory.

Modeling chemical reactions using Quantum Chemistry is a widely used predictive strategy capable to complement experiments in order to understand the intrinsic mechanisms guiding the chemicals towards the most favorable reaction products. However, at this purpose, it is mandatory to use reliable and computationally tractable theoretical methods. In this work, we focus on six Diels-Alder reactions of increasing complexity and perform an extensive benchmark of middle- to low-cost computational approaches to predict the characteristic reactions energy barriers. We found that Density Functional Theory, using the ωB97XD, LC-ωPBE, CAM-B3LYP, M11 and MN12SX functionals, with empirical dispersion corrections coupled to an affordable 6-31G basis set, provides quality results for this class of reactions, at a small computational effort. Such efficient and reliable simulation protocol opens perspectives for hybrid QM/MM molecular dynamics simulations of Diels-Alder reactions including explicit solvation.

Thermochemical Studies of Ni(II) and Zn(II) Ternary Complexes using Ion Mobility-Mass Spectrometry.

This article describes an experimental protocol using electrospray-ion mobility-mass spectrometry (ES-IM-MS) and energy-resolved threshold collision-induced dissociation (TCID) to measure the thermochemistry of the dissociation of negatively-charged [amb+M(II)+NTA]- ternary complexes into two product channels: [amb+M(II)] + NTA or [NTA+M(II)]- + amb, where M = Zn or Ni and NTA is nitrilotriacetic acid. The complexes contain one of the alternative metal binding (amb) heptapeptides with the primary structures acetyl-His1-Cys2-Gly3-Pro4-Tyr5-His6-Cys7 or acetyl-Asp1-Cys2-Gly3-Pro4-Tyr5-His6-Cys7, where the amino acids' Aa1,2,6,7 positions are the potential metal-binding sites. Geometry-optimized stationary states of the ternary complexes and their products were selected from quantum chemistry calculations (presently the PM6 semi-empirical Hamiltonian) by comparing their electronic energies and their collision cross-sections (CCS) to those measured by ES-IM-MS. From the PM6 frequency calculations, the molecular parameters of the ternary complex and its products model the energy-dependent intensities of the two product channels using a competitive TCID method to determine the threshold energies of the reactions that relate to the 0 K enthalpies of dissociation (ΔH0). Statistical mechanics thermal and entropy corrections using the PM6 rotational and vibrational frequencies provide the 298 K enthalpies of dissociation (ΔH298). These methods describe an EI-IM-MS routine that can determine thermochemistry and equilibrium constants for a range of ternary metal ion complexes.

Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories.

In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.

Implementing Historical Controls in Oncology Trials.

Drug development in oncology has broadened from mainly considering randomized clinical trials to also including single-arm trials tailored for very specific subtypes of cancer. They often use historical controls, and this article discusses benefits and risks of this paradigm and provide various regulatory and statistical considerations. While leveraging the information brought by historical controls could potentially shorten development time and reduce the number of patients enrolled, a careful selection of the past studies, a prespecified statistical analysis accounting for the heterogeneity between studies, and early engagement with regulators will be key to success. Although both the European Medicines Agency and the U.S. Food and Drug Administration have already approved medicines based on nonrandomized experiments, the evidentiary package can be perceived as less comprehensive than randomized experiments. Use of historical controls, therefore, is better suited for cases of high unmet clinical need, where the disease course is well characterized and the primary endpoint is objective. IMPLICATIONS FOR PRACTICE: Incorporating historical data in single-arm oncology trials has the potential to accelerate drug development and to reduce the number of patients enrolled, compared with standard randomized controlled clinical trials. Given the lack of blinding and randomization, such an approach is better suited for cases of high unmet clinical need and/or difficult experimental situations, in which the trajectory of the disease is well characterized and the endpoint can be measured objectively. Careful pre-specification and selection of the historical data, matching of the patient characteristics with the concurrent trial data, and innovative statistical methodologies accounting for between-study variation will be needed. Early engagement with regulators (e.g., via Scientific Advice) is highly recommended.

Unimolecular Fragmentation Properties of Thermometer Ions from Chemical Dynamics Simulations.

Thermometer ions are widely used to calibrate the internal energy of the ions produced by electrospray ionization in mass spectrometry. Typically, benzylpyridinium ions with different substituents are used. More recently, benzhydrylpyridinium ions were proposed for their lower bond dissociation energies. Direct dynamics simulations using M06-2X/6-31G(d), DFTB, and PM6-D3 are performed to characterize the activation energies of two representative systems: para-methylbenzylpyridinium ion (p-Me-BnPy+) and methyl,methylbenzhydrylpyridinium ion (Me,Me-BhPy+). Simulation results are used to calculate rate constants for the two systems. These rate constants and their uncertainties are used to find the Arrhenius activation energies and RRK fitted threshold energies which give reasonable agreement with calculated bond dissociation energies at the same level of theory. There is only one fragmentation mechanism observed for both systems, which involves C-N bond dissociation via a loose transition state, to generate either benzylium or benzhydrylium ion and a neutral pyridine molecule. For p-Me-BnPy+ using DFTB and PM6-D3 the formation of tropylium ion, from rearrangement of benzylium ion, was observed but only at higher excitation energies and for longer simulation times. These observations suggest that there is no competition between reaction pathways that could affect the reliability of internal energy calibrations. Finally, we suggest using DFTB with a modified-Arrhenius model in future studies.

Exploring reactivity and product formation in N(4S) collisions with pristine and defected graphene with direct dynamics simulations.

Atomic nitrogen is formed in the high-temperature shock layer of hypersonic vehicles and contributes to the ablation of their thermal protection systems (TPSs). To gain atomic-level understanding of the ablation of carbon-based TPS, collisions of hyperthermal atomic nitrogen on representative carbon surfaces have recently be investigated using molecular beams. In this work, we report direct dynamics simulations of atomic-nitrogen [N(4S)] collisions with pristine, defected, and oxidized graphene. Apart from non-reactive scattering of nitrogen atoms, various forms of nitridation of graphene were observed in our simulations. Furthermore, a number of gaseous molecules, including the experimentally observed CN molecule, have been found to desorb as a result of N-atom bombardment. These results provide a foundation for understanding the molecular beam experiment and for modeling the ablation of carbon-based TPSs and for future improvement of their properties.

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.

Chemical dynamics simulations are performed to study the collision induced gas phase unimolecular fragmentation of a model peptide with the sequence acetyl-His1-Cys2-Gly3-Pro4-Tyr5-His6-Cys7 (analogue methanobactin peptide-5, amb5) and in particular to explore the role of zinc binding in reactivity. Fragmentation pathways, their mechanisms, and collision energy transfer are discussed. The probability distributions of the pathways are compared with the results of the experimental IM-MS, MS/MS spectrum and previous thermal simulations. Collisional activation gives both statistical and non-statistical fragmentation pathways with non-statistical shattering mechanisms accounting for a relevant percentage of reactive trajectories, becoming dominant at higher energies. The tetra-coordination of zinc changes qualitative and quantitative fragmentation, in particular the shattering. The collision energy threshold for the shattering mechanism was found to be 118.9 kcal mol-1 which is substantially higher than the statistical Arrhenius activation barrier of 35.8 kcal mol-1 identified previously during thermal simulations. This difference can be attributed to the tetra-coordinated zinc complex that hinders the availability of the sidechains to undergo direct collision with the Ar projectile.

Solvation effects drive the selectivity in Diels-Alder reaction under hyperbaric conditions.

High pressure effects on the Diels-Alder reaction in condensed phase are investigated by means of theoretical methods, employing advanced multiscale modeling approaches based on physically grounded models. The simulations reveal how the increase of pressure from 1 to 10 000 atm (10 katm) does not affect the stability of the reaction products, modifying the kinetics of the process by lowering considerably the transition state energy. The reaction profile at high pressure remarkably differs from that at 1 atm, showing a submerged TS and a pre-TS structure lower in energy. The different solvation between endo and exo pre-TS is revealed as the driving force pushing the reaction toward a much higher preference for the endo product at high pressure.

Infrared-Assisted Synthesis of Prebiotic Glycine.

A novel approach has been developed to synthesize complex organic molecules (COMs) relevant to prebiotic chemistry, using infrared (IR) radiation to trigger the reaction. An original laboratory reactor working at low gas density and using IR irradiation was developed. In this way, glycine, the simplest brick of life, has been synthesized by assisting ion-molecule reaction with IR laser light. The ion-molecule complex constituted by acetic acid and hydroxylamine was formed in a mass spectrometer reactor and then irradiated with IR photons. As photoproducts, we obtained both glycine structures and some of its isomers. Anharmonic vibrational frequency calculations and fragmentation dynamics simulations allow for a better interpretation of the experimental data. This novel approach can be now extended to study other new synthetic pathways responsible for the formation of further COMs also with potential prebiotic relevance.

Clustered allocation as a way of understanding historical controls: Components of variation and regulatory considerations.

There has been increasing interest in recent years in the possibility of increasing the efficiency of clinical trials by using historical controls. There has been a general recognition that in replacing concurrent by historical controls, the potential for bias is serious and requires some down-weighting to the apparent amount of historical information available. However, such approaches have generally assumed that what is required is some modification to the standard inferential model offered by the parallel group trial. In our opinion, the correct starting point that requires modification is a trial in which treatments are allocated to clusters. This immediately shows that the amount of information available is governed not just by the number of historical patients but also by the number of centres and of historical studies. Furthermore, once one accepts that external patients may be used as controls, this raises the issue as to which patients should be used. Thus, abandoning concurrent control has implications for many aspects of design and analysis of trials, including (a) identification, pre-specification and agreement on a suitable historical dataset; (b) an agreed, enforceable and checkable plan for recruiting the experimental arm; (c) a finalised analysis plan prior to beginning the trial and (d) use of a hierarchical model with sufficient complexity. We discuss these issues and suggest approaches to design and analysis making extensive reference to the partially randomised Therapeutic Arthritis Research and Gastrointestinal Event Trial study. We also compare some Bayesian and frequentist approaches and provide some important regulatory considerations. We conclude that effective use of historical data will require considerable circumspection and discipline.

Fragmentation Spectra Prediction and DNA Adducts Structural Determination.

In this work, chemical dynamics simulations were optimized and used to predict fragmentation mass spectra for DNA adduct structural determination. O6-methylguanine (O6-Me-G) was used as a simple model adduct to calculate theoretical spectra for comparison with measured high-resolution fragmentation data. An automatic protocol was established to consider the different tautomers accessible at a given energy and obtain final theoretical spectra by insertion of an initial tautomer. In the work reported here, the most stable tautomer was chosen as the initial structure, but in general, any structure could be considered. Allowing for the formation of the various possible tautomers during simulation calculations was found to be important to getting a more complete fragmentation spectrum. The calculated theoretical results reproduce the experimental peaks such that it was possible to determine reaction pathways and product structures. The calculated tautomerization network was crucial to correctly identifying all the observed ion peaks, showing that a mobile proton model holds not only for peptide fragmentation but also for nucleobases. Finally, first principles results were compared to simple machine learning fragmentation models.