Expanding horizons in conceptual density functional theory: Novel ensembles and descriptors to decipher reactivity patterns.

Conceptual density functional theory (CDFT) and the quantum reactivity descriptors stemming from it have proven to be valuable tools for understanding the chemical behavior of molecules. This article is presented as being intrinsically of dual character. In a first part, it briefly reviews, in a deliberately didactical way, the main ensembles in CDFT, while the second half presents two additional ensembles, where the chemical hardness acts as a natural variable, and their respective reactivity descriptors. The evaluation of these reactivity descriptors on common organic chemical reagents are presented and discussed.

Predicting the activity of methoxyphenol derivatives antioxidants: II-Importance of the nature of the solvent on the mechanism, a DFT study.

The various mechanisms of primary antioxidant action of a series of 2-Methoxyphenols are investigated in the present work. The electronic properties have just been studied in a joint article, so that we focus here on Hydrogen Atom Transfer (HAT), Single Electron Transfer-Proton Transfer (SET-PT) and Sequential Proton-Loss Electron-Transfer (SPLET) mechanisms. The two computational approaches used in the previous study of the structure and reactivity of these compounds [Computational and Theoretical Chemistry 1229 (2023) 114287] provide identical mechanisms trends in gas phase. In gas media, hydrogen atom transfer (HAT) is more favored. On the contrary, the solvent polarity has a significant effect on the mechanism of the antioxidant activity, since the polarity of the solvent increases the contribution of the SPLET mechanism.

Determining the Factors accounting for Reaction Selectivity A Relative Energy Gradient - Interacting Quantum Atoms and Natural Bonding Orbitals Study.

Identifying the main physicochemical properties accounting for the course of a reaction is of utmost importance to rationalize chemical syntheses. To this aim, the relative energy gradient (REG) method is an appealing approach because it is an unbiased and automatic process to extract the most relevant pieces of energy information. Initially formulated within the interacting quantum atoms (IQA) framework for a single reaction, here we extend the REG method to natural bond orbitals (NBO) analysis and to the case of two competitive processes. This development enables the determination of the driving forces of any chemical selectivity. We illustrate the extended REG method on the case study of ring opening in cyclobutenes, which is an important instance of the so-called torquoselectivity.

Understanding the mechanism and regio- and stereo selectivity of [3 + 2] cycloaddition reactions between substituted azomethine ylide and 3,3,3-trifluoro-1-nitroprop-1-ene, within the molecular electron density theory.

The selectivity and the nature of the molecular mechanism of the [3 + 2] cycloaddition (32CA) reaction between 2-(dimethylamino)-1H-indene-1,3(2H)-dione (AY11) and trans(E)-3,3,3-trifluoro-1-nitroprop-1-ene(FNP10) has been studied, in which the molecular electron density theory using density functional theory methods at the MPWB1K/6-31G(d) computational level was used. Analysis of the global reactivity indices permits us to characterize FNP10 as a strong electrophile and AY11 as a strong nucleophile. Four reactive pathways associated with the ortho/meta regioselective channels and endo/exo stereoselective approaches modes have been explored and characterized in the gas phase and in the benzene solvent. The analysis of the relative energies associated with the different reaction pathways indicates that the 32CA reactions of the azomethine ylide (AY) with the nitroalkene (FNP) is meta regioselective with high endo stereoselectivity. This result is in good agreement with the experimental observations. electron localization function topological analysis of the most favored reactive pathways allows for characterizing the mechanism of this 32CA reactions as a non-concerted two-stage one-step mechanism. Finally, non-covalent interactions and quantum theory of atoms in molecule analyses at the meta/endo transition state structure indicate that the presence of different several weak interactions, namely, CF and NH contributed in favoring the formation of a meta-endo cycloadduct.

Insight into the Varying Reactivity of Different Catalysts for CO2 Cycloaddition into Styrene Oxide: An Experimental and DFT Study.

The cycloaddition of CO2 into epoxides to form cyclic carbonates is a highly sought-after reaction for its potential to both reduce and use CO2, which is a greenhouse gas. In this paper, we present experimental and theoretical studies and a mechanistic approach for three catalytic systems. First, as Lewis base catalysts, imidazole and its derivatives, then as a Lewis acid catalyst, ZnI2 alone, and after that, the combined system of ZnI2 and imidazole. In the former, we aimed to discover the reasons for the varied reactivities of five Lewis base catalysts. Furthermore, we succeeded in reproducing the experimental results and trends using DFT. To add, we emphasized the importance of non-covalent interactions and their role in reactivity. In our case, the presence of a hydrogen bond was a key factor in decreasing the reactivity of some catalysts, thus leading to lower conversion rates. Finally, mechanistically understanding this 100% atom economy reaction can aid experimental chemists in designing better and more efficient catalytic systems.

