The inclusion of Quantum Mechanical (QM) effects such as zero point energy (ZPE) and tunneling in simulations of chemical reactions, especially in the case of light atom transfer, is an important problem in computational chemistry. In this respect, the hydrogen exchange reaction and its isotopic variants constitute an excellent benchmark for the assessment of approximate QM methods. In particular, the recently developed ring polymer molecular dynamics (RPMD) technique has been demonstrated to give very good results for bimolecular chemical reactions in the gas phase. In this work, we have performed a detailed RPMD study of the H + H₂ reaction and its isotopologues Mu + H₂, D + H₂ and Heμ + H₂, at temperatures ranging from 200 to 1000 K. Thermal rate coefficients and kinetic isotope effects have been computed and compared with exact QM calculations as well as with quasiclassical trajectories and experiment. The agreement with the QM results is good for the heaviest isotopologues, with errors ranging from 15% to 45%, and excellent for Mu + H₂, with errors below 15%. We have seen that RPMD is able to capture the ZPE effect very accurately, a desirable feature of any method based on molecular dynamics. We have also verified Richardson and Althorpe's prediction [J. O. Richardson and S. C. Althorpe, J. Chem. Phys., 2009, 131, 214106] that RPMD will overestimate thermal rates for asymmetric reactions and underestimate them for symmetric reactions in the deep tunneling regime. The ZPE effect along the reaction coordinate must be taken into account when assigning the reaction symmetry in the multidimensional case.

Cold collisions of light molecules are often dominated by a single partial wave resonance. For the rotational quenching of HD ($v=1$, $j=2$) by collisions with ground state para-${\mathrm{H}}_2$, the process is dominated by a single $L=2$ partial wave resonance centered around 0.1 K. Here, we show that this resonance can be switched on or off simply by appropriate alignment of the HD rotational angular momentum relative to the initial velocity vector, thereby enabling complete control of the collision outcome.

RAF kinases are key players in the MAPK signaling pathway and are important targets for personalized cancer therapy. RAF dimerization is part of the physiological activation mechanism, together with phosphorylation, and is known to convey resistance to RAF inhibitors. Herein, molecular dynamics simulations are used to show that phosphorylation of a key N‐terminal acidic (NtA) motif facilitates RAF dimerization by introducing several interprotomer salt bridges between the αC‐helix and charged residues upstream of the NtA motif. Additionally, we show that the R‐spine of RAF interacts with a conserved Trp residue in the vicinity of the NtA motif, connecting the active sites of two protomers and thereby modulating the cooperative interactions in the RAF dimer. Our findings provide a first structure‐based mechanism for the auto‐transactivation of RAF and could be generally applicable to other kinases, opening new pathways for overcoming dimerization‐related drug resistance.

Protein phosphorylation is one of the most common post-translational modifications in cell regulatory mechanisms. Dimerization plays also a crucial role in the kinase activity of many kinases, including RAF, CDK2 (cyclin-dependent kinase 2) and EGFR (epidermal growth factor receptor), with heterodimers often being the most active forms. However, the structural and mechanistic details of how phosphorylation affects the activity of homo- and hetero-dimers are largely unknown. Experimentally, synthesizing protein samples with fully specified and homogeneous phosphorylation states remains a challenge for structural biology and biochemical studies. Typically, multiple changes in phosphorylation lead to activation of the same protein, which makes structural determination methods particularly difficult. It is also not well understood how the occurrence of phosphorylation and dimerization processes synergize to affect kinase activities. In the present article, we review available structural data and discuss how MD simulations can be used to model conformational transitions of RAF kinase dimers, in both their phosphorylated and unphosphorylated forms.

A detailed study of the proton exchange reaction H⁺ + D₂(v = 0, j = 0) → HD + D⁺ on its ground 1¹A′ potential energy surface has been carried out using ‘exact’ close-coupled quantum mechanical wavepacket (WP-EQM), quasi-classical trajectory (QCT), and statistical quasi-classical trajectory (SQCT) calculations for a range of collision energies starting from the reaction threshold to 1.3 eV. The WP-EQM calculations include all total angular momenta up to J_max = 50, and therefore the various dynamical observables are converged up to 0.6 eV. It has been found that it is necessary to include all Coriolis couplings to obtain reliable converged results. Reaction probabilities obtained using the different methods are thoroughly compared as a function of the total energy for a series of J values. Comparisons are also made of total reaction cross sections as function of the collision energy, and rate constants. In addition, opacity functions, integral cross sections (ICS) and differential cross sections (DCS) are presented at 102 meV, 201.3 meV and 524.6 meV collision energy. The agreement between the three sets of results is only qualitative. The QCT calculations fail to describe the overall reactivity and most of the dynamical observables correctly. At low collision energies, the QCT method is plagued by the lack of conservation of zero point energy, whilst at higher collision energies and/or total angular momenta, the appearance of an effective repulsive potential associated with the centrifugal motion “over” the well causes a substantial decrease of the reactivity. In turn, the statistical models overestimate the reactivity over the whole range of collision energies as compared with the WP-EQM method. Specifically, at sufficiently high collision energies the reaction cannot be deemed to be statistical and important dynamical effects seem to be present. In general the WP-EQM results lie in between those obtained using the QCT and SQCT methods. One of the main, unexpected, conclusions of this work is that an accurate description of the reaction and of its various dynamical features requires a computationally expensive, accurate quantum mechanical treatment.