Reaction Electronic Flux Perspective on the Mechanism of the Zimmerman Di-π-methane Rearrangement

How does dopamine convert into norepinephrine? Insights on the key step of the reaction

CONTEXT

Dopamine β -monooxygenase (D β M) is an essential enzyme in the organism that regioselectively converts dopamine into R-norepinephrine, the key step of the reaction, studied in this paper, is a hydrogen atom transfer (HAT) from dopamine to a superoxo complex on D β M, forming a hydroperoxo intermediate and dopamine radical. It was found that the formation of a hydrogen bond between dopamine and the D β M catalyst strengthens the substrate-enzyme interaction and facilitates the HAT which takes place selectively to give the desired enantiomeric form of the product. Six reactions leading to the hydroperoxo intermediate were analyzed in detail using theoretical and computational tools in order to identify the most probable reaction mechanism. The reaction force analysis has been used to demonstrate that the nature of the activation energy is mostly structural and largely due to the initial approach of species in order to get closer to each other to facilitate the hydrogen abstraction. On the other hand, the reaction electronic flux revealed that electronic activity driving the reactions is triggered by polarization effects and, in the most probable reaction among the six studied, it takes place in a concerted and non-spontaneous way. Chemical events driving the reaction have been identified and the energy absorbed or delivered by each one was quantified in detail.

METHODS

The dopamine and a computational model of the copper superoxo complex on D β M were optimized at B3LYP-D3(BJ)/6-311 G(d,p) level theory in the Gaussian 16 software package. Optimization and IRC calculations were performed in the gas phase and through the PCM solvation model to mimic the protein medium. Non-covalent interactions were plotted using the NCI-plot software.

Towards a complete description of the reaction mechanisms between nitrenium ions and water

CONTEXT

Nitrenium ions are intermediates in the metabolic routes producing the highly carcinogenic nitrosamines and binding to DNA molecules. The reaction mechanism of nitrenium molecules with explicit water molecules is sensibly dependent on the number of waters: when a second molecule is involved, it acts as a catalyst for the reaction, lowering intrinsic activation barriers regardless of the substituent. For all cases, the reaction force constants and reaction electron flux indicate highly synchronous reactions for n = 1 . Conversely, for n = 2 highly non-synchronous reactions are obtained, involving two separate proton transfers happening early and late in the reaction path. As a test case, for the simplest[ NH 2 ] + + 2 H 2 O reactions, orbital interactions within the NBO paradigm, bond orders, and their derivatives indicate that each individual proton transfer is highly synchronous.

METHODS

Molecular geometries were optimized and characterized at the B3LYP/6-311++G(d, p) level. Intrinsic reaction coordinates were calculated. CCSD(T) single point energies with the same basis were computed on all stationary points. The reaction force, reaction force constant, and reaction electron flux are used to study the evolution of the reacting systems. Natural bond orbitals are used to understand the primitive changes driving the reaction.

The transition state region in nonsynchronous concerted reactions

The critical and vanishing points of the reaction force F(ξ) = -dV(ξ)/dξ yield five important coordinates (ξ_R, ξ_R^* , ξ_TS, ξ_P^* , ξ_P) along the intrinsic reaction coordinate (IRC) for a given concerted reaction or reaction step. These points partition the IRC into three well-defined regions, reactants (ξ_R→ξ_R^* ), transition state (ξ_R^* →ξ_P^* ), and products (ξ_P^* →ξ_P), with traditional roles of mostly structural changes associated with the reactants and products regions and mostly electronic activity associated with the transition state (TS) region. Following the evolution of chemical bonding along the IRC using formal descriptors of synchronicity, reaction electron flux, Wiberg bond orders, and their derivatives (or, more precisely, the intensity of the electron activity) unambiguously indicates that for nonsynchronous reactions, electron activity transcends the TS region and takes place well into the reactants and products regions. Under these circumstances, an extension of the TS region toward the reactants and products regions may occur.

How Does Electronic Activity Drive Chemical Reactions? Insights from the Reaction Electronic Flux for the Conversion of Dopamine into Norepinephrine

Hydrogen atom transfer (HAT) is a crucial step in the physiological conversion of dopamine into norepinephrine catalyzed by dopamine β-monooxygenase. The way the reaction takes place is unclear, and a rational explanation on how the electronic activity drives the HAT seems to be necessary. In this work, we answer this question using the reaction electronic flux (REF), a DFT-based descriptor of electronic activity. Two reaction mechanisms will be analyzed using the REF's decomposition in polarization and electron transfer effects. Results show that both mechanisms proceed as follows: (1) polarization effects initiate the reactions producing structural distortions; (2) electron transfer processes take over near the transition states, triggering specific chemical events such as bond forming and breaking which are responsible to push the reactions toward the products; (3) after passing the transition state, polarization shows up again and drives the relaxation process toward the product. Similar polarization effects were observed in both reactions, but they present an opposite behavior of the electronic transfer flux disclosing the fact that electron transfer phenomena govern the reaction mechanisms.

Symmetry-Adapted Perturbation Theory Decomposition of the Reaction Force: Insights into Substituent Effects Involved in Hemiacetal Formation Mechanisms

The decomposition of the reaction force based on symmetry-adapted perturbation theory (SAPT) has been proposed. This approach was used to investigate the substituent effects along the reaction coordinate pathway for the hemiacetal formation mechanism between methanol and substituted aldehydes of the form CX₃CHO (X = H, F, Cl, and Br), providing a quantitative evaluation of the reaction-driving and reaction-retarding force components. Our results highlight the importance of more favorable electrostatic and induction effects in the reactions involving halogenated aldehydes that leads to lower activation energy barriers. These substituent effects are further elucidated by applying the functional-group partition of symmetry-adapted perturbation theory (F-SAPT). The results show that the reaction is largely driven by favorable direct noncovalent interactions between the CX₃ group on the aldehyde and the OH group on methanol.

Theoretical insights into the E1cB/E2 mechanistic dichotomy of elimination reactions

Modeling the E1cB, E2, and E1cB/E2 borderline mechanisms in terms of carbanion stabilization.