Reactivity and a Charge-Transfer Model Analysis in Aminopolycarboxylic-Metal Complexes.

In the present work, we have calculated several density functional theory (DFT) reactivity descriptors for the aminopolycarboxylate (APC) acids at the B3LYP/6311++G (d,p) levels of theory, aiming to analyze their reactivity. Reactivity descriptors such as ionization energy, molecular hardness, electrophilicity, and condensed Fukui function local indices have been determined to predict the reactivity of APCs. The influence of the solvent was taken into account by employing the CPCM model. The results indicate that the solvation slightly modifies the tendency of the reactivity of the APCs studied. On the other hand, we applied a global and local charge-transfer partitioning model, which introduces two charge-transfer channels [one for accepting electrons (electrophilic) and another for donating one (nucleophilic)] to the complexation reaction of a set of APC acids with transition metals (Mn, Co, and Ni targets enlarged by Fe, Cu, and Zn). The correlation between the charges obtained for the interaction between APC acids and transition metal stability constants provides support for their interpretation as measures of the electrophilicity and nucleophilicity of a chemical species and, at the same time, allows one to describe the donation and back-donation processes in terms of the DFT of chemical reactivity. Also, the application of dual descriptors for these acids provides valuable information concerning the atoms in the reactants playing the most important roles in the reaction, thus helping to improve our understanding of the reaction under study.

Formation Mechanism of Inter-Crosslink in DNA by Nitrogen Oxides Pollutants through A Diazonium Intermediate.

Outdoor air pollution is a mixture of multiple atmospheric pollutants, among which nitrogen oxide (NOx) stands out due to its association with several diseases. NOx reactivity can conduct to DNA damage as severe as interstrand crosslinks (ICL) formation, that in turn is able to block DNA replication and transcription. Experimental studies have suggested that the ICL formation due to NOx is realized through a diazonium intermediate (DI). In this work, we have modeled the DI structure, including a DNA double-strand composed of two base pairs GC/CG, being diazotized as one of the guanine nucleotides. The structural stability of DNA with DI lesion was essayed through 500 ns molecular dynamics simulations. It was found that the DNA structure of the oligonucleotide is stable when the DI is present since the loss of a Guanine-Cytosine hydrogen bond is replaced by the presence of two cation-π interactions. Additionally, we have studied the mechanism of formation of a crosslink between the two guanine nucleobases from the modeled DI by carrying out DFT calculations at the M06-L/DNP+ level of theory. Our results show that the mechanism is thermodynamically favored by a strong stabilization of the ICL product, and the process is kinetically viable since its limiting stage is accessible.

Understanding the intermolecular Diels-Alder cycloaddition promotion: Activation strain model/energy decomposition analysis model and conceptual density functional theory viewpoints.

The present work reports the computational study of the major Diels-Alder reaction between 2-bromocycloalkenone and a variety of mono- and di-substituted dienes. Through density functional theory (DFT) calculations and subsequent activation strain model/energy decomposition analysis/conceptual DFT (C-DFT) analyses, the key factors governing the activation barriers heights, and thus reactivity, are characterized. In contrast with a previous study, steric effects do not appear to control reactivity. Conversely, in all presented cases, a subtle interplay between deformation and interaction energies is evidenced at transition states. In the end, neither term alone is enough to explain or predict reactivity. Yet a simple C-DFT descriptor allows to predict with a reasonable efficiency the activation barriers: the excitation energy needed to observe a charge transfer from the diene to the dienophile. Theoretical elements are provided to support the use of this descriptor.

The density polarization reveals directions of electron displacements due to the substituent effect: Analysis performed on a metal-organic Mo-Oxo catalyst.

