More than 100 years since its discovery, the mechanism of the Grignard reaction remains unresolved. Ambiguities arise from the concomitant presence of multiple organomagnesium species and the competing mechanisms involving either nucleophilic addition or the formation of radical intermediates. To shed light on this topic, quantum-chemical calculations and ab initio molecular dynamics simulations are used to study the reaction of CH3MgCl in tetrahydrofuran with acetaldehyde and fluorenone as prototypical reagents. All organomagnesium species coexisting in solution due to the Schlenk equilibrium are found to be competent reagents for the nucleophilic pathway. The range of activation energies displayed by all of these compounds is relatively small. The most reactive species are a dinuclear Mg complex in which the substrate and the nucleophile initially bind to different Mg centers and the mononuclear dimethyl magnesium. The radical reaction, which requires the homolytic cleavage of the Mg-CH3 bond, cannot occur unless a substrate with a low-lying π*(CO) orbital coordinates the Mg center. This rationalizes why a radical mechanism is detected only in the presence of substrates with a low reduction potential. This feature, in turn, does not necessarily favor the nucleophilic addition, as shown for the reaction with fluorenone. The solvent needs to be considered as a reactant for both the nucleophilic and the radical reactions, and its dynamics is essential for representing the energy profile. The similar reactivity of several species in fast equilibrium implies that the reaction does not occur via a single process but by an ensemble of parallel reactions.

We describe the mechanism of self-aggregation of α-tocopherol transfer protein into a spherical nano-cage employing by Monte Carlo simulations. The protein is modelled by a patchy coarse-grained representation, where the protein-protein interfaces, determined in the past by x-ray diffraction, are represented by simplified two-body interaction potentials. Our results show that the oligomerization kinetics proceeds in two steps, with the formation of meta-stable trimeric units, and the subsequent assembly into the spherical aggregates. Data are in agreement with experimental observations regarding the prevalence of different aggregation states at specific ambient conditions. Finally, our results indicate a route for the experimental stabilization of the trimer, crucial for the understanding of the physiological role of such aggregates in vitamin E body trafficking.

The Schlenk equilibrium is a complex reaction governing the presence of multiple chemical species in solution of Grignard reagents. The full characterization at the molecular level of the transformation of CH3MgCl into MgCl2 and Mg(CH3)2 in tetrahydrofuran (THF) by means of ab initio molecular dynamics simulations with enhanced-sampling metadynamics is presented. The reaction occurs via formation of dinuclear species bridged by chlorine atoms. At room temperature, the different chemical species involved in the reaction accept multiple solvation structures, with two to four THF molecules that can coordinate the Mg atoms. The energy difference between all dinuclear solvated structures is lower than 5 kcal mol-1. The solvent is shown to be a direct key player driving the Schlenk mechanism. In particular, this study illustrates how the most stable symmetrically solvated dinuclear species, (THF)CH3Mg(μ-Cl)2MgCH3(THF) and (THF)CH3Mg(μ-Cl)(μ-CH3)MgCl(THF), need to evolve to less stable asymmetrically solvated species, (THF)CH3Mg(μ-Cl)2MgCH3(THF)2 and (THF)CH3Mg(μ-Cl)(μ-CH3)MgCl(THF)2, in order to yield ligand exchange or product dissociation. In addition, the transferred ligands are always departing from an axial position of a pentacoordinated Mg atom. Thus, solvent dynamics is key to successive Mg-Cl and Mg-CH3 bond cleavages because bond breaking occurs at the most solvated Mg atom and the formation of bonds takes place at the least solvated one. The dynamics of the solvent also contributes to keep relatively flat the free energy profile of the Schlenk equilibrium. These results shed light on one of the most used organometallic reagents whose structure in solvent remains experimentally unresolved. These results may also help to develop a more efficient catalyst for reactions involving these species.