1
|
Forero-Girón AC, Toro-Labbé A. How does dopamine convert into norepinephrine? Insights on the key step of the reaction. J Mol Model 2025; 31:32. [PMID: 39751655 DOI: 10.1007/s00894-024-06256-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2024] [Accepted: 12/14/2024] [Indexed: 01/04/2025]
Abstract
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.
Collapse
Affiliation(s)
- Angie Carolay Forero-Girón
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avenida Vicuña Mackenna 4860, 7820436, Santiago de Chile, Chile.
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avenida Vicuña Mackenna 4860, 7820436, Santiago de Chile, Chile
| |
Collapse
|
2
|
Balhara R, Chatterjee R, Jindal G. Mechanism and stereoselectivity in metal and enzyme catalyzed carbene insertion into X-H and C(sp 2)-H bonds. Chem Soc Rev 2024; 53:11004-11044. [PMID: 39392229 DOI: 10.1039/d4cs00742e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/12/2024]
Abstract
Constructing highly proficient C-X (X = O, N, S, etc.) and C-C bonds by leveraging TMs (transition metals) (Fe, Cu, Pd, Rh, Au, etc.) and enzymes to catalyze carbene insertion into X-H/C(sp2)-H is a highly versatile strategy. This is primarily achieved through the in situ generation of metal carbenes from the interaction of TMs with diazo compounds. Over the last few decades, significant advancements have been made, encompassing a wide array of X-H bond insertions using various TMs. These reactions typically favor a stepwise ionic pathway where the nucleophilic attack on the metal carbene leads to the generation of a metal ylide species. This intermediate marks a critical juncture in the reaction cascade, presenting multiple avenues for proton transfer to yield the X-H inserted product. The mechanism of C(sp2)-H insertion reactions closely resembles those of X-H insertion reactions and thus have been included here. A major development in carbene insertion reactions has been the use of engineered enzymes as catalysts. Since the seminal report of a non-natural "carbene transferase" by Arnold in 2013, "P411", several heme-based enzymes have been reported in the literature to catalyze various abiological carbene insertion reactions into C(sp2)-H, N-H and S-H bonds. These enzymes possess an extraordinary ability to regulate the orientation and conformations of reactive intermediates, facilitating stereoselective carbene transfers. However, the absence of a suitable stereochemical model has impeded the development of asymmetric reactions employing a lone chiral catalyst, including enzymes. There is a pressing need to investigate alternative mechanisms and models to enhance our comprehension of stereoselectivity in these processes, which will be crucial for advancing the fields of asymmetric synthesis and biocatalysis. The current review aims to provide details on the mechanistic aspects of the asymmetric X-H and C(sp2)-H insertion reactions catalyzed by Fe, Cu, Pd, Rh, Au, and enzymes, focusing on the detailed mechanism and stereochemical model. The review is divided into sections focusing on a specific X-H/C(sp2)-H bond type catalyzed by different TMs and enzymes.
Collapse
Affiliation(s)
- Reena Balhara
- Department of Organic Chemistry, Indian Institute of Science, Bengaluru-560012, Karnataka, India.
| | - Ritwika Chatterjee
- Department of Organic Chemistry, Indian Institute of Science, Bengaluru-560012, Karnataka, India.
| | - Garima Jindal
- Department of Organic Chemistry, Indian Institute of Science, Bengaluru-560012, Karnataka, India.
| |
Collapse
|
3
|
Barrales-Martínez C, Durán R, Caballero J. Shannon entropy variation as a global indicator of electron density contraction at interatomic regions in chemical reactions. J Mol Model 2024; 30:371. [PMID: 39382590 DOI: 10.1007/s00894-024-06171-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2024] [Accepted: 10/04/2024] [Indexed: 10/10/2024]
Abstract
CONTEXT The negative of the Shannon entropy derivative is proposed to account for electron density contraction as the chemical bonds are breaking and forming during a chemical reaction. We called this property the electron density contraction index, EDC, which allows identifying stages in a reaction that are dominated by electron contraction or expansion. Four different reactions were analyzed to show how the EDC index changes along the reaction coordinate. The results indicate that the rate of change of Shannon entropy is directly related to the rate of change of the electron density at the bond critical points between all the atomic pairs in the molecular systems. It is expected that EDC will complement the detailed analysis of reaction mechanisms that can be performed with the theoretical tools available to date. METHODS Density functional theory calculations at the B3LYP/6-31G(d,p) level of theory were carried out using Gaussian 16 to analyze the reaction mechanisms of the four reactions studied. The reaction paths were obtained via the intrinsic reaction coordinate method, which served as the reaction coordinate to obtain the reaction force and the EDC profiles in each case. Shannon entropy and electron density at the bond critical points were calculated using the Multiwfn 3.7 package.
Collapse
Affiliation(s)
- César Barrales-Martínez
- Instituto de Investigación Interdisciplinaria (I3), Vicerrectoría Académica, Universidad de Talca, Campus Talca, Talca, Chile.
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Talca, Talca, Chile.
| | - Rocío Durán
- Facultad de Ciencias, Departamento de Química Ambiental, Universidad Católica de la Santísima Concepción, Concepción, Chile.
| | - Julio Caballero
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Talca, Talca, Chile
| |
Collapse
|
4
|
Barrales-Martínez C, Durán R, Jaque P. Transition from synchronous to asynchronous mechanisms in 1,3-dipolar cycloadditions: a polarizability perspective. J Mol Model 2024; 30:355. [PMID: 39347808 DOI: 10.1007/s00894-024-06161-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2024] [Accepted: 09/23/2024] [Indexed: 10/01/2024]
Abstract
CONTEXT This study investigates the energetic and polarizability characteristics of three 1,3-dipolar cycloaddition reactions between diazene oxide and substituted ethylenes, focusing on the transition from synchronous to asynchronous mechanisms. Synchronicity analysis, using the reaction force constant, indicates that the bond evolution process becomes increasingly decoupled as the number of cyano groups increases. Polarizability analysis reveals that isotropic polarizability reaches its maximum near the transition state in all cases, while anisotropy of polarizability shifts from the transition state toward the product direction as asynchronicity increases. The larger the shift, the more asynchronous the mechanism, as reflected by the weight of the transition region. A detailed examination of the parallel and perpendicular polarizability components to the newly formed sigma bonds shows that the evolution of the parallel component is closely aligned with the energetic changes along the reaction coordinate, particularly in the synchronous reaction. We have also identified a relationship between the displacement in the maximum state of the parallel component from the transition state and the synchronicity of the mechanism. The larger the displacement, the more asynchronous the mechanism. These findings suggest that asynchronous 1,3-dipolar cycloaddition mechanisms are characterized by a decoupling of isotropic and anisotropic polarizabilities and a shift in the maximum polarizability state of the parallel component toward the product direction. METHODS Density functional theory calculations were performed at the B3LYP/6-311 + + G(d,p)//B3LYP/6-31G(d,p) level of theory. The polarizability was calculated at each point of the reaction path, obtained using the intrinsic reaction coordinate method, as implemented in Gaussian 16.
Collapse
Affiliation(s)
- César Barrales-Martínez
- Instituto de Investigación Interdisciplinaria (I3), Vicerrectoría Académica, Universidad de Talca, Campus Talca, Talca, Chile.
- Facultad de Ingeniería, Centro de Bioinformática, Simulación y Modelado (CBSM), Universidad de Talca, Campus Talca, Talca, Chile.
| | - Rocío Durán
- Departamento de Química Ambiental, Facultad de Ciencias, Universidad Católica de La Santísima Concepción, Concepción, Chile.
| | - Pablo Jaque
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
- Centro de Modelamiento Molecular, Biofísica y Bioinformática, CM2B2, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
| |
Collapse
|
5
|
Durán R, Barrales-Martínez C, Matute RA. Hidden intermediate activation: a concept to elucidate the reaction mechanism of the Schmittel cyclization of enyne-allenes. Phys Chem Chem Phys 2023; 25:6050-6059. [PMID: 36458512 DOI: 10.1039/d2cp04684a] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
The mechanistic paradigm in which the Schmittel cyclization transitions from one-step to stepwise has been investigated through the stabilization of a full hidden intermediate in the framework of the Diabatic Model of Intermediate Stabilization. Hidden intermediate activation was studied in silico employing quasi-classical trajectories and the Electron Localization Function. The stabilization of hidden intermediates achieved by substituting enyne-allenes with cyano and nitro groups generates the appearance of a partially hidden and an explicit intermediate, leading to one-step asynchronous biradical and stepwise biradical/zwitterionic mechanisms, respectively. The mechanistic feature associated with the activation level of the hidden intermediate arises from the Thornton effect and non-RRKM dynamics, where in the case of the CN-substituted system, despite having a single transition state, 54% of the effective trajectories remain in the intermediate zone after 540 fs, indicating that a mixture of mechanisms is observed.
Collapse
Affiliation(s)
- Rocío Durán
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Lircay, Talca 3460000, Chile. .,Instituto de Investigación Interdisciplinaria (I3), Universidad de Talca, Campus Lircay, Talca 3460000, Chile
| | - César Barrales-Martínez
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago 8380492, Chile.,Centro Integrativo de Biología y Química Aplicada (CIBQA), Universidad Bernardo O'Higgins, General Gana 1702, Santiago 8370854, Chile.
| | - Ricardo A Matute
- Centro Integrativo de Biología y Química Aplicada (CIBQA), Universidad Bernardo O'Higgins, General Gana 1702, Santiago 8370854, Chile. .,Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
6
|
Forero-Girón AC, Toro-Labbé A. How Does Electronic Activity Drive Chemical Reactions? Insights from the Reaction Electronic Flux for the Conversion of Dopamine into Norepinephrine. J Phys Chem A 2022; 126:4156-4163. [PMID: 35748576 DOI: 10.1021/acs.jpca.2c01469] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Angie Carolay Forero-Girón
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avenida Vicuña Mackenna 4860, Santiago de Chile, 7820436, Chile
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Avenida Vicuña Mackenna 4860, Santiago de Chile, 7820436, Chile
| |
Collapse
|
7
|
Barrales-Martínez C, Jaque P. A deeper analysis of the role of synchronicity on the Bell-Evans-Polanyi plot in multibond chemical reactions: a path-dependent reaction force constant. Phys Chem Chem Phys 2022; 24:14772-14779. [PMID: 35686531 DOI: 10.1039/d2cp01460b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The role of the degree of synchronicity in the formation of the new single-bonds in a large set of 1,3-dipolar cycloadditions and its relation in the fulfilment of the classical Bell-Evans-Polanyi principle and Hammond-Leffler postulate are deeply investigated. Our results confirm that asynchronicity is an important path-dependent factor to be taken into account: (i) the Bell-Evans-Polanyi is fulfilled as the degree of (a)synchronicity is quite the same, and a linear relationship between reorganisation energy and asynchronicity is found; (ii) the asynchronicity is the origin of deviations of this classical principle of chemical reactivity since any decrease of the energy barrier is due to an increase of asynchronicity at the same exothermicity; and (iii) the less exothermic the reaction is, the more asynchronous the mechanism is, at the same energy barrier. Thus, this implies that TS imbalance decreases the reorganisation energy, consequently affecting the reaction exothermicity as well.
Collapse
Affiliation(s)
- César Barrales-Martínez
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile. .,Centro de Modelamiento Molecular, Biofísica y Bioinformática, CM2B2, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
| | - Pablo Jaque
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile. .,Centro de Modelamiento Molecular, Biofísica y Bioinformática, CM2B2, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
| |
Collapse
|
8
|
Durán R, Herrera B. Theoretical study of the substituent effect on the O–H insertion reaction of copper carbenoids. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02876-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
9
|
Barrales-Martínez C, Gutiérrez-Oliva S, Toro-Labbé A, Pendás ÁM. Interacting Quantum Atoms Analysis of the Reaction Force: A Tool to Analyze Driving and Retarding Forces in Chemical Reactions. Chemphyschem 2021; 22:1976-1988. [PMID: 34293240 DOI: 10.1002/cphc.202100428] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Revised: 07/21/2021] [Indexed: 12/22/2022]
Abstract
The analysis of the reaction force and its topology has provided a wide range of fruitful concepts in the theory of chemical reactivity over the years, allowing to identify chemically relevant regions along a reaction profile. The reaction force (RF), a projection of the Hellmann-Feynman forces acting on the nuclei of a molecular system onto a suitable reaction coordinate, is partitioned using the interacting quantum atoms approach (IQA). The exact IQA molecular energy decomposition is now shown to open a unique window to identify and quantify the chemical entities that drive or retard a chemical reaction. The RF/IQA coupling offers an extraordinarily detailed view of the type and number of elementary processes that take reactants into products, as tested on two sets of simple reactions.
Collapse
Affiliation(s)
- César Barrales-Martínez
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
| | - Soledad Gutiérrez-Oliva
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago, Chile
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago, Chile
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, 33006, Oviedo, Spain
| |
Collapse
|