Some Mo-oxo complexes bearing pyridine rings have the capability for dihydrogen production from water. However, energy barrier and overall energy vary depending on the effect exerted by several substituent groups located at different positions around one or more pyridine rings which are ligands of these compounds. Based on the Karunadasa and coworkers investigation where the para-position was experimentally tested in compounds derivatised from the 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine oxo-molybdenum complex synthesized (Karunadasa et al., Nature, 2010, 464, 1329), we tested the combined effect of electron-withdrawing and electron-donating groups simulated as perturbations represented by point-charges. Then, we used the density polarization concept, δρ(r), a local reactivity descriptor corresponding to the partially integrated linear response function, χ(r, r') (a non-local reactivity descriptor), which is able to reveal different displacements of π-electrons on molecular structures. We perturbed the para-positions in the pentadentate ligand 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine in the Mo-based complex by means of point-charges. They were located in three different configurations of the organic ligand (trans, geminal, and cis) which could help to explain energy barriers and overall energy of reactions catalyzed by this type of Mo-complexes. Our results indicate that the trans configuration of point-charges induces the most amount of fraction of electron shifted on the complex. A Mo-based complex bearing the same trans configuration for electron-withdrawing and electron-donating substituent groups (cyano and amino, respectively), leads to a kinetically more favorable H2 release than the cis or geminal configuration of the substituent groups aforementioned.

A computational investigation of the selectivity and mechanism of the Lewis acid catalyzed oxa-Diels-Alder cycloaddition of substituted diene with benzaldehyde.

The selectivity and the mechanism of the uncatalyzed and AlCl3 catalyzed hetero-Diels-Alder reaction (HDR) between ([E]-4-methylpenta-2,4-dienyloxy)(tert-butyl)dimethylsilane 1 and benzaldehyde 2 have been studied using density functional theory at the MPWB1K/6-31G(d) level of theory. The uncatalyzed HDR between diene 1 and alkene 2 is characterized by a polar character and proceeds via an asynchronous one-step mechanism for the meta paths and synchronous for the ortho ones. In the presence of AlCl3 catalyst, the mechanism changes to be stepwise, while the first step is the rate-determining step. The activation energies widely decrease, and the polar character increases dramatically. A large analysis of the mechanism is performed using the activation strain model/energy decomposition analysis (ASM/EDA) model, the natural bond orbital (NBO) and state specific dual descriptors (SSDDs). The obtained results indicate that the combined interaction energy associated with the distortion of the reactants in these HDR are at the origin of the observed kinetics. NBO analyses were applied to estimate the Lewis-acid catalyst donor-acceptor interaction with the molecular system. The SSDD analysis shed light into the orientation effects on the reaction kinetics by providing important information about charge transfer interactions during the chemical reaction. It indicates that the more favorable HDR pathway have the lowest excitation energies, facilitating the interaction between diene 1 and benzaldehyde 2 moieties. Non-covalent interaction (NCI) and QTAIM analyses of the meta-endo structure indicate that the presence of several weak NCIs formed at this approach is at the origin of the meta-endo selectivity.

Polarisation of Electron Density and Electronic Effects: Revisiting the Carbon-Halogen Bonds.

Electronic effects (inductive and mesomeric) are of fundamental importance to understand the reactivity and selectivity of a molecule. In this article, polarisation temperature is used as a principal index to describe how electronic effects propagate in halogeno-alkanes and halogeno-alkenes. It is found that as chain length increases, polarisation temperature decreases. As expected, polarisation is much larger for alkenes than for alkanes. Finally, the polarisation mode of the carbon-fluorine bond is found to be quite different and might explain the unusual reactivity of fluoride compounds.

A statistical thermodynamics view of electron density polarisation: application to chemical selectivity.

A fundamental link between conceptual density functional theory and statistical thermodynamics is herein drawn, showing that intermolecular electrostatic interactions can be understood in terms of effective work and heat exchange. From a more detailed analysis of the heat exchange in a perturbation theory framework, an associated entropy can be subsequently derived, which appears to be a suitable descriptor for the local polarisability of the electron density. A general rule of thumb is evidenced: the more the perturbation can be spread, both through space and among the excited states, the larger the heat exchange and entropy.

Substituent Effect on the Himbert Intramolecular Arene/Allene Diels-Alder Reaction: NBO Analysis and State Specific Dual Descriptors.

Conceptual density functional theory has been applied to study the Himbert intramolecular arene/allene Diels-Alder reaction. The effect of substitutions at different positions on the kinetics of these reactions has been analyzed. Therefore, from the calculation of the activation energies of more than 27 reactions involving concerted mechanisms, the selectivity of these reactions can be predicted and rationalized with the aid of conceptual DFT descriptors. An application of the two concepts, natural population analysis (NBO) and the state-specific dual descriptor (SSDD) for evaluating substituent effects, allows the investigation of the different interactions that promote a reaction compared to another. The SSDDs computed for the transition state structures provide important information about charge transfer interactions during the chemical reaction. In our case, the SSDD results show that the substituents promoting Himbert reaction have the lowest excitation energies, a fact which facilitates the allene/arene interaction. The NBO results show that according to the nature of the substituent, the Himbert reaction stands as a normal-electron demand or reverse. Thus, the interactions favoring each reaction are mentioned. The geometric deformation observed in the case of OCH3 is at the origin to the emergence of other low interactions between diene and dienophile as well as a strong electronic delocalization stabilizing the arene moiety. The calculated synchronicity indexes show that the Himbert intramolecular Diels-Alder reactions are very synchronous.

Interstrand cross-linking implies contrasting structural consequences for DNA: insights from molecular dynamics.

Oxidatively-generated interstrand cross-links rank among the most deleterious DNA lesions. They originate from abasic sites, whose aldehyde group can form a covalent adduct after condensation with the exocyclic amino group of purines, sometimes with remarkably high yields. We use explicit solvent molecular dynamics simulations to unravel the structures and mechanical properties of two DNA sequences containing an interstrand cross-link. Our simulations palliate the absence of experimental structural and stiffness information for such DNA lesions and provide an unprecedented insight into the DNA embedding of lesions that represent a major challenge for DNA replication, transcription and gene regulation by preventing strand separation. Our results based on quantum chemical calculations also suggest that the embedding of the ICL within the duplex can tune the reaction profile, and hence can be responsible for the high difference in yields of formation.

Conformational polymorphism or structural invariance in DNA photoinduced lesions: implications for repair rates.

DNA photolesions constitute a particularly deleterious class of molecular defects responsible for the insurgence of a vast majority of skin malignant tumors. Dimerization of two adjacent thymines or cytosines mostly gives rise to cyclobutane pyrimidine dimers (CPD) and pyrimidine(6-4)pyrimidone 64-PP as the most common defects. We perform all-atom classical simulations, up to 2 µs, of CPD and 64-PP embedded in a 16-bp duplex, which reveal the constrasted behavior of the two lesions. In particular we evidence a very limited structural deformation induced by CPD while 64-PP is characterized by a complex structural polymorphism. Our simulations also allow to unify the contrasting experimental structural results obtained by nuclear magnetic resonance or Förster Resonant Energy Transfer method, showing that both low and high bent structures are indeed accessible. These contrasting behaviors can also explain repair resistance or the different replication obstruction, and hence the genotoxicity of these two photolesions.

Unexpected Structure of a Helical N4 -Schiff-Base Zn(II) Complex and Its Demetallation: Experimental and Theoretical Studies.

A new Zn-N4 -Schiff base L=((±)-trans-N,N'-Bis(2-aminobenzylidene)-1,2-diaminocyclohexane) complex was synthesized and fully characterized, showing an unexpected self-assembled double-stranded helicate structure. The X-ray crystal analysis of the Zn2 L2 complex ((C40 H44 N8 Zn2 ,CH2 Cl2 , a=14.2375(3) Å, b=16.7976(4) Å, c=16.1613(4) Å, monoclinic, P21 /n, Z=4) shows a centrosymmetrical structure in which zinc atoms are in distorted tetrahedral environments, revealing an M- (R, R) left-handed helicity in its asymmetric unit. However, it was observed that this dinuclear complex is thermodynamically unstable in the presence of small water amounts and undergoes demetallation into free N4-Schiff base ligand and ZnO nanoparticles. This hydrolysis process was thoroughly identified and monitored through detailed 1 H NMR, DOSY NMR analysis. The reaction mechanism of this demetallation event was elucidated by using the DFT method, involving an activation energy smaller than 13 kcal/mol. Besides, a theoretical mechanism of the demetallation process is given for the first time.

Fluorine substituent effect on the stereochemistry of catalyzed and non-catalyzed Diels-Alder reactions. The case of R-butenone with cyclopentadiene: a computational assessment of the mechanism.

The present work studies theoretically the mechanisms involved in the fluorine substituent effect on the stereochemistry of Diels-Alder reactions. The case of R-butenone with cyclopentadiene is used for the purpose of modelling more general α-fluoro-α,ß-unsaturated carbonyl compounds, in catalyzed and uncatalyzed cases. A thorough analysis of the mechanism is performed using energy decomposition analysis (EDA) and conceptual DFT tools. It is shown that the endo conformation is privileged in all the studied cases with the exception of the α-fluorinated ketone. It is found that the endo selectivity of the non-fluorinated reactions is only due to the decrease of dispersion energy. On the other hand, the presence of a fluorine atom in the dienophile moieties increases remarkably the magnitude not only of the interaction energy between the reactants but that of the strain energy as well. Moreover, it is the strong destabilization strain energy occurring at the transition state of the endo pathway of the reaction cyclopentadiene/3-fluorobutenone that is mainly responsible for the exo selectivity. The effect of a Lewis acid catalyst on these reactions is also studied. The Lewis acid catalyst affects the activation energy of the studied Diels-Alder reactions but not their stereoselectivity. Furthermore, the dual descriptor results shed light onto the mechanism. Besides, natural bond orbital analysis (NBO) and determination of the condensed values of the state-specific dual descriptors (SSDD) are carried out to evaluate the donor-acceptor properties in these reactions. For the first time, a semiquantitative prediction of stereoselectivity due to substitutions of dienophile is obtained, thus complementing the previous interpretations (R. Hoffmann and R. B. Woodward, J. Am. Chem. Soc., 1965, 87, 4388; I. Fernández and F. M. Bickelhaupt, Chem. Soc. Rev., 2014, 43, 4953). Finally, since the role of dispersion forces is evidenced in some cases, a comparison between some popular exchange-correlation functionals is presented, assessing the performance of some standard functionals besides functionals with explicit dispersion corrections.

Probing the reactivity of singlet oxygen with purines.

The reaction of singlet molecular oxygen with purine DNA bases is investigated by computational means. We support the formation of a transient endoperoxide for guanine and by classical molecular dynamics simulations we demonstrate that the formation of this adduct does not affect the B-helicity. We thus identify the guanine endoperoxide as a key intermediate, confirming a low-temperature nuclear magnetic resonance proof of its existence, and we delineate its degradation pathway, tracing back the preferential formation of 8-oxoguanine versus spiro-derivates in B-DNA. Finally, the latter oxidized 8-oxodGuo product exhibits an almost barrierless reaction profile, and hence is found, coherently with experience, to be much more reactive than guanine itself. On the contrary, in agreement with experimental observations, singlet-oxygen reactivity onto adenine is kinetically blocked by a higher energy transition state.

Correlation of bistranded clustered abasic DNA lesion processing with structural and dynamic DNA helix distortion.

Clustered apurinic/apyrimidinic (AP; abasic) DNA lesions produced by ionizing radiation are by far more cytotoxic than isolated AP lesion entities. The structure and dynamics of a series of seven 23-bp oligonucleotides featuring simple bistranded clustered damage sites, comprising of two AP sites, zero, one, three or five bases 3' or 5' apart from each other, were investigated through 400 ns explicit solvent molecular dynamics simulations. They provide representative structures of synthetically engineered multiply damage sites-containing oligonucleotides whose repair was investigated experimentally (Nucl. Acids Res. 2004, 32:5609-5620; Nucl. Acids Res. 2002, 30: 2800-2808). The inspection of extrahelical positioning of the AP sites, bulge and non Watson-Crick hydrogen bonding corroborates the experimental measurements of repair efficiencies by bacterial or human AP endonucleases Nfo and APE1, respectively. This study provides unprecedented knowledge into the structure and dynamics of clustered abasic DNA lesions, notably rationalizing the non-symmetry with respect to 3' to 5' position. In addition, it provides strong mechanistic insights and basis for future studies on the effects of clustered DNA damage on the recognition and processing of these lesions by bacterial or human DNA repair enzymes specialized in the processing of such lesions.

Molecular Dynamics Insights into Polyamine-DNA Binding Modes: Implications for Cross-Link Selectivity.

Biogenic polyamines, which play a role in DNA condensation and stabilization, are ubiquitous and are found at millimolar concentration in the nucleus of eukaryotic cells. The interaction modes of three polyamines-putrescine (Put), spermine (Spm), and spermidine (Spd)-with a self-complementary 16 base pair (bp) duplex, are investigated by all-atom explicit-solvent molecular dynamics. The length of the amine aliphatic chain leads to a change of the interaction mode from minor groove binding to major groove binding. Through all-atom dynamics, noncovalent interactions that stabilize the polyamine-DNA complex and prefigure the reactivity, leading to the low-barrier formation of deleterious DNA-polyamine cross-links, after one-electron oxidation of a guanine nucleobase, are unraveled. The binding strength is quantified from the obtained trajectories by molecular mechanics generalized Born surface area post-processing (MM-GBSA). The values of binding free energies provide the same affinity order, Put